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1 Executive summary 

Recently, neutrino oscillation experiments have unequivocally demonstrated that neu- 
trinos have mass and that the neutrino mass eigenstates mix. These experiments have 
yielded valuable information on the mixing angles and on the mass differences of the 
three eigenstates. They cannot, however, determine the scale of the neutrino mass, 
which is fixed by the lightest neutrino mass eigenvalue. This can only be directly deter- 
mined by beta decay end point spectral shape measurements, or in the case of Majorana 
neutrinos, by the observation and measurement of the neutrinoless double-beta decay 
half-life. The oscillation experiments yield values for the mixing angles and mass dif- 
ferences accurate enough to allow the prediction of a range of values of the effective 
mass of the Majorana electron neutrino. The CUORE experiment is designed with a 
sensitivity capable of probing all but a small portion of this range. 

CUORE is a proposed array of 988, 750 g, Te02 bolometers operating at 8 to 10 
milli-degrees Kelvin, hanging on the bottom of the mixing chamber of a dilution re- 
frigerator. Natural tellurium has a 33.87% abundance of ^'^"^Te , which is an even-even 
nucleus that can only decay by double-beta decay. This large natural abundance elimi- 
nates the requirement for the very expensive isotopic enrichment required in all of the 
other proposed next generation experiments. The signal of neutrinoless double beta 
decay would be a sharp peak at 2528.8 keV in the bolometric energy spectrum. The 
proposed array is a cylindrical geometry of 19 towers of 52 detectors each. Each tower 
consists of thirteen, four-detector modules. One such tower has been successfully con- 
structed and is being operated in the Gran Sasso Laboratory as a test experiment, and 
also as an independent, first-step next generation double-beta decay experiment called 
CUORICINO. 

Thus far the stability, the energy resolution, as well as the background counting 
rates, imply that the expanded array, CUORE, will be a very competitive next generation 
experiment. In fact, operating CUORICINO for 3 years with the present background 
level will yield a half -life sensitivity for neutrinoless double-beta decay of 6.1x10^^ 
years, corresponding to an effective mass of the electron neutrino of the order of 0.3 
eV. This will be superior to the present upper bound on the effective electron-neutrino 
mass set by the ''^Ge experiments. The early CUORICINO data have already shown us 
what and where the sources of background are and how to reduce them. We conclude 
that the required reduction in the background is possible. This would allow the full 
CUORE array to achieve a 5-year sensitivity of the order of 30 meV to the effective 
neutrino mass. Another order of magnitude reduction will present a real challenge, but is 
possible. In that case, a 5-year sensitivity in the effective mass would be around 15 meV. 
In the least optimistic case above, and in the case of inverted hierarchy of the neutrino 
mass spectrum, a null experiment would strongly disfavour the Majorana character of 
neutrinos, whereas a definite observation would confirm that neutrinos are Majorana 
particles. Only double-beta experiments can achieve this goal. The technology is now 
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proven, and with adequate funding, there are no show-stoppers. It is proposed to build 
the CUORE experiment in the Laboratori Nazionali del Gran Sasso (LNGS) in Assergi, 
Italy, the present site of CUORICINO. CUORE is by far the most cost effective next 
generation double-beta decay experiment proposed. 
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2 CUORE Science Motivation 

The neutrino mass range of interest favoured by the results of neutrino oscillation ex- 
periments is now within the grasp of well-designed PP{Ou) decay experiments. A well 
established technique has been developed by the MIBETA and CUORICINO exper- 
iments that made dramatic improvements in cryogenic spectroscopy. The realization 
that the technology is available to achieve such a fundamental physics goal provides the 
basic motivation for the CUORE experiment. 

To convey the importance of neutrino mass and of the CUORE experiment, we 
present a theoretical motivation and a brief recapitulation of past double-beta decay 
experiments, drawing heavily on our own completed work. We present our new techno- 
logical capabilities to show how the CUORE Collaboration is positioned to make rapid 
strides toward neutrino mass discovery, as well as measuring the Majorana mass scale. 
The CUORE apparatus also has application in searches for Cold Dark Matter and Solar 
Axions and more in general for rare nuclear events. 



Parameter Best fit Scr 

5.4-9.5 
1.4-3.7 

< 0.054 



Table 1: Best-fit values and 3(t intervals for the three-flavour neutrino oscillation 
parameters from global data including solar (SK and SNO), atmospheric (SK and 
MACRO), reactor (KamLAND and CHOOZ) and accelerator (K2K) experiments lE^ 

(Smjj = mf — m'j). 



5ml-^ [10-5e1/2j 6 9 

5m|2 [IQ-^el/^j 2.6 

sin^^i 0.52 0.31-0.72 

sin^ $2 0.006 

sin^ ^3 0.30 0.23-0.39 



2.1 Theoretical Motivation of Neutrinoless Double-Beta Decay Ex- 
periments 

Neutrinoless double-beta decay is by now an old subject |[D|2l|3lSl|5l|6l|7|. What is new 
is the fact that positive observation of neutrino oscillations in atmospheric |l9l [TOl HU 
fT2t . solar jnHnHISlIISlini, reactor CHI US) and accelerator [W] neutrinos gives new 
motivation for more sensitive searches. In fact, recently published constraints on the 
mixing angles of the neutrino-mixing matrix [,2L ,22 , 23 , 24 . 25 1 make a strong case that 
if neutrinos are Majorana particles, there are many scenarios in which next generation 
double-beta decay experiments should be able to observe the phenomenon and measure 
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the effective Majorana mass of the electron neutrino, | (mj^) | , which would provide a 
measure of the neutrino mass scale m. One fact is clear: neutrino oscillation experi- 
ments can only provide data on the mass differences of the neutrino mass-eigenstates. 
The absolute scale can only be obtained from direct mass measurements, end point 
measurements for example or in the case of Majorana neutrinos, more sensitively 
by neutrinoless double-beta decay. The time for large, next generation PP -decay ex- 
periments has arrived, for if the mass scale is below ~ 0.35 eV, PP -decay may be the 
only hope for measuring it. 



2.2 Neutrino Mixing Matrix 

Using the Chau and Keung parametrization of the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) neutrino mixing matrix 1.27 J : 

Z/^ = f/y Z/2 = -S3C1 - 0381826''^ C3C1 - S3SiS2e'^ S1C2 X 

Ur J \ ^3 / V S3S1 - 0301526*"^ -C3S1 - S3CiS2e*'' C1C2 / 



1 














gj02/2 





II 

















(1) 



where q = cos 6i, Si = sin 6i, V is a diagonal matrix containing Majorana CP phases 
that do not appear in neutrino oscillations. While this looks very complicated and pop- 
ulated with many unknowns, neutrino oscillation data ifTTl (151 fTTl fT^ fTHl [T^ EOl have 
constrained all three of the angles and squared neutrino mass differences as shown in 
Tabled] 

Considering the values found in TablelU we find it is possible to make the approxi- 
mation that 02 = and assume maximal mixing for the atmospheric neutrino oscillation 
data, or 6i = 45°. Accordingly, 

/ C3 S3 

U = -S3C1 C3C1 si 

\ ■S3S1 — C3S1 Ci 

where C3 = S3 = | and Ci = Si = 4j were used in the second matrix. 



2 

]_ 

\ 2V2 



\ 



1 

2 

VI 

V3 J_ I 

2V2 V2 / 



1 

V2 



(2) 



2.3 Neutrinoless Double-Beta Decay 

The decay rate for the process involving the exchange of a Majorana neutrino can be 
expressed as follows: 

l^l/2j - ^2 
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where 

Fjv = G°lM°n^ (4) 

with G^'^ being the two-body phase- space factor including coupling constants and Af "'^ 
the /3/3(Oz/) nuclear matrix element. The quantity | (m,^) | is the effective Majorana elec- 
tron neutrino mass given by: 

|(m,)| = llf/fipmi + \U^^\''m2e"t'' + If/fgl W^'l, (5) 

where e"^^ and e**^^ are the Majorana CP phases (±1 for CP conservation) and mi 2,3 
are the mass eigenvalues. The measured values of 5m|^ (Sntg solar) and ^mgg {5m\j, 
atmospheric) given in tabled motivate the pattern of masses in two possible hierarchy 
schemes shown in Figure [B 




Figure 1: Normal and inverted mass hierarchies schemes. Quasi-degenerate hierarchy 
corresponds to the case mi >> 5m aT » 6ms. 

In the case of inverted hierarchy the correct expression for \ {m^) \ can be obtained 
by interchanging the first and third columns of UV, i.e. by interchanging the roles of mi 
and m3. 

In general, prior to the approximation 62 = = S2, in the normal hierarchy case we 
have: 

\{m,)\NH = Icjclmi + slcle'^^m2 + s^e^'^^ms]. (6) 
while in the inverted hierarchy we have: 

\{m,)\iH = Is^e'-^^mi + s^c^e'^^ma + c^c^msj. (7) 
With the values and errors {3a) from Table[Tl these becomes 

\{m,)\NH = |(0.69lHt)^i + (0.30^°:°8^)e*^^m2 + (< 0.054)e*<^^m3|. (8) 

\{m,)\iH = |(< 0.054)e^^^mi + (0.30l°:°De'^^m2 + (0.69lHt)m3|. (9) 
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From the relations in Figured we can write m2 
in the case of normal hierarchy and m2 



2 

AT 



a/ 5m| + mf and = 6m\rp + 6m 
m\ and = 5m| + 5m^y 



in the case of inverted hierarchy. Equations Q and Q can be therefore expressed in 
terms of mixing angles, (5m|, 6m\j; and CP phases as follows l.28>.2LI : 



IK) I 



NH 



2 2 I 2 2 itAo 



(5m| + ml + sle'^'' 



5m\rp + (5m| + mf 



(10) 



mi.)|/H= |s2e*'^='mi + S3C2e''^S/5m^7, + mf + C3C2W5m| + 5m^j, + mf^ (11) 



With the approximation 6*2 = and the further approximation of 5m| << 5m\rp, 
equations (flOb and (fTTTl are rewritten as follows: 



|(m^)| = mi 



(5m? 



mt 



(12) 



= \ SttiIt + ml\sle"^^ + cl\. 



(13) 



Actual values of | (m^) \ for different hierarchies have been recently evaluated by 
many authors on the base of the latest results of neutrino oscillation experiments. The 
results of the analysis of Pascoli et al. are summarized for reference in tableElL2_lJ. 



sin^^ 






|(m,)|Sf^ 




0.0 


3.7 


8.7 


50.6 


47.9 


0.02 


4.6 


8.6 


49.6 


42.8 


0.04 


5.3 


9.9 


48.6 


45.4 



Table 2: Present constraints on | (m^) \ . The maximal values of | {m^) \ (in units of 10^^ 
eV) for the normal (NH) and inverted (IH) hierarchy cases, and the minimal values of 
I {m^) I (in units of 10^^ eV) for the IH and QD (quasi-degenerate) cases, for the 90% 
C.L. allowed values of 6ms, (^s and 6mAT- 



2.4 Experimental Prospects 

It is interesting to compare the values given in Table|2lwith the projected sensitivity of 
the CUORE experiment for ^^°Te /5/3(Oz/) . The projected 5 years sensitivity, discussed 
in detail later in this proposal, is 2.1x10^^ y. Using eq. |3] and El we can obtain the 
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corresponding sensitivity on | {rrii,) \ . A list of values of Fat calculated in the framework 
of different nuclear models are summarized in table |3l 

Using the average value of these values, we find that the CUORE experiment pre- 
dicts a sensitivity of ~30 meV. This value probes a significant portion of the range of 
values of mi given in the table El We note that some of the proposed experiments al- 
ready consider a possible extension to even larger masses. This is for instance the case 
for the full ten tons proposal of GENIUS. The cost of this experiment would amount 
to two billion dollars and would require a very long delay for the enrichment of such a 
large mass of ^^Ge. A similar extension to ten tons of CUORE would cost at least an 
order of magnitude less, and would avoid the long delay required by enrichment. The 
option for an "enriched CUORE" will be considered later. 

Since Srrig « Sm\j. and 62 is small (Tab.[T]and Fig.[T]), equations (fT2ll and (fT3t can 
be used to constrain the sum of the masses of all the neutrino eigenstates (S = S^mj) 
||281 

2|K)| + ± S<r < ^ < ^ ^' '''' 

where the plus is for normal hierarchy and the minus for inverted hierarchy. Equation 
(fT4ll can be simplified significantly for values of | {m^) \ achievable in next generation 
experiments. 

When Sm\rp << eq. (fT4l) transforms into 

IK)|<|<2|K)| (15) 

which, using the sensitivity projection given above, gives 90 meV < E < 180 meV. 
CUORE will therefore improve the WMAP result (S < 700 meV pT^) further con- 
straining the Majorana neutrinos contribution to the Cosmic Microwave Background 
fluctuations. 

In a similar manner, CUORE would be very complementary to the proposed KA- 
TRIN '^H neutrino mass experiment. KATRIN has a projected sensitivity of 200 meV 
ll26l . therefore should it observe a positive effect, and CUORE place a bound | (mj,) | « 
'rriKATRiN, it would be clear evidence that neutrinos are Dirac particles. 

In conclusion, according to most theoretical analyses of of present neutrino experi- 
ment results, next-generation DBD experiments with mass sensitivities of the order of 
10 meV may find the Majorana neutrino with a non-zero effective electron neutrino 
mass, if the neutrino is self-conjugate and the neutrino mass spectrum is of the quasi- 
degenerate type or it has inverted hierarchy. Majorana massive neutrinos are common 
predictions in most theoretical models, and the value of a few 10~^ eV predicted for 
its effective mass, if reached experimentally will test its Majorana nature. DBD ex- 
periments with even better sensitivities (of the order of meV) will be essential to fix 
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the absolute neutrino mass scale and possibly to provide information on CP violation. 
It is therefore evident that next generation neutrinoless double-beta decay experiments 
are the next important step necessary for a more complete understanding of the physics 
of neutrinos. In this proposal we describe the CUORE experiment, show how it could 
reach the required sensitivity, and give details of the work breakdown, cost estimate, 
and main project milestones. 
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Table 3: /9/?(0z/) nuclear factors of merit Fjv for ^^oxe and ''^Ge according to different 
evaluation methods (QRPA: Quasi Random Phase Approximation, SM: Shell Model and 
OEM: Operator Expansion Method) and authors. The foreseen P(3{0i') half-lifetime for 
i3°Te (I {rriu) \ =1 eV) is also reported. 
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3 Status of Neutrinoless Double-Beta Decay Experiments 

In the Standard Model, neutrinos and their antiparticles, antineutrinos, are different par- 
ticles, but no experimental proof has been provided so far. Nuclear double beta decay 
addresses both the question of the neutrino nature (Dirac or Majorana) and of its pos- 
sible Majorana mass. The standard nuclear Double Beta Decay is a process where a 
nucleus (A,Z) transforms into its isobar (A,Z± 2) with the emission of two electrons (or 
positrons) and their corresponding neutrinos (z/g or TQ, say: (A,Z)— >(A,Z-i-2)-i-2e^-i-2l^. 
The most direct way to determine if neutrinos are Majorana particles is to explore, 
in potential nuclear double beta emitters (A,Z), if they decay without emitting neutri- 
nos (A,Z) {A,Z+2) + 2e~, violating the lepton number conservation. For this non- 
standard /3/3(Oz/) process to happen, the emitted neutrino in the first neutron decay must 
be equal to its antineutrino and match the helicity of the neutrino absorbed by the second 
neutron. 

Phenomenologically that implies the presence of a mass term or a right-handed cou- 
pling. Well-known arguments of Kayser, Petcov, Rosen and Schechter and Valle ["481 
show that in the context of any gauge theory, whatever mechanism be responsible for the 
neutrinoless decay, a Majorana neutrino mass is required. In fact [49|, the observation 
of a /5/?(0z/) decay implies a lower bound for the neutrino mass, i.e. at least one neutrino 
eigenstate has a non-zero mass. 

Another form of neutrinoless decay, (A, Z) — > (A, Z-i-2) + 2e~ + x may reveal 
also the existence of the Majoron (x), the Goldstone boson emerging from the sponta- 
neous symmetry breaking of B-L (Baryon minus Lepton number), of most relevance in 
the generation of Majorana neutrino masses and of far-reaching implications in Astro- 
physics and Cosmology. 

Double Beta Decay (DBD) is a classical research topic, which current neutrino 
physics results have put on the front-line as an unique probe to explore and elucidate, 
in the foreseeable future, important open questions left unanswered by the oscillation 
experiments, namely the determination of the absolute neutrino mass scale. 

The DBD are extremely rare processes. In the two neutrino decay mode their half- 
lives range from T1/2 ~ 10^^ y to 10^^ y. On the other hand, the non-standard neutrino- 
less decay has not yet been found. The experimental investigation of these phenomena 
requires a large amount of DBD emitter, in low-background detectors with the capabil- 
ity for selecting reliably the signal from the background. In the case of the neutrinoless 
decay searches, the detectors should have a sharp energy resolution, or good tracking of 
particles, or other discriminating mechanisms. There are several natural and enriched 
isotopes that have been used in experiments with tens of kilograms. Some of them could 
be produced in amounts large enough to be good candidates for next generation exper- 
iments. The choice of the emitters should be made also according to its two-neutrino 
half-life (which could limit the ultimate sensitivity of the neutrinoless decay), according 
also to its nuclear factor-of-merit and according to the experimental sensitivity that the 
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Table 4: Predicted theoretical half-lives in some representative nuclear models. 



detector can achieve in the case that the emitter is, at the same time, the detector, as is 
the case of CUORE. As explained at length later on, the CUORE experiment has cho- 
sen as its (313 emitter natural tellurium oxide, as the natural isotopic abundance of^'^'^Te is 
about 34%. Finally the extremely low background required for high sensitivity double 
beta decay searches requires the development of techniques for identifying, reducing 
and suppressing the background of all types and origins in the detectors and their envi- 
ronments . All these conditions set the strategies to search for the neutrinoless double 
beta decay. 

The expected signal rate depends on the nuclear matrix element, but the dispersion 
of results of the current calculations (tab. El) makes uncertain the interpretation of the 
experimental output. 

Improvements in the precision of the theoretical evaluations of the matrix elements 
are essential. Experimental studies of nuclear structures relevant to DBD will help cal- 
culations of the matrix elements. The exploration of the conventional two-neutrino 
double beta decay of several potential double beta emitters and its comparison with 
theory must serve to help in determining the most suitable nuclear model. These im- 
provements will hopefully provide accurate nuclear matrix elements for the neutrinoless 
DBD which are crucial for extracting the effective Majorana mass parameter. Table El 
gives a summary of theoretical predictions for (3l3{2u) half-lives of several emitters in 
some representative nuclear models. 

Due to the extraordinary importance of the neutrinoless double beta decay in ex- 
tracting information on the neutrino mass scale, the nuclear matrix element issue has 
become of paramount importance both through direct calculation and through informa- 
tion obtained from nuclear reactions (of type (p,n), for instance), or through "virtual" 
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(3~ and processes which may give informations on the two branches connecting the 
relevant nuclear states in the double beta transition through the intermediate states. 

In the case of the two-neutrino double beta decay, the half-lives are customarily 
expressed as [Ty^ (0+ 0+)"^ = G2,y \ where is an integrated kinemat- 

ical factor ll56l and M^,^ the nuclear Gamow Teller matrix element. Analogously, 
the neutrinoless decay half-life (as far as the mass term contribution is concerned) is 
expressed by means of equations © and Q through the nuclear factor-of-merit F^v, 
the nuclear matrix-element M^'^, the effective neutrino mass parameter | {m,^) \ and the 
integrated kinematic factor Gq,^ f5^. 

Concerning the neutrino mass question, the discovery of a /?/5(0z/) decay will show 
that at least one of the Majorana neutrino eigenstates has a mass equal or larger than 
I {rriiy) I = rrie/ {F^Ty^Y^"^ eV, where Ty^ is the neutrinoless half life. On the contrary, 
when only a lower limit of the half-life is obtained, one gets only an upper bound on 
I (m^) I , but not an upper bound on the masses of any neutrino. In fact, | {m^) \ can 
be much smaller than the actual neutrino masses. The | {m^) \ bounds depend on the 
nuclear model used to compute the PP{Ou) matrix element. The PP{2u) decay half- 
lives measured till now constitute bench-tests to verify the reliability of the nuclear 
matrix element calculations which, obviously, are of paramount importance to derive 
the Majorana neutrino mass upper limit. As stated above the wide range of the nuclear 
matrix elements calculations for a given emitter is the main source of uncertainty in the 
derivation of the neutrino mass bound. 

3.1 Main strategies for Double Beta Decay searches 

The experimental signatures of the nuclear double beta decays are in principle very 
clear: in the case of the neutrinoless decay, one should expect a peak (at the Qf^p value) 
in the two-electron summed energy spectrum, whereas two continuous spectra (each 
one of well-defined shape) will feature the two-neutrino and the Majoron-neutrinoless 
decay modes (the first having a maximum at about one third of the Q value, and the latter 
shifted towards higher energies). In spite of such characteristic imprints, the rarity of 
the processes under consideration make their identification very difficult. Such remotely 
probable signals have to be disentangled from a background due to natural radioactive 
decay chains, cosmogenic-induced activity, and man-made radioactivity, which deposit 
energy on the same region where the PP decays do it but at a faster rate. Consequently, 
the main task in -decay searches is to diminish the background by using the state- 
of-the-art ultra-low background techniques and, hopefully, identifying the signal. 

To measure double beta decays three general approaches have been followed: geo- 
chemical, radiochemical, and direct counting measurements. In the geochemical ex- 
periments, isotopic anomalies in daughters of decaying nuclei over geological time 
scales are investigated. Determination of the gas retention age of the ore is important. 
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They are inclusive 2u + Ou measurements, not distinguishing 2u from Ou modes. How- 
ever, when T^j2,exp'~^ ^ T^J^Zpfound' ^ost of the decay is through 2u mode. The 
finite half-lives measured geochemically in the cases of ^^Se, ^^Zr, i28,i30-pg ^j^^^ radio- 
chemically in the case of ^'^^U can be regarded as I3[3{2u) half-life values. Also the 
Ti/2 values measured geochemically can be taken as a bound for T^j2^ because T^J^ (or 
the half-life of whatever decay mode) cannot be shorter than Ti/2,exp.- Consequently, 
bounds on | {m„) \ can be derived from geochemical half-life measurements. 

Another way to look for double beta decay are the radiochemical experiments, by 
noticing that when the daughter nuclei of a double beta emitter are themselves radioac- 
tive, they can be accumulated, extracted and counted. If the daughter has a half-life 
much smaller than lO^y and has no other long-lived parents, its presence can be only 
due to PP . Noticeable examples are that of '^^^U "^^^Pu (88 y, a decay) and 
244 ^ 244(-.^ ^ decay). 

Most of the recent activity, however, refers to direct counting experiments, which 
measure the energy of the PP emitted electrons and so the spectral shapes of the 2z/, 
Ou, and Oux modes of double beta decay. Some experimental devices track also the 
electrons (and other charged particles), measuring the energy, angular distribution, and 
topology of events. The tracking capabilities are useful to discriminate the PP signal 
from the background. 

The types of detectors currently used are lISTI : 

• Calorimeters where the detector is also the PP source ||5^ (Ge diodes, scintillators 
- CaF2, CdW04 -, thermal detectors -Te02 -, Xe ionization chambers, . . . ). They 
are calorimeters which measure the two-electron sum energy. Notable examples 
are IGEX, Heidelberg-Moscow, Milano-Gran Sasso (MiDBD), CUORICINO and 
the proposed CUORE and Majorana. 

• Tracking detectors of the inhomogeneous (sourceT^detector) type (Time Projec- 
tion Chambers TPC, drift chambers, electronic detectors). In this case, the PP source 
plane(s) is placed within the detector tracking volume, defining two -or more- de- 
tector sectors. Leading examples of tracking devices are the Irvine TPC's and the 
NEMO and ELEGANTS series. 

• Tracking calorimeters They are tracking devices where the tracking volume is 
also the pp source, for example a Xenon TPC (CALTECH/PSI/Neuchatel) and 
the future EXO. 

Well-known examples of PP emitters measured in direct counting experiments are 
^^Ca, ^^Ge, 96zr, ^^Se, lo^Mo, i^^cd, i^oxe, '^^Xe, i^ONd. 

The strategies followed in the PP searches are different. Calorimeters of good en- 
ergy resolution and almost 100% efficiency (e.g. Ge-detectors and tellurite bolometers) 
are well suited for PP{Ou) searches. They lack, obviously, the tracking capabilities to 
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identify the background on an event-by-event basis but they have, in favour, that their 
sharp energy resolutions do not allow the leakage of too many counts from ordinary 
double beta decay into the neutrinoless region, and so the sensitivity intrinsic limitation 
is not too severe. On the contrary, the identification capabilities of the various types of 
chambers make them very well suited for PP{2u) and /3/3(0t/, x) searches. However, 
their energy resolution is, at least presently, rather modest and the efficiency is only of a 
few percent. Furthermore, the ultimate irreducible background source in these devices 
when looking for /3/3(Oz/) decay will be that due to the standard /3/3(2z/) decay. The 
rejection of background provided by the tracking compensates, nevertheless, the figure 
of merit in PP{Ou) searches. 

Modular calorimeters can have reasonable amounts of PP emitters (Heidelberg/Moscow, 
IGEX, MiDBD and CUORICINO experiments) or large quantities (like CUORE, Ma- 
jorana and GENIUS). Tracking detectors, instead, cannot accommodate large amounts 
of PP emitters in the source plate. Recent versions of tracking devices have 10 kg and 
more (NEM03). On the other hand, TPC devices are planning to reach one ton and 
more of Xenon (EXO). 

The general strategy followed to perform a neutrinoless double beta decay experi- 
ment is simply dictated by the expression of the half-life 

^ /n2 X ^ (16) 

where N is the number of PP emitter nuclei and S the number of recorded counts 
during time t (or the upper limit of double beta counts consistent with the observed 
background). In the case of taking for S the background fluctuation (when that might 
be possible) one has the so-called detector factor-of-merit or neutrinoless sensitivity 
(Icr) which for homogeneous (source=detector) devices - for which the background 
rate scales with the detector mass - reads 

Fd = 4.17 X 10'^%f/A){Mt/BTy^^er years (17) 

where B is the background rate (c/(keV kg y)), M the mass of PP emitter (kg), €r the 
detector efficiency in the energy interval T around Q^j^ (F = FWHM) and t the running 
time measurement in years (f is the isotopic abundance and A the mass number). A 
slightly different behaviour is expected when the background rate is very low ("zero 
background" or BTMt « I hypothesis) 

Fd ~ fMter/A years (18) 

The other guideline of the experimental strategy is to choose, according to equation 
dSj), a PP emitter of large nuclear factor of merit F^ = Gqu \ M^" where the 
kinematical factor qualifies the goodness of the Q^^ value and M'^^ the likeliness of the 
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transition. Notice that the upper limit on | {m^) \ is given by | (m^) \ < nie/ {FdFnY/"^ , 
or in terms of its experimental and theoretical components 

|(m,)| < 2.5 X 10-^ X {A/ff/^ x {BT/Mtf^ x e,^'^ x G^l'^x \ M""" ["^ eV (19) 

3.2 Overview of main experiments and synopsis of results 

It is out of the scope of this Proposal to fully review the experimental status of the field. 
We refer the reader to various recent or completed reviews, also including comparison 
with theory IDEIEISIIIIEIIIIHIEIIEI. We will sketch in the following a briefing 
of the main techniques and results and summarize the current experimental situation in 
Tables |5] and |6l We will obviously devote more attention to the bolometric experiments 
with Te02, MiDBD and CUORICINO as the antecedents of CUORE. 

The Germanium case, as prototype of the calorimetric approach where the PP source 
is at the same time detector medium, has a special meaning for this Proposal. The search 
for the nuclear double beta decay of the isotope ^^Ge started as earlier as in 1967 by the 
Milano group f60', "61 1 which is now leading the CUORE Collaboration. Other groups 
both from USA (USSB-LBL-UCB, Caltech, PNL-USC) and Europe (Zaragoza, PSI- 
Neuchatel, H/M) continued this search with germanium detectors. The mastering of the 
techniques in the fabrication of Ge detectors and the continuous reduction of radioactive 
background achieved steadily since the first prototypes (PNL-USC and USSB-LBL) 
made the Ge detectors a reference in the field of double beta decay and other rare event 
physics phenomena. The effective neutrino mass bound obtained from the neutrinoless 
half life limits achieved with Ge detectors has been steadily improving along the last 
decades and now, the current most stringent bounds are led by the results obtained with 
Germanium experiments (IGEX, H/M). As they are now classical reference results for 
neutrinos mass bounds, we will describe them even briefly. 

Two recent experiments have looked for the double beta decay of ^^Ge. They both 
employ several kilograms of enriched ''^Ge (86%) in sets of detectors: the IGEX Col- 
laboration experiment (a set of three detectors of total mass 6.0 kg) in the Can- 
franc Underground Laboratory (Spain) and the Heidelberg/Moscow Collaboration ex- 
periment (63] (a set of five detectors amounting to 11 kg) in Gran Sasso and 2.1 kg in 
the Baksan Neutrino Observatory (Russia). Both experiments were designed to get the 
highest possible sensitivity for neutrinoless double beta decay of ^^Ge: a large amount 
of the ^^Ge isotope, in detectors of good energy resolution and very low radioactive 
background. Very efficient shieldings, both active and passive, were used in the two 
experiments. As a result, in the case of IGEX the background recorded in the energy 
region between 2.0 and 2.5 MeV is about 0.2 c/keV kg y prior to Pulse Shape Dis- 
crimination and most of the background in the relevant Qpp region of 2039 keV is 
accounted for by cosmogenic activated nuclei (^^Ge and ^°Co). Background reduction 
through Pulse Shape Discrimination successfully eliminates multisite events, charac- 
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teristic of non-f3f3 events, leading to less than ~ 0.07c/keV.kg.y. IGEX data corre- 
sponding to 8.87 kg y in ^^Ge provided (using the statistical estimator recommended 
by the Particle Data Group) a half-life lower bound of T^J^ > 1.57 x lO'^^y at 90% 
C.L. for the complete data set with PSD ll64l l65l l66l . Accordingly, the limit on the 
neutrino mass parameter is 0.33-1.31 eV. The uncertainties originate from the spread 
in the calculated nuclear structure parameters. Data from one of the IGEX detectors, 
RG-3 -which went underground in Canfranc several years ago- corresponding to 291 
days, were used to set a value for the 2z/-decay mode half-life by simply subtracting 
the MC-simulated background. The best fit to the stripped data corresponds to a half- 
life T^^2 = (1-45 ± 0.20) X 10^^ y. The Heidelberg/Moscow experiment operates five 
p-type HPGe detectors of enriched ^^Ge (86%) with a total active mass of 10.96 kg, cor- 
responding to 125.5 mol of ^^Ge. The recorded background spectrum of the detectors 
of the H/M experiment corresponding to an exposure of 47.4 kg.y substracted by the 
MonteCarlo simulation of the background, provides a half-life for the /3/3(2z/) -decay at 
68% C.L. T^^^ = {1.55 ± OM{stat)toilisyst)} x lO^^y in agreement with the IGEX 
result. To derive the neutrinoless decay half-life limit from the H/M experiment, the raw 
data of all five detectors as well as data with pulse shape analysis are considered. No 
indication for a peak at the Q-value/3/3(Oz/) -decay is seen in none of the two data sets 
(the first 200 d of measurement of each detector were suppressed, because of possible 
interference with the cosmogenic ^^Co). Data prior to PSD show a background around 
the Q-value of (0.19 ±0.1) c/keV.kg.y similar to that of IGEX. Both background figures, 
in IGEX and H/M, together with nthe recent results of CUORICINO to be reported later 
are the lowest ever obtained without using discrimination mechanisms. The set of H/M 
data analyzed with PSD corresponds to 35.5 kg.y of exposure and its background in the 
energy region between 2000 and 2080 keV is (0.06 ± 0.01)counts/(keV.kg.y). Follow- 
ing the method proposed by PDG, the limit on the half-life is |l67| T^J^ > 1.9 x 10^^ 
y at 90% C.L. No evidence for neutrinoless DBD has been reported so far ll2l [8l fTl l58l 
with the exception of the claimed discovery of the decay of ''^Ge reported by a sub- 
set of the Heidelberg-Moscow collaboration [68 J. This claim has been contested by 
various authors ll72l l73l and also by other members of the same Heidelberg-Moscow 
Collaboration ll74l . A new analysis in favour of the previous claim has however been 
published recently 111591 l75l . After a new statistical analysis of the same set of data, 
some of the authors of the H/M collaboration conclude [68l|69l|70l|71 1 that there exists 
evidence of a neutrinoless double beta peak. This result has been very widely contested 
ll23lEai73l|74|. 

Another technique, this time of the detector-different-from- source type, employs 
Time Projection Chambers. An example of a TPC, successfully operated along many 
years, was the Time Projection Chamber of the UC Irvine group which played an his- 
torical role in the discovery of the conventional double beta decay. It was a rectangular 
box filled with helium and located underground at 290 m.w.e. (Hoover Dam). A cen- 
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tral /3(3 source plane divides the volume into two halves, whereas a magnetic field is 
placed perpendicular to the source plane. Electrons emitted from the source follow 
helical trajectories from where the momentum and the angles of the /5-particles are de- 
termined. The PP signal is recognized as two electrons emitted from a common point 
in the source, with no other associated activity, during some time before and after the 
event. The PP source is thin enough (few mg/cm^) to allow a-particles to escape and 
be detected for tagging the background. The UCI TPC series of experiments, finished 
more than five years ago, had the merit of having obtained the first direct observation of 
the double beta decay (^^Se). They also successfully measured the two-neutrino double 
beta decay of ^°°Mo, ^^^Nd and ^^Ca. 

The NEMO series of detectors are iTT^rTTHTHI electron tracking devices (with open 
Geiger cells) filled with helium gas. An external calorimeter (plastic scintillator) covers 
the tracking volume and measures the P energies and time of flight. The PP source 
is placed in a central vertical plane and is divided in two halves, one enriched and 
another of natural abundance (of about 150 grams each), to monitor and subtract the 
background. To identify a PP signal, one should have a 2e-track with a common ver- 
tex {cosa < 0.6) in the source plus two fired plastic scintillators (E deposition>200 
keV each). The two-electron events are selected by time of flight analysis (in the en- 
ergy range of PP ). For instance, NEMO 2 has been operating for several years at the 
Modane Underground Laboratory (Frejus Tunnel) at 4800 m.w.e and has measured the 
PP{2iy) decays of i°°Mo, "*^Cd, ^^Se and ^^Zr. A new, bigger detector of the NEMO 
series, NEMO 3, has recently started its data taking with sources of up to 10 kg of ^°°Mo 
and other emitters [79|. 

The ELEGANTS V detector of the University of Osaka (placed successively in 
Kamioka and Oto) is an electron tracking detector which consists of two drift cham- 
bers for /3-trajectories, sixteen modules of plastic scintillators for P energies and timing 
measurement, and twenty modules of Nal for X- and 7-rays identification. The PP sig- 
nals should appear as two tracks in the drift chamber with the vertex in the source plus 
two signals from two plastic scintillator segments. Both enriched and natural sources (of 
about 100 grams) are employed in the detector for background monitoring and subtrac- 
tion. This detector has measured the PP{2u) decay of ^^^Cd, ^°°Mo, and several other 
emitters [fSDI . A new variant of ELEGANTS is searching for the double beta decay of 
^^Ca. 

The Caltech/PSI/Neuchatel Collaboration lEH [82| has investigated the double beta 
decay of ^^^Xe in the Gotthard Tunnel (3000 m.w.e.) by using a time projection chamber 
where the Xenon, enriched up to 62.5% in ^^^Xe, with a total mass of m=3.3 kg, was 
both the source and the detector medium, i.e. a calorimeter plus a tracking device. It 
had a cylindrical drift volume of 180 fiducial litres at a pressure of 5 atm. The PP signal 
appears as a continuous trajectory with distinctive end features: a large angle multiple 
scattering and increase charge deposition (charge "blobs") at both ends. As usual, the 
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(3(3 topology gives powerful background rejection, leading to a figure of B ~ 10 — 
10-2 c/keV kg y (at 2480 keV). 

Other results worth mentioning are those of the ITEP and Florence-Kiev groups. 
The ITEP group measured JH the double beta decay of ^^°Nd (40 g) with a TPC of 
~ 300 litres filled with CH4 at atmospheric pressure, in a 700 gauss magnetic field. On 
the other hand, the group of INR at Kiev [84| is investigating the double beta decay of 
^^^Cd with cadmium tungstate (^^^CdW04) scintillator crystals. 

To conclude this synopsis. Table |5l summarizes the neutrinoless half-life limits de- 
rived in the various experiments. From such neutrinoless half-life limits, one can derive 
Majorana neutrino mass bounds according to the nuclear model chosen. Tabled shows 
the range of bounds derived from the most sensitive experiments. 
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Table 5: Limits on Neutrinoless Decay Modes 



3.3 Thermal detector approach to double beta decay 

A series of bolometer experiments has been carried out by the Milan group since 1989 
in the Gran Sasso Laboratory, searching for the double beta decay of ^^oxe (sec. |6l). The 
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Experiment \{^u)\ < (eV) 

IGEX enrich. ^^Ge (0.33 ~ 1.35) (6.0kg) 

H/M enrich. ^^Ge (0.35 ~ 1.05) (1 1kg) 

MiDBD nat. ^^^Te (0.62 ~ 3.1) (6.8kg) 

CUORICINO nat. ^^oxe (0.38 ~ 2.2) (40.7kg) 



Table 6: Current best constraints (upper limits) on | (mj,) | . 



working principle of these thermal detectors is explained in section |5] of this Proposal. 
This saga of experiments, their development and achievements, are the origin of the 
present CUORE proposal. The increase in temperature produced by the energy released 
in the crystal due to a nuclear event (i.e. ), is measured by means of a sensor in 
thermal contact with the absorber. The Milano group has been using Tellurium oxide 
crystals as absorbers, and glued NTD Ge thermistors as sensors. Natural Tellurium 
contains 34% of ^^^Te. After using successively Te02 crystals of 73 g and 334 g, as 
well as a set of four of these large crystals, a tower-like array of 20 crystals of 340 g in 
a copper frame at ~ 10 mK was installed and operated: that was the so-called MiDBD 
experiment which successfully got the second best neutrino mass bound (after the Ge 
resuks). An enlarged version of MiDBD, called CUORICINO (meaning small CUORE) 
is now taking data in Gran Sasso (see l6.2l and l6.3l) and the limits on the effective neutrino 
mass obtained in two months only are already the best except from those of the Ge 
experiments (Table |^. 

3.4 Next generation experiments 

The requirements for a next-generation experiment can easily be deduced by reference 
to equations ([T6b and (flTt . To improve the sensitivity of \{m,,)\ by a factor of 100, 
the quantity Nt/S must be increased by a factor of 10^. The quantity N can easily be 
increased by a factor of ~ 10^ over present experiments, so that t/S must also be im- 
proved by that amount. Since practical counting times can only be increased by a factor 
2-4, the background should be reduced by a factor 25-50 below the present levels. These 
are approximately the target parameters of the next-generation neutrinoless double-beta 
decay experiments. 

The CAMEO proposal involves placing isotopically enriched parent isotopes at the 
center of BOREXINO. One example given involves 65 kg of ^^^CdW04 scintillation 
crystals. The collaboration predicts a sensitivity of | (m,,) \ ~60 meV, and with 1000 kg 
the prediction is | {m^) \ ~20 meV [ 185) . 

As already stated in sec. |2l CUORE (sec. |71) is a proposed cryogenic experiment 
with 19 towers of 52 detectors, each one being a 750 g Te02 bolometer. It will be the 



3 STATUS OF NEUTRINOLESS DOUBLE-BETA DECAY EXPERIMENTS 



29 



only future detector not requiring isotopic enrichment. This means that it would utilize 
natural abundance Te, containing 33.8% ^^^Te. A detailed description of CUORE, its 
performances, sensitivity and physics prospects is given in the present Proposal. A pilot 
experiment CUORICINO comprising one CUORE tower, is now taking data at LNGS 
(section lOb . With equivalent background, CUORE would be as sensitive as 400-950 
kg of Ge enriched to 86% ^^Ge, depending on the nuclear matrix elements used to derive 
I {m^) I . It will be performed in the Laboratori Nazionali del Gran Sasso (LNGS). It has 
been already approved by the Scientific Committee of this laboratory and its location in 
Hall A of the same Laboratory has been established. 

EXO is a large proposed TPC, either high-pressure gas or liquid, of enriched ^'^^Xe. 
This novel technique involves schemes for locating, isolating, and identifying the daugh- 
ter ^^^Ba+ ion by laser resonance ionization spectroscopy. EXO is planned as a zero- 
background experiment and it is expected to be one of the relevant experiments of the 
future. Its plans include the use of 1 to 10 tons of Xenon. A program of research and 
development is underway at the Stanford Linear Accelerator (SLAC) |86|. 

GENIUS is a proposal to use between 1.0 and 10 tons of "naked" germanium detec- 
tors, isotopically enriched to 86% in ^^Ge, directly submerged in a large tank of liquid 
nitrogen functioning both as a cooling medium and a clean shield. Extensive studies 
were made based on certain assumptions, measurements, and Monte Carlo simulations. 
In Ref. Il87ll88l . the authors claim a sensitivity range of 1-10 meV for \ {m^)\ , using 
10"^-10^ kg of enriched Ge. A research and development program is underway in the 
Gran Sasso Laboratory to develop the techniques for cooling and operating "naked" Ge 
detectors in liquid nitrogen for extended periods I 89t.90il . 

The Majorana Project is a proposed significant expansion of the IGEX experiment, 
utilizing newly developed segmented detectors along with pulse-shape discrimination 
techniques. It proposes 500 fiducial kg of Ge isotopically enriched to 86% in ^^Ge 
in the form of 200-250 detectors. Each detector will be segmented into electrically 
independent volumes, each of which will be instrumented with the new PSD system. 

The MOON experiment is a proposed major extension of the ELEGANTS experi- 
ment. It will utilize between 1 and 3 tons of Mo foils, isotopically enriched to 85% in 
loojyjQ^ inserted between plastic scintillators. It will have coincidence and tracking capa- 
bilities to search for /3/3(Oz/) decay as well as solar neutrinos. This novel technique for 
detecting solar neutrinos depends on the special properties of the nuclear decay schemes 
of ^°°Mo and its daughters, allowing both event and background identification ll9Tl . 

This list of proposals should produce several experiments with sensitivities sufficient 
to actually observe /?/3(0z/) decay or obtain stringent upper bounds on | {rrijj) \ reaching 
the sensitivity range implied by recent neutrino oscillation results. The attainability of 
such a goal strongly depends on the true capability of these projects to reach the re- 
quired background levels in the /3/3(Oz/) region. An experimental confirmation of the 
(often optimistic) background predictions of the various projects (even if extrapolated 
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from the results of lower scale successful experiments) is therefore worthwhile and the 
construction of preliminary test setups is absolutely needed. 

Frequently appearing publications on the subject almost always refer to the impor- 
tance of conducting next-generation neutrinoless double-beta decay experiments. A 
complete understanding of the neutrino mass matrix depends on three types of data: 
neutrino oscillations, tritium beta-decay measurements, and neutrinoless double-beta 
decay. Each is analogous to one leg of a three-legged stool, and each is necessary for 
the complete picture. The case for a significant investment in next-generation experi- 
ments of all three types is compelling. 
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4 CUORE experimental prospects 

As discussed in the previous sections, a number of recent theoretical interpretations of 
atmospheric, solar and accelerator neutrino experiments imply that the effective Majo- 
rana mass of the electron neutrino, | {m^) \ , expressed in Eq. |5l could be in the range 
0.01 eV to the present bounds. 

Considering this range, could a next generation PP{Ou) - decay experiment detect 
it? If so, what technique would be the best for a possible discovery experiment? We will 
address these questions in an effort to demonstrate that CUORE, an array of 988, 750 g 
Te02 bolometers, is one of the best approaches presently available. It can be launched 
without isotopic enrichment nor extensive Research and Development (R&D), and it 
can achieve next generation sensitivity. 

The CUORE project originates as a natural extension of the succesfuU MiDBD and 
CUORICINO ^-^^Te experiments, described before, where for the first time large arrays 
of bolometers were used to search for /3/3{0iy) - decay. The good results obtained so 
far prove that the bolometric technique although novel is competitive and alternative to 
the traditional calorimetric Ge technique. The pertinent details of the detectors of the 
CUORE array and its operation at temperature of ~ 10 niK, as well as the background 
issues, electronics, DAQ and data-analysis are discussed in other parts of this proposal. 
In the following we will discuss the Physics potential of CUORE regarding Double Beta 
Decay, Dark Matter, Solar Axions searches, and other rare processes. 

Rare event hysics is playing a significant role in Particle Astrophysics and Cosmol- 
ogy. Examples of such rare phenomena could be the detection of non-baryonic particle 
dark matter (axions or WIMPs), supposedly filling a substantial part of the galactic ha- 
los, nucleon decay or neutrinoless double beta decay. These rare signals, if detected, 
would be important evidence of new physics beyond the Standard Model of Particle 
Physics, and would have far-reaching consequences in Cosmology. The experimental 
achievements accomplished during the last decade in the field of ultra-low background 
detectors have led to sensitivities capable of searching for such rare events. Dark matter 
detection experiments have largely benefited from the techniques developed for dou- 
ble beta decay searches. We expect a similar symbiosis also in the case of CUORE 
and CUORICINO for which Dark Matter detection will be therefore another important 
scientific objective. 

Due to the low rates of the proposed processes, the essential requirement of these 
experiments is to achieve an extremely low radioactive background and large mass. Ac- 
cordingly the use of radiopure detector components and shieldings, the instrumentation 
of mechanisms of background identification, the operation in an ultra-low background 
environment, in summary, the use of the state-of-the-art of low background techniques 
is mandatory. In some of these phenomena, as in the case of the interaction of a particle 
dark matter with ordinary matter, a very small amount of energy is deposited, and the 
sensitivity needed to detect events within such a range of energies, relies on how low 



4 CUORE EXPERIMENTAL PROSPECTS 



32 



the energy threshold of detection is. In addition, to increase the chances of observing 
such rare events, a large amount of detector mass is in general needed, and so most 
of the experiments devoted to this type of searches are planned to have large detector 
mass, while keeping the other experimental parameters (background, energy thresholds 
and resolutions) optimized. On the other hand, the recent development of cryogenic 
particle detection has led to the extended use of thermal detectors to take 

advantage of the low energy threshold and good energy resolution required for the theo- 
retically expected thermal signals. This detection technology has also the advantage of 
enlarging the choice of materials which can be used, either as DM targets or PP decay 
emitters. After a long period of R&D to master the techniques used in cryogenic particle 
detectors, low temperature devices of various types are now applied to the detection of 
double beta decay or particle dark matter ll95l . A major example of this development 
were the MiDBD (Milano Double Beta Decay) experiment (sec. 16.21) . which success- 
fully operated a 20 Te02 crystal array of thermal detectors of a total mass of 6.8 kg and 
CUORICINO discussed below. With the objective of going to larger detector masses 
and of improving the sensitivity achieved in the smaller arrays, the CUORE {Cryogenic 
Underground Observatory for Rare Events) project 1%", W] was bom some years ago 
as a substantial extension of MiDBD. The objective of CUORE is to construct an array 
of 988 bolometric detectors with cubic crystal absorbers of tellurite of 5 cm side and of 
about 750 g of mass each. The crystals will be arranged in a cubic compact structure and 
the experiment will be installed in the Gran Sasso Underground Laboratory (LNGS). 

The material to be employed first is Te02, because the main goal of CUORE is to 
investigate the double beta decay of ^'^°Te, although other absorbers could also be used 
to selectively study several types of rare event phenomenology as it will be discussed 
later. Apart from the MiDBD experiment, already completed, a smaller and intermediate 
stage of CUORE, CUORICINO is presently taking data at LNGS. It consists of 44 cubic 
crystals slighly larger with respect to those foreseen for CUORE (790 g) and 18 smaller 
crystals (330 g) for a total mass of ~ 40.7 kg. CUORICINO is, by far, the largest 
cryogenic detector on the stage. Its early results are excellent (see section l63t and 
represent the most eloquent demonstration of the CUORE feasibility. 

In the following, the Physics Potential of the CUORE experiment with respect to its 
double beta decay discovery potential and its detection capability of WIMP and axions 
is presented. For this purpose, we will rely on the CUORE background expected sensi- 
tivity described in other sections of this proposal (l6.3l andfT4ll and tested with the results 
already obtained in CUORICINO. 

In CUORICINO a background improved by a factor 1.5 with respect to the one 
measured in MiDBD-II has been observed in the /3/5(Oz/) region. The analysis of the 
low energy region, relevant for dark matter search, has not been completed but we 
expect a similar improvement there also. For what concerns CUORE on the other 
hand, MiDBD and CUORICINO have shed light on how the background contami- 
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nation dominating the two regions of interest is mostly the intrinsic background due 
to bulk and surface contaminations of the constructing materials. With the presently 
achieved quality of low contamination materials and considering the worst possible con- 
dition for bulk contaminations (i.e. all the contamination equal to the present 90%/,C.L. 
measured upper limits) we proved that the corresponding contribution to the CUORE 
background will be ~0.004 counts/(keV kg y) at the PP{Ou) transition energy and 
~0.025 counts/(keV kg d) near threshold ("section (141). Even at these conservative val- 
ues therefore, bulk contaminations will not represent a problem for CUORE sensitivity. 
A larger background contribution is however expected when extrapolating to CUORE 
the surface contributions measured in CUORICINO. A dedicated R&D aiming at re- 
ducing the radioactive surface contaminations of the detector structure (mainly copper) 
is therefore required. A reduction by an order of magnitude of the surface contributions 
seems attainable in a reasonable time (less than 2 years) by simply improving the sur- 
face cleaning procedures used for CUORICINO (e.g. by using cleaner solutions). This 
would be enough to reduce the /?/9(0z/) background contributions to the level of few 
counts/keV/ton/y. Further improvements are still possible (e.g. active surface bolome- 
ters) and are considered by the CUORE collaboration as part of a second R&D phase 
aiming at an even better sensitivity both for /3/3(Oz/) and dark matter searches (sec.fTSt. 
In particular, the goal of this CUORE technical R&D will be the reduction of back- 
ground to the level of ~0.001 counts/(keV kg y) at the (3(3{0u) transition energy and 
to ~0.01 counts/(keV kg d) at threshold. The potential of the experiment in this back- 
ground configuration will also be discussed. 

Regarding the expected threshold and resolution, energy thresholds of ~5 keV and 
energy resolutions lower than 1 keV at the 46 keV line of ^^°Pb were obtained in some 
test measurements carried out at LNGS ("Hall C" setup) for the early CUORE R&D 
program We will assume therefore conservative values of 10 keV for the energy 

threshold and of 1 keV for the energy resolutions at threshold, as observed in MiDBD 
and CUORICINO. As far as the energy resolutions obtained in the double beta decay 
region, values of 3 keV at 2615 keV were achieved in some of the above mentioned 
test measurements. On the other hand, a 9 keV average value has been measured in 
the CUORICINO calibration spectra. A possible energy resolution interval of 5-10 keV 
(taking into account possible improvements in the foreseen CUORE R&D, sec.fT5l will 
be therefore assumed for CUORE. 

Taking into account these expectations, we discuss in the following the prospects 
of CUORE for double beta decay searches, for WIMP detection, for solar axion explo- 
ration and for searches on rare nuclear events. 

4.1 Double beta decay prospects 

The main scientific objective of the CUORE detector is the search for the neutrinoless 
double-beta decay of the ^^°Te isotope contained in the (natural) Te02 crystals. 
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In fact, cryogenic thermal detectors provide several double-beta nuclei to be ex- 
plored in "active" source/detector calorimeters as it will be discussed later. Some of 
them have been tested and others are already in running detectors, like ^^Ca in CaF2, 
^^°Te in Te02, and ^^^Cd in CdW04. As far as the Tellurium Oxide is concerned, the 
^^°Te isotope is a good candidate for double beta decay searches: its isotopic content in 
natural Telluriumis 33.87%, and its PP Q-value (Q/3/3 = 2528. 8±1. 3 keV) is reasonably 
above the main radioimpurity lines when looking for a neutrinoless signal. Moreover, 
this Q-value happens to be between the peak and the Compton edge of the 2615 keV 
line of ^°^T1, which leaves a clean window to look for the signal. Finally, it has a fairly 
good neutrinoless nuclear factor-of-merit F^v (sec. 12. 31) . It is in fact apparent (tab.|3l) that, 
no matter which nuclear model is used to compute the neutrinoless double-beta decay 
matrix elements, the nuclear factors of merit of ^•^°Te are a factor 5-10 more favorable 
than those of ^^Ge (the emitter where the best neutrinoless double beta decay half-life 
limits have been achieved so far). This translates into a factor 2 to 3 better as far as the 
I {rriiy) I (Majorana neutrino mass parameter) bounds are concerned. 

The detector factor-of-merit F^^ (eq. (fTTl ). or detection sensitivity, introduced earlier 
by Fiorini, provides an approximate estimate of the neutrinoless half-life limit (la) 
achievable with a given detector. In the case of a 5x5x5 cm'^ Te02 bolometers (for 

which the /3/?(0i/) detector efficiency is 0.86) we have F^j^ ~7.59x lO^^^M^ with M 
the crystal mass in kg and B the background in counts per keV per year and per kg of 
detector mass. 

To demonstrate the full CUORE physics potential, let us evaluate F^ according to 
our conservative projection for the CUORE array (750 kg of Te02 ). As described 
above, a background of B=0.01 c/keV/kg/y would be achievable with a slight improve- 
ment of the current available material selection and cleaning techniques (sections 1 12. 11 
fT^andfT^t and exploting the possibility to build a dedicated cryostat with low activ- 
ity materials and effective shields (sec. El). Assuming finally an energy resolution r(2.5 
MeV)=5 keV, we would get ¥^j^ of 9.4x lO^^^t years (6.5 x lO'^^^/i years for r=10 keV), 
which in t years of statistics would provide | {m^) \ bounds in the range 0.036-0.2 t~^/^ 
eV (according to QRR\ models Fjv predictions, tab. 13). However, the R&D to be car- 
ried out in CUORE, if successful, would provide a value of B~0.001 c/keV/kg/y, i.e. a 
detection sensitivity of Fi:)~2.96x 10^^\/t years (2.1 x 10^^ years for r=10 keV), or 
I {m^) I bounds in the range ~0.02-0.1 1 1"^/"^ eV. Te02 crystals made with ^^°Te enriched 
material have been already operated in MiDBD and CUORICINO, making an enriched 
CUORE a feasible option (secfTTt. Assuming a 95% enrichment in ^^°Te and a back- 
ground level of b=0.001 c/keV/kg/y, the sensitivity would become F£)~8.32x lO^^v^ 
years. For an exposure of 5 years, the corresponding | {m^) \ bounds would range from 
8 meV to 45 meV depending on the nuclear matrix element calculations. 
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4.2 Other physics with CUORE 

Even if double beta decay, particularly in its neutrinoless channel, is the main scientific 
aim of CUORE, other interesting searches can be tackled with this massive detector. For 
instance, as stressed above, dark matter detection experiments have largely benefitted in 
the past from the techniques developed for double beta decay searches and we have valid 
motivations to expect such a symbiosis also in the case of CUORE and CUORICINO. 
A number of other searches for rare processes (rare nuclear a and P decays, electron 
and nucleon stability, etc.) could however take advantage of the large mass and low 
background of CUORE. We will discuss in the following, the potential of CUORE for 
such searches. 

4.2.1 WIMP detection 

Recent cosmological observations flOOl provide compelling evidence for the existence 
of an important component of non-baryonic cold dark matter in the Universe. Among 
the candidates to compose this matter. Weakly Interacting Massive Particles (WIMPs) 
and axions are the front runners. The lightest stable particles of supersymmetric theo- 
ries, like the neutralino HlOll . constitute a particular class of WIMPs. 

Under the hypothesis that WIMPS are the main component of the dark matter, these 
particles should fill the galactic halos and explain the flat rotation curves which are ob- 
served in spiral galaxies. The detection of such particles could be attempted both by 
means of direct and indirect methods. The direct detection of WIMPs relies on the mea- 
surement of their elastic scattering off the target nuclei of a suitable detector [.95J . The 
non relativistic and heavy (GeV - TeV) WIMPs could hit a detector nucleus producing 
a nuclear recoil of a few keV. Because of the small WIMP-matter interaction cross sec- 
tions the rate is extremely low. In the case of SUSY WIMPs, most of the cross section 
predictions HI 021 1 1031 1 10411 (derived using MSSM as the basic frame implemented with 
different unification hypothesis) encompass a range of values several orders of magni- 
tude (the so-called scatter plots) providing rates rangingfrom 1 c/kg/day down to 10~^ 
c/kg/day according to the particular SUSY model. 

It is well known that the predicted signal for the WIMP elastic scattering has an ex- 
ponentially decaying energy dependence, hardly distinguishable from the background 
recorded in the detector. The simple comparison of the theoretical WIMP spectrum 
with the one experimentally obtained, provides an exclusion curve (at a given confi- 
dence level), as dark matter component of the halo, of those WIMPs with masses (m) 
and cross sections on nucleons (a) which yield spectra above the measured experimen- 
tal rate. To claim a positive identification of the WIMP, however, a distinctive signature 
is needed. The only identification signals of the WIMP explored up to now are provided 
by the features of the Earth's motion with respect to the dark matter halo. In particular, 
the annual modulation 11051 is originated by the combination of the motion of the solar 
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system in the galactic rest frame and the rotation of the Earth around the Sun. Due to 
this effect, the incoming WIMP velocities in the detector rest frame change continu- 
ously during the year, having a maximum in summer and a minimum in winter in the 
northern hemisphere. Therefore the total WIMP rate changes in time with an oscillating 
frequency which corresponds to an annual period and a maximum around the beginning 
of June. 

The relative annual variation of the signal is small (a few percent) so in order to 
detect it one needs large detector masses to increase statistics and several periods of 
exposure to minimize systematics. Several experiments have already searched for this 
effect 111061 1 1071 1 10811 and since 1997 the DAMA group has reported a positive signal 
HI 091 which has been appearing throughout four yearly periods. Recently, the DAMA 
experiment has extended up to seven yearly periods the annual modulation effect, which 
stands now at 6.3 a level [ 1 10|. The present situation is no doubt exciting. On one hand 
that result has triggered an intense activity in the field; on the other, the experimental 
sensitivities of various types of underground detectors are entering the supersymmetric 
parameter space [ 103J and in particular several of them 11111 11121 11131 11141 are ex- 
cluding, to a larger or lesser degree (within standard framework models), the region of 
mass and cross-section where the reported WIMP is supposed to exist. This is partic- 
ularly true after the recent results of the CDMS experiment 1 1 15| whose disagreement 
with DAMA is statistically more significant. It is then compulsory to explore this effect 
with other nuclear targets, and that will be one of the purposes of CUORE. We will 
discuss in the following the capabilities of CUORE to exclude or detect WIMPs using 
the total time-integrated experimental rate and comparing it with the predicted nuclear 
recoil rate. To look for the annual modulation signal in CUORE, which in principle 
has enough mass to be sensitive to it, one needs to know its stability performance. The 
analysis of the CUORE potential for annual modulation searches will be performed - 
following statistical consideration- (see Ref. 111611 ). with the proviso that systematic 
uncertainties are under control. Data on the stability of CUORICINO will be crucial to 
assess such hypotheses. 

To calculate the theoretical WIMP rate, standard hypothesis and astrophysical pa- 
rameters are assumed, i.e., that the WIMPs form an isotropic, isothermal, non-rotating 
halo (the isothermal sphere model) of density p = 0.3 GeV/cm'^, which has a maxwellian 
velocity distribution with Vj-ms = 270 km/s (with an upper cut corresponding to an es- 
cape velocity of 650 km/s), and a relative Earth-halo velocity of Vr = 230 km/s). Other, 
more elaborated halo models, which have been considered recently 11171 would lead to 
different results. The same applies when other astrophysical parameters are employed 
or when uncertainties in the halo WIMPs velocity distribution are included 1111811 . The 
theoretical predicted rate is expressed in terms of the mass and cross-section of the 
WIMP-matter interaction. The cross sections are normalized per nucleon assuming a 
dominant scalar interaction, as is expected, for instance, for one of the most popular 
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Figure 2: Exclusion projected for 1 year of CUORE assuming a threshold of 10 keV, 
a low energy resolution of 1 keV, and low energy background levels of 0.05 and 0.01 
c/keV/kg/day respectively. The dashed line corresponds to the MiDBD result. 



dark matter candidates, the neutralino: 



0'Nx = O-nx^ — (20) 

where A is the target (oxygen and tellurium) mass number, /i^jv' '^he WIMP-nucleus 
reduced mass, and the WIMP-nucleon reduced mass. The Helm parameterization 
HI 191 is used for the scalar nucleon form factor. The (m, a) exclusion plot is then 
derived by requiring the theoretically predicted signal for each m and a in each energy 
bin to be less than or equal to the (90% C.L.) the Poisson upper limit of the recorded 
counts. The bin width is assumed to be equal to the detector resolution. 

In Fig. 121 the exclusion plots for coherent spin-independent WIMP-matter interac- 
tion are shown for two possible values of the background of CUORE: 0.05 and 0.01 
c/keV/kg/dayl and 0.1 which we expect to be attainable after a dedicated R&D directed 
to improve the detector performance in the low energy region. 

Notice, moreover, that values of a few 0.01 c/keV/kg/day have been obtained above 
10 keV in the raw spectra of Germanium experiments (e.g. IGEX II113I ) without using 
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mechanisms of background rejection, and so it does not seem impossible to achieve such 
equivalent small values in crystal thermal detectors of tellurium (with only phonons, 
and no discrimination mechanism). To draw the two exclusion contours of Fig. |2l a low 
energy resolution of 1 keV, an energy threshold of 10 keV, the two quoted values of the 
background of CUORE (0.05 and 0.01 c/keV/kg/day) and an exposure of 1 year (750 
kg-year) have been assumed. 

As previously noted, CUORE is characterized by a detector mass large enough to 
search for the annual modulation signal. As it is well known, an essential requirement to 
estimate the prospects of any detector to search for annual modulation is to have a superb 
control of systematic errors and to assure that the stability of the various experimental 
parameters, which might mimic periodic variations of the signals, are kept within a 
small fraction of the (already small) expected signal. The various changes of the set- 
up, crystals and shielding of the MiDBD experiment have not provided a definitive 
estimation of the long-term stability parameters of MiDBD; it will be however measured 
in CUORICINO. Possible instabilities are that of the electronic gain and the ensuing 
time fluctuation of the energy scale (both in energy thresholds and energy resolutions), 
the temperature variations, the possible fluctuation in time of the efficiency with which 
the triggered noise is rejected. 

The fact that we are dealing with a signal depending on time, which typically amounts 
to a fraction between 1% and 7% of the average count rates, reinforces the need for a 
control of the stability of the experiment well below that range over long periods of 
time. 

When assuming that all these fluctuations are controlled well below the levels needed 
(<1%), then one can proceed to analyze the sensitivity of CUORE to the annual mod- 
ulation signal on purely statistical grounds. This has been first attempted in HI 201 and 
II121I . but a more extensive and rigorous approach is followed in ref. II116II where sen- 
sitivity plots for several types of detectors (and experimental parameters) are presented. 

Following the guidelines of that reference, the sensitivity of a given experimental 
device to the annual modulation signal can be precisely quantified by means of the 5 
parameter, defined from the likelihood function or, equivalently, from the function of 
the cosine projections of the data (for further details see ref. II116II '): 



This parameter measures the statistical significance of the modulation signal de- 
tected in an experimental set of data. However, for a given (m, a) and a given exper- 
iment, the expected value (5^) can be estimated using the expression derived in ref. 



5^ = via = 0)-y., 



mm — 



(21) 




(22) 
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Figure 3: The solid lines represent the sensitivity plot in the (m, cr) plane for CUORE, 
assuming a threshold of 10 keV, two years of exposure (1500 kg year) and flat back- 
grounds of 0.05 and 0.01 c/keV/kg/day. It has been calculated for (5^) = 5.6 (see the 
text). The sensitivity curve has been also calculated for a possible threshold of 5 keV 
with a background of 0.01 c/keV/kg/day (dashed line). 

where Sm,k and So^k are the modulated and non-modulated parts of the WIMP signal in 
the kth energy bin of AEf^ width, is the background in that energy bin and MTa the 
effective exposure, being a a coefficient accounting for the temporal distribution of the 
exposure time around modulation maxima and minima (a = ^/nJ2^^i cos^ to{ti — to) 
for n temporal bins). 

Using this equation we have estimated the region that could be within reach for 
CUORE with the above mentioned assumptions on the background levels. We have 
fixed a value of 5.6 for (5^) that corresponds to 50% probability of obtaining a positive 
result at 90% C.L. In figure El curves are shown obtained for a threshold of 10 keV, two 
years of exposure and two assumed flat backgrounds of 0.05 and 0.01 c/keV/kg/day 
(solid lines). The possibility of a lower threshold of 5 keV with a background of 0.01 
c/keV/kg/day is also shown (dashed line). 

In conclusion, CUORE will be able to explore (and/or exclude) WIMPs lying in 
large regions of their parameter space. 

The capability of CUORE to investigate the DAMA region through the exclusion 
plot (time integrated method) relies on getting a background of 0.1 c/keV/kg/day from 
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10 keV onwards, independently of more elaborated time modulation methods which 
require an exhaustive control of the stability of the experiment. However, CUORE 
could also attempt to look for annual modulation of WIMP signals provided that the 
stability of the experiment is sufficient. 

We would like to note in addition that the nuclei of Tellurium are totally different 
from those of Germanium and close in atomic number to Iodine. The ground state of 
93% of the Te nuclei is 0+. Comparison among the results from these nuclei seems 
therefore very interesting. 

4.2.2 Solar axion detection 

Axions are light pseudoscalar particles which arise in theories in which the Peccei- 
Quinn U(l) symmetry has been introduced to solve the strong CP problem [jl22|. They 
could have been produced in early stages of the Universe being attractive candidates for 
the cold dark matter (and in some particular scenarios for the hot dark matter) respon- 
sible to 1/3 of the ingredients of a flat universe. Dark matter axions can exist in the 
mass window lO^^*^"^^ eV < rua < 10^^ eV, but hadronic axions could exist with masses 
around the eV. 

Axions could also be copiously produced in the core of the stars by means of the 
Primakoff conversion of the plasma photons. In particular, a nearby and powerful source 
of stellar axions would be the Sun. The solar axion flux can be easily estimated II123L 
1 12411 within the standard solar model, resulting in an axion flux of an average energy of 
about 4 keV that can produce detectable X-rays when converted in an electromagnetic 
field. We would like to stress that, although we focus on the axion because its special 
theoretical motivations, all this scenario is also valid for any generic pseudoscalar (or 
scalar) particle coupled to photons |ll23. Needless to say that the discovery of any type 
of pseudoscalar or scalar particle would be extremely interesting in Particle Physics. We 
will keep our discussion, however, restricted to the case of solar axions. 

Single crystal detectors provide a simple mechanism for solar axion detection HI 241 
I126II . Axions can pass in the proximity of the atomic nuclei of the crystal where the 
intense electric field can trigger their conversion into photons. The detection rate is en- 
hanced if axions from the Sun coherently convert into photons when their incident angle 
with a given crystalline plane fulfills the Bragg condition. This induces a correlation of 
the signal with the position of the Sun which can be searched for in the data and allows 
for background subtraction. 

The potentiality of Primakoff conversion in crystals relies in the fact that it can 
explore a range of axion masses (m^ > 0.1 keV) not accessible to other direct searches. 
Moreover it is a relatively simple technique that can be directly applied to detectors 
searching for WIMPs. 

Primakoff conversion using a crystal lattice has already been employed in two ger- 
manium experiments: SOLAX II127II and COSME-II II128II with the ensuing limits for 
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axion-photon coupling ga-y-y < 2.7 x 10~^ GeV~^ and ga^^ < 2.8 x 10~^ GeV~^ re- 
spectively. Also the DAMA collaboration has analyzed 53437 kg-day of data of their 
Nal set up ['1291, in a search for solar axions, following the techniques developed in ref. 
lll30l . where a calculation of the perspectives of various crystals detectors (including 
Nal) for solar axion searches has been made. The DAMA result ga-y-y < 1.7 x 10^^ 
GeV^^ improves slightly the limits obtained with other crystal detectors II1271I1281 and 
agrees with the result predicted in ref. 1 130|. These "crystal helioscopes" constraints are 
stronger than that of the Tokyo axion helioscope I.131J for rua > 0.26 eV and do not 
rely on astrophysical considerations (i.e. on Red Giants or HB stars dynamics [132|). 
The orientation of the crystal was not known so that the data were analyzed taking the 
angle corresponding to the most conservative limit. 

It has been noted that the model that yields the solar axion fluxes used to calculate the 
expected signals is not compatible with the constraints coming from helioseismology if 
ga^yy > 10~^ GeV^^ [ 133 1. This would imply a possible inconsistency for solar axion 
limits above that value, and sets a minimal goal for the sensitivity of future experiments. 

The use of CUORE to search for solar axions via Bragg scattering should have some 
advantages with respect to germanium detectors, because of the higher Z, larger mass 
and known orientation of the crystals. Since the cross-section for Primakoff conversion 
depends on the square of the atomic number, Te02 will be a better candidate than Ger- 
manium. Needless to say that a low energy threshold is mandatory because the expected 
signal lies in the energy region 2 keV < E < 10 keV and is peaked at ~ 4 keV. As 
for dark matter searches, this will of course require specific R&D directed to a better 
performance in the low energy region. 

A detailed analysis has been performed HI 301 for a Te02 crystal (which has a tetrag- 
onal structure II134I ') assuming different values for the experimental parameters. As it 
is shown in Ref. II130I . the bound on axion-photon coupling which a given experiment 
can achieve can be estimated through the expression: 



where k depends on the crystal structure and material, as well as on the experimental 
threshold and resolution. For the case of Te02 and a threshold of 5 keV, k has been 
calculated to be A; = 2.9 assuming an energy resolution of 1 keV. The computation of 
this expression for some assumed values of the experimental parameters is shown in 
table|7lfor CUORE. Flat backgrounds and 2 years of exposure are assumed. 

It is worth noticing the faible dependence of the ultimate achievable axion-photon 
coupling bound on the experimental parameters, background and exposure MT: the 1/8 
power dependence of gayy on such parameters softens their impact in the final result. 
The best limit shown in table [7\ is in fact only one order of magnitude better than the 
present limits of SOLAX and COSME-II. The ga-y-y bound that CUORE could provide 




(23) 
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Figure 4: Best bound attainable with CUORE (straight Une labelled "CUORE") com- 
pared with others limits. 
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Mass 


Resolution 


Threshold 


Background 


9a^^ (2 years) 


(kg) 


(keV) 


(keV) 


(c/kg/keV/day) 


(GeV-i) 


750 


1 


5 


0.01 


6.5x10-^0 


750 


1 


5 


0.05 


8.0x10^10 



Table 7: Expected limits on the photon-axion coupling for 2 years of exposure of 
CUORE assuming the quoted values for the experimental parameters 

The limit which can be expected from CUORE is better than the helioseismological 
bound mentioned before (see Figure |4| and Table |7l). Notice that in both cases an energy 
threshold of Ethr ~ 5 keV (and resolution of ~1 keV) has been assumed. As was de- 
scribed at length in Ref. II130L the crucial parameters for estimating the perspectives on 
solar axion detectors with crystals rely on the energy threshold and resolution (appear- 
ing in k), and the level of background achieved (although the influence of this parameter 
is damped by a factor 1/8). In particular, a threshold of 5-8 keV would loose most of 
the axion signal. Other crystal detectors with, say, Ge or Nal (GENIUS, MAJORANA, 
GEDEON, DAMA, LIBRA, ANAIS, . . .) could surpass CUORE as axion detectors 
because the energy thresholds of these projects are supposed to be significantly lower. 
Also the background is expected to be better. 

Nevertheless, it should be stressed that the bounds on gajj obtained with this tech- 
nique in the various proposed crystal detector arrays stagnate at a few xlO~^° GeV^^, 
not too far from the goal expected for CUORE, as has been demonstrated in [|130|. 
There are no realistic chances to challenge the limit inferred from HB stars counting 
in globular clusters II132I and a discovery of the axion by CUORE would presumably 
imply either a systematic effect in the stellar-count observations in globular clusters or 
a substantial change in the theoretical models that describe the late-stage evolution of 
low-metallicity stars. To obtain lower values of ga-y-y one should go to the magnet he- 
lioscopes like that of Tokyo II131I and that of CERN (CAST experiment II135I currently 
taking data). 

In particular, until recently the best experimental bound of ga^^ published came 
from the Tokyo helioscope: gayy < 6 x 10"^° GeV"^ for nia < 0.03 eV and ga-y-y < 
6.8 - 10.9 X 10-^° GeV"^ form^ ~ 0.05-0.27 eV. The sensitivity of CAST is supposed 
to provide a bound gayy < 5 x 10^^^ GeV^^ or even lower. A recent run of CAST has 
improved the Tokyo limit with a bound gayy < 3 x 10^^° for nia < 0.05 eV. 

4.2.3 Searches for other rare processes 

Searches for rare atomic or nuclear events have provided and are going to provide a 
very rich harvest of important results, some of them probably beyond our imagination. 
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We believe that the technique of massive thermal detectors is particularly suitable to 
search for many of them, but we present here only a few examples just to stimulate our 
theoretical friends to suggest others. 

Electron stability Searches on this decay II136II , can represent an important by-product 
for CUORE. An inclusive experiment on this process can be carried out by detecting the 
disappearance of an atomic electron followed by the consequent decay (or series of de- 
cays) of the outer-shell electrons to fill the vacancy left by this disappearance. These 
searches, which have relevance also in connection with charge conservation and Pauli 
principle validity, can be carried out by looking for sharp lines corresponding to decays 
in the atomic tellurium structure. The exclusive charge violating decay of an electron in 
a neutrino and a photon can be revealed by the presence of peak at 255 keV, spread how- 
ever by the Doppler broadening of the atomic electron motion. Even if not fully exploit- 
ing the excellent energy resolution of CUORE, this decay can be efficiently searched, 
due to the large mass and low background expected for this detector. 

Nucleon stability CUORE cannot obviously compete with the various experiments 
specifically designed to investigate specific channels of this decay II137II . It can how- 
ever provide competitive results on the inclusive decays of a proton or neutron through 
the search for delayed coincidences between the nucleon decay signature (e.g. high 
energy releases by nucleon decay products or prompt 7's from nuclear re- arrangement 
processes) and possible unstable produced nuclei. Possible nucleon decays leading to 
the production of unstable nuclei emitting gamma rays are listed in Table [8]. Delayed 
gamma rays emitted in these secondary decays could be detected in nearby crystals. If 
the half-lifetime of the unstable produced nuclei would be reasonably short an effective 
coincidence selection could be applied. 

Unusual decay modes (e.g. n — > 3z/) have been recently suggested as possible dominant 
modes II138II . The CUORE unique capabilities for the detection of prompt low energy 
7-rays following such "invisible" decay channels could provide competitive sensitiv- 
ities for the search of these decay modes which may have been missed by previous 
dedicated experiments. In particular, results from CUORE can be complementary to 
those recently obtained by the SNO collaboration. 

Rare beta processes 

^^^Cd An experiment carried out with a large CdW04 crystal HI 391 has determined the 
transition energy and lifetime of this nucleus and measured for the first time the 
shape of the beta spectrum. The study of this transition was also repeated later 
by exploiting the scintillating properties of CdW04 [ 140 1. These results, can be 
further improved if a large CdW04 would be inserted into the CUORE detector. 



4 CUORE EXPERIMENTAL PROSPECTS 45 

Table 8: Possible nucleon decays leading to unstable nuclei characterized by the emis- 
sion of gamma rays. 



Nucleus Abundance % Decay Product Lifetime Gamma Rays (keV) 





0.096 


p 


ii^Sn 


38.1 h 


24 






n 


ii9Te 


16 h 


644. 








119m'-pg 


4.69 d 


154, 1213,271 


i22Te 


2.6 


n 


i2iTe 


16.8 d 


573, 508 








121mj'g 


154 d 


1102, 37 


i23Te 


0.908 


p 


122Sb 


2.72 d 


564 








122msb 


4.191 m 


61,76 


i24Te 


4.816 


n 


123m'-pg 


119.7d 


159 


i25xe 


7.14 


p 


124Sb 


60.2 d 


603, 1691,773 








124ms|, 


20.2 m 


603, 946, 498 


126jg 


18.95 


n 


125m'-|-'g 


58 d 


36 


i28Te 


31.69 


p 


127Sb 


3.84 d 


686, 474, 784 






n 


i27Te 


9.4 h 


418,380 








127mrpg 


109 d 


58 




33.8 


P 


129Sb 


4.4 h 


813,915,545 








129msb 


17.7 m 


760, 658, 434 






n 


129jg 


1.16h 


460, 28 








129mrpg 


33.6 d 


696 



Te Evidence for the existence of the second forbidden electron capture of this nu- 
cleus has been reported by various authors. An experiment recently carried out 
by the Milano group 11411 has shown that these evidences could be simply mim- 
icked by other background effects and that this decay has actually never been 
found. The reported values of the half-lifetime for this decay have therefore to 
be still considered only as lower limits. The discovery of this decay in CUORE 
would have relevant implications in nuclear structure theory 11421 . 



Rare alpha decays 

^°^Bi A search for the possible decay of this naturally occurring isotope of Bismuth, 
commonly regarded as the heavier stable nucleus, was proposed in our previous 
letter of intent. Recently, however a bolometric experiment |143| has detected 
the a-decay of this nucleus with a lifetime of (1.9 ± 0.2) x 10^^ years. Better 
measurements and a better study of this decay could be carried out in CUORE. 

180^ This nucleus can energetically a-decay into ^^^Hf. Indication of this rare decay 
was obtained with conventional detectors 114411 but unambiguously detected for 
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the first time thanks to the simuhaneous measurement of phonon and light signals 
with the CRESST cryogenic detectors lfH5J (T1/2 = (1.8 ± 0.2) xlO^^ y, Q = 
(2516.4 ±1.1 (stat.) ±1.2 (sys.)) keV ). This and other rare nuclear alpha decays 
could be investigated with great accuracy in CUORE. 
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5 Principles of operation of Te02 bolometers 

In a very naive approach a thermal detector is a sensitive calorimeter which measures the 
energy deposited by a single interacting particle through the corresponding temperature 
rise. 

This is accomplished by using suitable materials (dielectric crystals, superconduc- 
tors below the phase transition, etc.) and by running the detector at very low tem- 
peratures (usually below 100 mK) in a suitable cryostat (e.g. dilution refrigerators). 
According to the Debye law, in fact, the heat capacity of a single dielectric and diamag- 
netic crystal at low temperature is proportional to the ratio (T/Qd)^ {Qd is the Debye 
temperature) so that, providing that the temperature is extremely low, a small energy 
release in the crystal results in a measurable temperature rise. 

This temperature change can be measured by means of a proper thermal sensor. 

A low temperature detector (LTD or "bolometer") consists therefore of three main 
components (Fig|5](a)). 

Particle absorber : it is the sensitive mass of the device where the particles deposit 
their energy. The absorber material can be chosen quite freely, the only require- 
ments being, in fact, a low heat capacity and the capability to stand the cooling in 
vacuum. The absorber can therefore be easily realised with materials containing 
any kind of unstable isotopes and many interesting searches are therefore possible 
(e.g. beta decay spectroscopy, neutrinoless double beta decay and dark matter). 
The absorber mass can range from few micrograms to almost one kilogram. 

Temperature sensor : it represents a crucial element of the thermal detector needed 
to measure the temperature rise induced by the particle interaction. Among the 
variety of thermometers so far developed, doped semiconductor thermistors are 
preferred because of their ease of use and the reproducibility with which they can 
be produced in large quantities. 

Thermal link : it mainly depends on the way the massive crystal is held to the structure 
(heat sink). It influences the time response of the device. 

Absorber masses usually characterize LTD's applications and can range from few 
hundred micrograms (X-ray spectroscopy or neutrino mass direct measurements) to sev- 
eral hectograms (Dark Matter or DBD searches). 

5.1 Basic principles of thermal particle detectors 

The basic idea behind the thermal detector is very simple II146I . Let us assume that the 
detector has a heat capacity C (absorber) in thermal connection with a heat sink (kept at 
temperature Tq) through a thermal conductance G (fig. 13 a)). 
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Figure 5: Simplified model of a thermal detector (a) and schematic read-out of a resistive 
thermometer (b) 

Let T{t) be the temperature of the absorber as a function of time and let us assume 
that AT = |T(t) - To I < Tq for all times t, so that C and G can be treated as con- 
stants. Then an instantaneous deposition of an energy AE in the absorber gives rise to 
a temperature pulse whose time development is described by: 

AE / t\ C 
AT{t) = -^expi--j with r = - (24) 

As Tmax is given by AE / C, one immediately sees that, in order to have a significant 
percentage increase of the temperature (for an energy release in the keV-MeV range) it 
is necessary to have a very small heat capacity. As already stated above, this can be 
achieved only operating at very low temperatures. 

Very small specific heats can be attained by exploiting the strong decrease of the 
specific heat with temperature. The specific heat of a dielectric and diamagnetic crystal 
is in fact due only to lattice contributions [93|, with a temperature dependance following 
the Debye law 

C{T) = P (^] T<Qd (25) 



e 



D 



where (3 = 1994 J i^'^^mol^^ and Oz) is the Debye temperature of the crystal. 

In a very simplified model HI 471 in which all primary phonons undergo the thermal- 
ization process, an evaluation of the energy resolution can be obtained. In this case in 
fact the energy resolution of the detector is only limited by the thermodynamic fluctua- 
tions in the number of thermal phonons exchanged with the heat sink, which produces 
random fluctuations in the energy content of the absorber. Since the mean energy of 
thermal phonons is ~kT (k is the Boltzman constant), the average total number of ther- 
mal phonons N is proportional to C(T)/k and, assuming a Poisson statistics, this imply 
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a contribution to the energy resolution given by 



N = , AN = y/N ^ (Xe = AN ■ kT = J ^^kT = ^JkC(T)T^ (26) 

k \ k 

where ^ is a dimensionless factor (generally of the order of unity) whose actual value de- 
pends on the details of the temperature sensor, of the thermal link and of the temperature 
dependance of C(T). It should be stressed that, according to eq.|26l aE is independent 
of E, with impressive "theoretical" expectations even in the case of a massive crystal: 
a Te02 crystal with a mass of 1 kg operating at 10 mK could energy depositions of the 
order of few MeV with a resolution of about 20 eV. It should be also stressed, however, 
that this "thermodynamic limit" doesn't take into account other sources of energy fluc- 
tuations (e.g. metastable electron-hole states or long-lived non thermal phonons) which 
participate to the energy deposition process and whose contributions to the energy res- 
olution worsening should be therefore minimized. One should in addition take into 
account fluctuations in the gain due to temperature instabilities which could normally 
worsen the resolution. 



5.2 Semiconductor thermistors 

Thermistors are heavily doped high resistance semiconductors, with an impurity con- 
centration slightly below the metal-insulator transition II148II . Good thermistors require 
a higly homogeneous doping concentration (see se d7.3l) which is achieved, in the case 
of NTD (Neutron Transmutation Doped) thermistors, by means of thermal neutron ir- 
radiation throughout the entire volume |149| (sec. 17.21) . The electrical conductivity of 
these heavily doped semiconductors strongly depends on the temperature according to 
the hopping and variable range hopping (VRH) mechanisms II150I . In general the re- 
sistivity varies with temperature as 

P = Poe^^^ (27) 

where pq, Tq and 7 depends on the doping concentration. In the case of VRH, the value 
of 7 is 1/2. 

Thermistors are usually parameterized by their sensitivity A, defined as 

= 7(^y withi?(T) = p(T)^ = i?oexp(^^y (28) 

where I and s are the distance between the contacts and the cross section of the thermis- 
tor, respectively. The value of the sensitivity is usually in the range 2^10. 



A{T) 



d\nR 
d\nT 
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Figure 6: Typical load curve for a thermistor at T= 8 mK) (a) and resistance-power 
curves for a thermistor at different base temperatures (b). 

NTD germanium thermistors developed and built at LBNL have been used in all the 
large mass LTD's developed by the Milano group. They are made by exposing a ultra- 
pure germanium single crystal to a nuclear reactor neutron flux. They are characterized 
by an excellent reproducibility even if they are in part made by hand. The values of Rq, 
To and 7 must be experimentally measured for each thermistor (characterization pro- 
cess). This measurement is performed by coupling the thermistors to a low temperature 
heat sink using a high conductivity epoxy. The temperature Tb (base temperature) of 
the heat sink is then varied (15^ 50 mK) while a steady current / (bias current) flows 
trough the thermistors and a voltage V = IR appears across them. This produces a 
power dissipation which raises the temperature and acts back on the resistance, until an 
equilibrium is reached. This phenomenon makes the V — I relation deviate from linear- 
ity. This characteristic behaviour of bolometers is often referred to as "electrothermal 
feedback". The static resistance is simply the ratio V/I while the dynamic resistance 
(the only one actually involved in the preamplifier matching requirements) is the tan- 
gent at the V — I curve. By further increasing the bias current the dynamic resistance 
crosses the so called inversion point (where it vanishes) and becomes then negative. For 
semiconductor thermistors, a typical V — I curve, usually referred to as load curve, 
is represented in Fig|6l(a). By a combined fit to a set of load curves at different base 
temperatures, all the thermistor parameters are evaluated (see Fig|6l(b)) 

5.3 Detector operation 

As mentioned in section I5TT1 in first approximation, the thermal pulse produced by an 
energy release in the absorber is characterized by a very fast rise time and an exponential 
decay time determined by the physical characteristics of the individual detector. Such 
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a thermal pulse is then converted into an electrical one by means of the thermistor bias 
circuit and detector optimization is equivalent to determining the best configuration for 
a maximum pulse amplitude. Considering the basic circuit shown in Fig. |5tb), it's 
straightforward to obtain the relation between the maximum potential drop across the 
thermistor (Vb), the thermistor parameters and the deposited energy: 



where E is the deposited energy, A is the sensitivity given above in equation (ESI . 
and ATfe are the operating temperature and the temperature increase of the bolometer, 
C is the heat capacity of the absorber and P is the electrical power dissipated on the 
thermistor by the Joule effect. This expression vanishes both in the limit P ^ and 
P ^ 00. When the thermistor is operated as a detector sensor, a steady current / flows 
through it and, in a steady condition, its electrical and thermal parameters are described 
by a point on the detector V — I curve. This point is usually named working point, 
or optimum point when corresponding to the maximum detector response (maximum 
amplitude of the voltage pulse). The optimum point of each detector must be determined 
experimentally. This is accomplished by means of a fixed resistance {heater) directly 
glued onto the crystal through which it is possible to dissipate a known joule power into 
the crystal in a manner perfectly equivalent to a particle energy deposit. The optimum 
point correponds to the bias voltage for which the amplitude of the heater pulses is 
maximum. A typical bolometer response as a function of detector bias is shown in 



In fact, the actual optimum point is chosen by investigating the signal/noise ratio. 
Since the electronic sources of noise decrease with the resistance of the bolometer the 
best working point corresponds to a bias voltage slightly higher than the one correspond- 
ing to the optimum point (see Fig|7l). The use of the heater also plays a fundamental 
role in the stabilization of the bolometer response which changes with the temperature 
of the heat sink [11511 . 

5.4 Bolometer applications 

Since much information on the initial particle interaction is washed out by the thermal- 
ization process (position, momentum and interaction type), true cryogenic detectors are 
sensitive only to the total deposited energy which, however, can be measured very ac- 
curately (sec. 15.11) . In fact, the intrinsic energy resolution of such devices can be as low 
as a few eV. Main advantages of thermal detectors with respect to conventional ones 
are however represented by their sensitivity to low or non-ionizing events (Dark Matter 
(DM) searches) and material choice flexibility. Time response is sometimes poor (large 
mass bolometer pulse rise/decay times are of the order of few tens/hundreds msec or 




(29) 



FiglTl 
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Figure 7: left panel: evaluation of the optimum point of a bolometer; right panel: eval- 
uation of operation point which maximizes the signal/noise ratio. 

more) but, due to the very low counting rates, this is not a limitation at all for DBD 
(or DM) searches. It should be stressed that the flexibility in the choice of materials, 
besides allowing an experimental investigation on any DBD active isotope, means also 
the possibility of selecting materials characterized by a large natural isotopic abundance 
(a crucial parameter showing a linear dependence in the DBD lifetime sensitivity). Be- 
sides prohibitive costs, in fact, isotopic enrichment is not always easily attainable and 
the possibility to choose the DBD candidate element characterized by the highest natural 
isotopic abundance represents a big advantage. Because of the calorimetric approach, 
the thermal detectors detection efficiency is generally of the order of unity. Concerning 
energy resolution, up to now the actual detectors are still far from reaching the limits im- 
posed by the thermodynamic fluctuations of the system internal energy (eq.l^. Actual 
values are in the 1-10 keV range (FWHM) and seem to scale little with the detector mass 
(similar values were obtained in detectors whose masses scaled from few grams to few 
hectograms) while showing a quasi-linear dependence on the energy. Energy resolution 
is therefore still limited by other contributions (e.g. electronic and microphonic noise, 
vibrations of the system, etc.) and a good detector design (thermistor choice, mounting 
system, etc.), a low noise read-out electronics and an effective data acquisition and data 
analysis systems are of crucial importance for a well performing detector. So far, most 
of the R&D has been devoted to optimizing the detector performance. Although energy 
resolution doesn't seem to scale with detector mass, there is a limit to the single detector 
mass (imposed not only by the increase of the heat capacity but also by the technical 
feasibility of large perfect crystals). In the case of the Te02 , such a limit has been set so 
far to about one kilogram but larger masses could be possible in the near future (section 
[TSt . The only possibility for a detector mass in the ton range is therefore a large array 
of LTD's and the main issue is therefore reproducibility. This item has been faced and 
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solved few years ago by the Milano group with the realization of a Te02 20 bolome- 
ter array (340 g each) and, more recently, by the successful operation of CUORICINO 
(section lot . 

5.5 Detector model 

A rough approximation of the NTD thermistor voltage pulse amplitude AV5 produced 
by an energy deposition E into the Te02 absorber can be obtained by assuming RLoad >> 
Rb ineq. (Ell): 

AV, = V,-A-—^ = V,-A--— (30) 
Jfe o lb 

Using the heat capacity discussed earlier, a one M eV energy deposition into the Te02 ab- 
sorber results in AT;, = 7 X 10^^ K. Using equation (BUI) , this results in AV5 ~ 600 fiV, 
assuming typical electrical and thermal paramaters at the detector operating point. How- 
ever, the observed pulses are lower by a factor 2 to 3, therefore a more quantitative model 
is needed. 

A semi-quantitative Thermal Model (TM) has been developed which describes the 
detector in terms of three thermal nodes with finite heat capacities: the Te02 crystal ab- 
sorber, Ge thermistor lattice, and Ge thermistor electrons. For simplification, the lattice 
component of heat capacity is assumed to be negligible. The network describing the 
detector is shown in Fig. HI The two heat capacities, and four thermal conductances and 
the parasitic power from the vibrations and power dissipated in the thermistor electrons, 
are the model parameters. A typical energy deposition of 1 MeV in the crystal is used 
and it is assumed that it is instantaneously thermalized. 

A computer code was developed in Milan to determine the static behaviour. The re- 
quired input data are: the heat sink temperature, the temperature behaviour of the ther- 
mal conductances shown in Fig.|8](partially determined with ad hoc measurements), the 
parasitic power in the crystal, and a set of values of electrical Joule power dissipated in 
the thermistor electrons. The code determines a complete solution to the static problem, 
for each power level, and yields the temperature of the three thermal nodes, the resis- 
tance of the thermistor, and the current through it for each voltage-current (or power 
resistance) point of a load curve. When the static operation point is determined, the 
time evolution of the pulse is computed. The main aim of the TM is to have a reliable 
prediction of the detector figure of merit on the basis of the known thermal parameters. 
To address this question it is necessary to solve the non trivial problem of how to com- 
pare the performance of two detectors. A more sensitive thermistor would yield larger 
amplitude pulses, but they are also more sensitive to spurious noise, for example micro- 
phonics, crosstalk, and slow heat pulses from vibrations that excite the 1-5 Hz portion 
of the noise spectrum. A lower heat sink temperature yields larger pulses at the price of 
higher resistance and more spurious noise. 
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Figure 8: Network representing the detector model. 
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Recent experience shows that the best way to define a detector figure of merit is by 
fixing the heat sink temperature between 7 mK and 11 niK and by measuring a load 
curve. Using particle or heater pulses, the operation point on the load curve, correspond- 
ing to the highest signal is determined. For each heat sink temperature, one determines 
a pair of points, the thermistor resistance, R, and the signal amplitude, V, usually ex- 
pressed in microvolts per MeV, corresponding to the maximum pulse amplitude for a 
given heat sink temperature. One obtains a characteristic curve for that detector by plot- 
ting voltage vs resistance, which to a good approximation is a straight line on a log-log 
plot. This is called the Detector Merit Curve (DMC). Experience shows that a detector 
has constant signal-to-noise ratio along a DMC. No large differences in performance 
are noted for heat sink temperatures of 8 mK or 12 mK, because the signal level and 
the spurious noise offset one another. It is better to avoid the extremities of the DMC 
that lead to higher spurious noise or low signals. 

The best way to compare two different detectors is by using their DMCs. If the 
DMC of detector A lies systematically above that of detector B, the performance of 
detector A is superior. It provides the highest signals for a fixed thermistor resistance 
value. The main value of the TM is to predict the DMC given the temperature behaviour 
of the thermal parameters discussed above. 

As a test of the model, a typical detector of the MiDBD 20-element array was sim- 
ulated. These have different masses and geometries than the CUORE crystals. A DMC 
was then constructed using the model for this detector. It was then compared to the real 
DMC measured for the 20-detector array. A significant spread among the experimental 
points, for all 20-detectors was found. There were two reasons for this; intrinsic detec- 
tor differences and the presence of points not corresponding exactly to the maximum 
signal. Nevertheless the simulated DMC lies inside of the experimental points and has 
the correct slope (see Fig.|51). The detector simulations are not distinguishable from a 
typical detector of the array. 

The TM must correctly predict not only pulse amplitude, but also time constants. 
At a typical operating resistance (~ 100 MQ,) the predicted rise time (10%-90%) and 
decay time, are respectively 40 ms and 430 ms, and lie inside the distribution of these 
parameters for the 20-detector array. In Fig.[TOl the simulated pulse is compared with 
real pulses collected with four detectors, randomly chosen among the 20 elements. It is 
evident that the pulse-shape of the simulated detector is very similar to that of a typical 
detector of the array. From this it was concluded that the TM is capable of explaining 
semi-quantitative performance of the 20 - element array. 

The TM can be used in parallel with specific low temperature experimental tests, 
to define the optimum properties of the CUORE element, and the results thus obtained 
imply the following optimum parameters for the CUORE elements: 

1. Te02 absorber dimensions: 50 x 50 x 50 mm, with a mass of 750 g; 

2. NTD Thermistor # 31, Tq = 3.0 K, Rq = 1.5 n, dimensions 3x3x1 mm, with 
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Figure 9: Detector Merit Curve: simulation (solid curve) and experimental points for 
the 20 detector array. 
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Figure 10: Real pulses (full line) and simulated pulses (line+points). 
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a contact distance of 3 mm, and with a resistance of 50 MQ at T = 10 mK; 

3. absorber-thermistor coupling made with 9 epoxy spots, 0.5 to 0.8 mm in diameter, 
and 50 /im high; 

4. gold connecting wires 50 fim in diameter, and 15 mm in length; 

5. crystals mounted on PTFE blocks with a total PTFE-to-crystal contact area equal 
to about a few square mm. 

In spite of the reasonable agreement obtained between the TM predictions and the 
performance of the 340 g Te02 detectors, the comprehension of the detector operation 
principles is far to be complete. The proponents would like to stress in fact that this 
agreement was obtained with the average detector behaviour, while the TM seems to 
fail when one attempts to fully simulate an individual device. In particular, one observes 
a general trend according to which the simulated pulse amplitudes are smaller than 
the experimental ones. Inconsistency appears also in the pulse shape. This is true in 
particular for the CUORICINO 790 g detectors, as described below. 

In the CUORICINO case, the proponents have used an improved version of the code 
implementing the thermal model, capable not only of simulating but also of fitting to 
experimental data of the load curve of individual detectors. After fitting, the agree- 
ment between the simulated and the real curve is excellent. Of course the fit procedure 
changes the initial values of the thermal conductances, as measured in the specific exper- 
iments mentioned above. However, this modification is not drastic and can be justified 
with the inevitably not complete reproducibility of delicate manual operations like ther- 
mistor gluing and crystal insertion in the Teflon blocks. This proves that the static part 
of the model describes the detector satisfactorily. Problems show up when dealing with 
the dynamic behaviour. The simulated pulse, reconstructed using the thermal parame- 
ters extracted from the static behaviour of a specific detector and, in addition, the lattice 
heat capacity and the electron heat capacity (taken as a free parameter), is too small, 
typically a fraction 3 to 5, with respect to the experimental pulses collected for that spe- 
cific detector. The decay time is also in disagreement, being typically 3 to 5 times too 
fast. Note that this disagreement cannot be naively solved by postulating a crystal heat 
capacity lower than expected: in fact, this can solve the discrepancy between the pulse 
amplitudes, but at the price of making the disagreement between decay times even more 
severe. In addition, while it would not be surprising to find higher than expected spe- 
cific heats for the Te02 crystals at low temperatures, it is really difficult to justify lower 
values. The impression is that a model based on heat flows, which assumes thermal 
distribution of the involved excitations (phonons in the absorber, electrons in the ther- 
mistor), is not totally adequate for the description of the detector performance. Some 
attempts to modify the model by introducing athermal phonons directly absorbed in the 
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thermistor seem to improve our ability to describe the pulse shape and amplitude. How- 
ever, work is in progress and for the moment no conclusive statement can be made about 
the mechanism of pulse formation. 
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6 Previous results of MiDBD and CUORICINO 

The MiDBD experiment, carried out in Hall A of LNGS by the Milano group, was 
only the last of a series of experiments on neutrinoless DBD of ^^°Te, based on the 
bolometric technique. In fact, the improvements in the detector performance obtained 
over the course of years of measurements with detectors of constantly increasing mass, 
led to a large mass, low background experiment whose results finally demonstrated the 
competitiveness of the low temperature calorimeters. The idea of CUORE proceeds 
just from the successful results of MiDBD which can be therefore considered as a first 
step towards CUORE. The possibility to use the available MiDBD cryogenic system to 
test a single module of CUORE, thus checking the feasibility of this second generation 
experiment, finally lead to CUORICINO. Approved and funded in 2001, CUORICINO 
is actually the most sensitive experiment on neutrinoless DBD of ^^°Te, and is presently 
running in Hall A of LNGS, in the same cryostat used for MiDBD. 

6.1 LNGS cryogenic setups 

Two cryogenic setups were installed in the eighties at LNGS by the Milano group. The 
first, used in the past for the MiDBD experiment (Hall A) is now hosting the CUORI- 
CINO detector. The second (Hall C) has been dedicated to the research and development 
activities for CUORICINO and CUORE. Both cryogenic setups consist of dilution re- 
frigerators having powers of 1000 fiW and 200 fiW at 100 mK respectively. They are 
both equipped with a dedicated Helium liquefier, providing a substantial recovery of the 
Helium and preventing Helium contamination of the tunnel atmosphere, and are housed 
inside Faraday cages to suppress electromagnetic interference. A much larger experi- 
mental volume is available in the Hall A cryogenic setup. 

All the materials used for construction of MiDBD and CUORICINO were analysed 
to determine their radioactive contamination levels. These measurements were carried 
out by means of two large Ge detectors installed in the Gran Sasso underground Low 
Radioactivity Laboratory. The level of the radon contamination in the air of the Labo- 
ratory is continuously monitored. Both dilution refrigerators are equipped with heavy 
shields against environmental radioactivity. In particular, the Hall A dilution refriger- 
ator is shielded with two layers of lead of 10 cm minimum thickness each. The outer 
layer is of commercial low radioactivity lead, while the internal one is made with spe- 
cial lead with a ^^°Pb contamination of 16 ± 4 Bq/kg. The external lead shields are 
surrounded by an air-tight box flushed with fresh nitrogen from a dedicated evaporator 
to avoid radon contamination of the gas close to the cryostat. 

In order to shield the detectors against the unavoidable radioactive contamination of 
some fundamental components of the dilution refrigerator thick layers of Roman lead 
are placed inside the cryostat just around the detectors. A borated polyethylene neutron 
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Figure 1 1 : Scheme of the MiDBD detectors, 
shield (10 cm) was added in 2001 to the hall A cryostat. 
6.2 MiDBD experiment 

The main goal of MiDBD was the study of neutrinoless and two neutrino double beta 
decay of ^•^°Te to ^^°Xe using the thermal method. In addition, its results on the detector 
performance, and on the background, were demonstrated to be particularly useful for 
the CUORE and CUORICINO development. 

MiDBD consisted of an array of 20 crystals of Te02 of 3 x 3 x 6 cm^ each (340 g), for 
a total mass of 6.8 kg. It was operated in the Hall A dilution refrigerator until December 
2001. Sixteen MiDBD crystals were made of natural Te02. Of the remaining four, two 
were isotopically enriched at 82.3 % in ^^^Te and other two at 75.0 % in ^'^°Te. Mass 
spectrometer measurements determined that, due to the crystallization process, these 
enrichment factors were lower than the original ones (94.7 % and 92.8 %, respectively). 
The array, held in a copper tower-like frame mounted inside the inner vacuum chamber 
(IVC) of the dilution refrigerator, was maintained at a constant temperature of about 10 
mK. 

A Neutron Transmutation Doped (NTD) thermistor and a Si doped resistor were 
glued on the surface of each crystal. The former was used to read-out the thermal signal 
while the latter acted as a heater to generate a reference heat pulse in the crystal in order 
to monitor continuously the gain of each bolometer! 151 1. The electronic readout was 
accomplished with a room temperature low noise differential voltage-sensitive pream- 
plifier, followed by an amplifier and an antialiasing Bessel filter. The MiDBD detector 
was shielded with an internal layer of Roman lead of 1 cm minimum thickness (Fig. fTTTl . 
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Figure 12: An U+Th calibration, sum spectrum of the 20 detectors 
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Figure 13: Distribution of the single MiDBD detector energy resolutions (FWHM) on 
the^osTl 2615 keV line. 
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Figure 14: Background spectrum: the alpha particle region 
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In order to shield the detectors against the unavoidable radioactive contamination due 
to some fundamental components of the dilution refrigerator a layer of 10 cm Roman 
lead was placed above the tower of the array inside the cryostat. A similar layer, also 
inside the cryostat, was placed below the tower. The significant effort was carried out to 
check both the knowledge of the background sources and the performance of the new 
detector mounting system proposed for CUORICINO and CUORE. This led to a sec- 
ond run (MiDBD-II) in 2001 in which the detector structure was heavily modified (Fig. 
fTTl) . The MiDBD tower was completely rebuilt with a structure very similar to the one 
designed for CUORE/CUORICINO in which all the surfaces of the Te02 crystals and 
of the mounting structure were subject to a dedicated etching process aiming to reduce 
radioactive contaminants. On this occasion, due to the more compact structure of the 
tower, an enlargement of the internal Roman lead shields was also realized. An external 
neutron shield (10 cm of borated polyethylene) was added in June 2001. MiDBD data 
collection was divided in runs of about 15-20 days. At the beginning and at the end of 
each run the detector was calibrated with an exposure, lasting about 3 days, to a mixed 
U and Th source placed at the two opposite side of the 20 crystal tower, just in contact 
with the external cryostat thermal shields. 

The excellent reproducibility of the MiDBD detectors is shown in Fig. fT2l where 
the total calibration spectrum obtained by summing the spectra of the 20 detectors is 
presented. The overall energy resolution in the total calibration spectrum and in the 
total background spectrum are reported in Table |9l The single detector energy resolution 
distributions at the ^osjl 2615 keV 7-line are shown in Fig. [131 

The two MiDBD runs totalled ~3 1,508 (MiDBD-I) and ~5,690 (MiDBD-II) kg x hours 
of effective running time, respectively. Sum spectra were obtained both with no anti- 
coincidence cut and by operating each detector in anticoincidence with all the others 
("anticoincidence cut"). In all these spectra the main lines due to the natural activity 
of the 232Th and ^ssy chains, of ^^K and the lines at 1173 and at 1332 keV due to 
cosmogenic ^°Co were observed. 

Table 9: FWHM resolutions (keV) in a calibration spectrum and in the 84000 hours x 
crystal background spectrum. 
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Figure 15: Scatter plot of coincident events between all possible MiDBD detector pairs. 
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Figure 16: MiDBD-I (continuous line) and MiDBD-11 (dashed line)background spectra 
after anticoincidence cut. 

The sum spectra of the two runs obtained after applying the anticoincidence cut 
are shown in Fig. [161 A clear reduction in both continuum and 7/« peaks intensities 
is observed. This proves the effectiveness of the changes applied in MiDBD-II and 
provide information on the localization of the main background sources (see Table fTUI). 
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Table 10: Background reduction factors measured between the two runs of MiDBD. 

A background analysis, supported by the Monte Carlo simulation of the detector, 
allows us to disentangle the different sources responsible of the MiDBD background in 
the two runs and to understand which of the changes applied in MiDBD-II are respon- 
sible of the observed reductions. The conclusions thus obtained are the following: 
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• the 7 peaks above ~500 keV of the U and Th chains and the two Co 7 peaks at 
1.173 and 1.332 MeV are mainly due to sources outside the detector, these peaks 
are substantially reduced by the increased shield of Roman lead; 

• the 7 peak is probably due to multiple sources, its reduction being due both 
to the surface treatements and to the increased shield of Roman lead; 

• the a peaks of the ^^®U and ^^^Th chains, clearly visible in the background spec- 
trum above 4 MeV, are mainly due to surface contamination of the crystals, their 
reduction has to be accounted for by the surface treatment of the crystals; 

• ^^^U and ^^^Th 7 peaks below ~ 500 keV are attributed to a surface contamination 
of the copper mounting structure; 

• low intensities 7 lines of rather short living isotopes are explained by the cosmo- 
genic activation of both copper and tellurium. 

According to our Monte Carlo simulation of the MiDBD detectors, the foregoing 
sources can completely account for the observed background (peaks and continuum) 

measured in the two runs. 

A semi-quantitative analysis was carried out considering only MiDBD-I data for 
which the available statistics were more significant. In particular, the sources responsi- 
ble of the background measured in the P(3{0i') region were identified as follows: 

1. ^^^Th contaminations external to the detector structure (only the 7 emission from 
this chain can contribute to the background, mainly throught the 2615 keV line of 

2. ^^^U and ^^^Th surface contaminations of the crystals (a, (3 and 7 emissions from 
these chains contribute to the background) 

3. ^^^U and ^^^Th surface contamination of the copper mounting structure {a, (3 and 
7 emissions from these chains contribute to the background) 

The first source is identified from the clear 2615 keV ^osj^ ^jj^g gggj^ jj^g spectra 
of both MiDBD runs. Its localization is based on the analysis of the relative intensi- 
ties of the 7 lines (of different energies) belonging to the ^^^Th chain. Environmental 
7 flux in the underground laboratory (measured with a Ge detector just close to our 
experimental set-up) and the measured contamination of the external lead shields lead 
us to rule out these two possible sources for the ^osj^ ^jj^g -^he contamination seems 
therefore to be in the cryostat, and it could reasonably be the sum of the contamination 
of the various thermal shields that surround our detector. In order to extrapolate the 
contribution of this source to the background in the PP{Ou) region, we analyzed the cal- 
ibration spectra obtained with an external ^^^Th source located just outside the cryostat. 
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Figure 17: Scheme of a surface a decay. 

The shape of the calibration spectrum is in fact in good agreement with the shape of 
the simulated spectra for radioactive ^'^^Th contaminations distributed over the cryostat 
thermal shields. The observed ratio between the area of 20871 line at 2615 keV and the 
integral in the /3/3(Oz/) region is 14± 2. When compared with the 2615 keV rate ob- 
served in the background spectrum (129 counts/kg/y), this translates into a contribution 
of 0.1 1 counts/keV/kg/y which is ~20% of the total /3/5(Oz/) background rate. 

The second source can be easily recognized on the basis of the a peaks visible in 
the background spectra. Its location on the surface is deduced from the low energy tails 
of the peaks and by the coincidence scatter plots. These features exclude in fact the 
possibility of a bulk origin. Indeed the peak energies correspond to the transition en- 
ergies of the decays implying that both the a particle and the nuclear recoil contribute 
to the signal. The decay can only be "contained" in the crystal and the contamination 
producing the peaks is localized at a sufficient depth of the crystal surface. The low 
energy tail can be explained in terms of the fraction of the decays occurring very near 
to the surface, where part of the a or of the recoil energy is lost outside the crystal. 
Since the "lost" energy can be however observed by a second detector faced to the one 
in which the decay occurred (fig. (TTl. a straight line with a negative slope equal to -1 
is expected when plotting coincident events in different faced detectors as is shown in 
Fig. [151 This kind of analysis not only yields a further proof that a peaks are due to a 
surface contaminations of the crystals, but also proves that this contamination produces 
a continuum background extending far below the a decay transition energy. In particu- 
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lar this source can partially - but not completely - account for the rather flat continuous 
background measured between 3 and 4 MeV, a region where no known relevant 7 nor 
a lines from natural and cosmogenic radioactivity are present. By the analysis of the 
coincident events from a particles we obtained a rate Rd=0.13 counts/keV/kg/y, in the 
3-4 Mev interval, for the events with multiplicity 2 (i.e. events in which 2 crystals are 
simultaneously hit). The only possible origin for these events, is radioactive decay pro- 
duced in the fraction of crystal surface facing other crystals (S Active)- On the other hand, 
when the decay occurs in the fraction of crystal surface (S/nert) which faces inert mate- 
rials (e.g. copper), this contributes to the background remaining in the anticoincidence 
spectrum. By simply considering the geometric evaluation of the ratio Sjnertl^ Active 
(average Monte Carlo evaluation) we can obtain from R/j the expected contribution to 
the anticoincidence spectrum in the 3-4 MeV region: 0.1 6± 0.05 counts/keV/kg/y. In 
order to extrapolate this rate back to the /5/5(Oi/) region we need to know the shape of 
the background spectrum. In particular, various estimates of the ratio R3_4/R^^ (where 
R* are the counting rates in the energy regions of interest after the anticoincidence cut) 
were obtained on the basis of Monte Carlo simulations of crystal surface contamina- 
tions of different depths and density profiles. Unfortunately the ratio depends critically 
on these parameters and the average Monte Carlo evaluation of the ratio has therefore 
a large error. The extrapolated contribution to the PPiOu) region background is 0.37± 
0.1 counts/keV/kg/y corresponding to ~43± 10% of the measured rate. 

Concerning the third source, the prove of its relevance is less direct and is mainly 
based on the observation that the continuous background measured between 3 and 4 
MeV cannot be completely explained by the surface contamination of the crystals. 

Other possible sources for this background were taken into account and excluded as 
relevant for MiDBD: neutrons (ruled out by comparing the background rates measured 
with and without the neutron shield), muon induced neutrons (ruled out by simulation 
of the production rate of neutrons in the experimental set up), contamination of other 
components close to the detectors (such as the Teflon spacers placed between the copper 
frame and the crystals, the NTD thermistors, the Si heaters, for which the measured 
contamination is too low to account for the observed background). Finally, some tests 
done with dedicated detectors in the hall C set-up proved that the etching of the copper 
structure in a clean environment yields a reduction of the background in the 3-4 MeV 
region. This fact, and the agreement of Monte Carlo simulation with the measured 
spectra, indicate that ^^^U or ^^^Th surface contamination of the copper should be the 
missing background source. 

As in the case of the second background source mentioned above, the contributions 
of the third source to the background rate can be obtained by extrapolating the rate 
in the 3-4 MeV to the PP^Ou) region. The 3-4 MeV region rate (after the anticoin- 
cidence cut) can be evaluated by subtracting the crystal surface contribution from the 
total one: 0.33± 0.1 counts/keV/kg/y. Once again the extrapolation to the /9/3(0i/) re- 
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Figure 18: Total spectrum (in anticoincidence) in the region of neutrinoless DBD ob- 
tained with the twenty crystal array. The solid curves represent the best fit (lowest curve) 
and the 68 % and 90 % C.L. excluded signals. 

gion is model dependent (even if, according to our Monte Carlo evaluations, less so 
in this case) and the extrapolated contribution the /3/3(Oz/) region background is 0.2± 
0.1 counts/keV/kg/y (~23± 11%). The results of this analysis are summarized in ta- 
bleCT] 



Table 1 1 : Estimate of the weight of the different sources responsible for the background 
measured in the MiDBD-I run. 



Source 


208T1 pp{Ou) region 


3-4 MeV region 


Te02 ^^^U and ^^^Th surface contamination 


43± 13% 


33± 10% 


Cu ^^^U and ^'^^Th surface contamination 


~10% 23±11% 


67± 20% 


^^^Th contamination of cryostat Cu shields 


~90% 20± 5% 





The sum of the spectra obtained in anticoincidence in the two runs in the region 
above 2000 keV is shown in Fig. [TSl It corresponds to ~3.56 kg xyear of Te02 and 
to ~0.98 kg xyear of ^'^°Te . The clear peaks corresponding to the lines at 2104 keV 
(single escape of the 2615 keV ^o^Tl line), at 2204 keV (^^^Bi) and at 2615 keV (^o^Tl 
), confirm the reproducibility of the array during both runs. 
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Table 12: Half lifetime limits (90 % C.L.) on lepton violating and conserving channels 
deduced from the MiDBD data analysis. Eq is the energy analyzed to obtain T1/2, while 
a.c. indicates that the anticoincidence spectrum was used. 



Isotope 


Transition 


Cuts 


Eo 


Efficiency 


Tl/2 










(keV) 


(%) 


(years) 


130jg 


Oz/ : 0+ - 




a.c. 


2528.8 


84.5 


> 2.1 X 10^3 




Oz/ : 0+* - 


-^0+ 


none 


1992.8 


7.9 


> 3.1 X 10^2 


130jg 


Oz/ : 0+* - 


2+ 


none 


1992.8 


37.5 


> 1.4 X 10^3 


130jg 


2z/ : 0+ - 


-.0+ 


a.c. 






(6.1 ± 1.4+2|) X 10^0 


130jg 


IX : 0+ - 


-^0+ 


a.c. 






> 2.2 X 10^1 


130jg 


2x : 0+ - 


-^0+ 


a.c. 






> 0.9 X 10^1 


i28Te 


Oz/ : 0+ - 


-.0+ 


a.c. 


867.2 


97.9 


> 1.1 X 10^3 



No peak appeared in the region of /5/3(Oz/) of ^^^Te, where the rates were, respec- 
tively, of 0.59 ± 0.06 and 0.33 ±0.11 counts/keV/kg/y, when operated in anticoin- 
cidence. Also, no peak appeared at the energies corresponding to neutrinoless DBD 
of ^3°Te to excited levels of ^^^Xe, nor at the energy of 867 keV corresponding to 
/3/3(Oz/) of ^^^Te. Fit parameters and 90 % C.L. limits for the various decay pro- 
cesses (Tab. fT2ll were evaluated with a maximum likelihood procedure I1152J. Similar 
results were however obtained following the approach proposed by G.J. Feldman and 
R.D. Cousins lfT54t . 

6.3 CUORICINO experiment 

CUORICINO (Fig. fT9b is a tower made by 13 planes containing 62 crystals of Te02 . 
44 of them are cubes of 5 cm side while the dimension of the others is 3x3x6 cm'^ . 
All crystals are made with natural tellurite, apart two 3x3x6 cm'^ crystals which are 
enriched in ^^^Te and two others enriched in ^'^^Te with isotopic abundance of 82.3 % 
and 75 %, respectively. The total mass of CUORICINO is 40.7 kg, the largest by more 
than an order of magnitude with respect to any cryogenic detector. More details on 
the preparation of the crystals and on the mechanical structure of the array is reported 
elsewhere II153I . In order to shield against the radioactive contaminants from the mate- 
rials of the refrigerator, a 10 cm layer of Roman lead with ^^°Pb activity of <4 mBq/kg 
II151I is inserted inside the cryostat immediately above the CUORICINO tower. A 1.2 
cm lateral layer of the same lead is framed around the array to reduce the activity of 
the thermal shields. The cryostat is externally shielded by two layers of Lead of 10 cm 
minimal thickness each. While the outer is made by common Lead, the inner one has 
a ^^°Pb activity of 16 ± 4 Bq/kg. An additional layer of 2 cm of electrolytic Copper 
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is provided by the cryostat thermal shields. The background due to environmental neu- 
trons is reduced by a layer of Borated Polyethylene of 10 cm minimum thickness. The 
refrigerator operates inside a Plexiglass anti-radon box flushed with clean N2 and inside 
a Faraday cage to reduce electromagnetic interferences. Thermal pulses are recorded 
by means of Neutron Transmutation Doped (NTD) Ge thermistors thermally coupled 
to each crystal. Calibration and stabilization of each bolometer is performed by means 
of a resistor attached to each absorber and acting as an heater. The front-end electron- 
ics for all the 3x3x6 cm'^ and for 20 of the 5x5x5 cm"^ detectors are located at room 
temperature. In the so called cold electronics, applied to the remaining 24 detectors, the 
preamplifier is located in a box at ~ 100 Kelvin near the detector to reduce the noise due 
to microphonics, which would be very dangerous when searching for WIMPS. More 
details on the read-out electronics and DAQ are reported in [ll53i . CUORICINO is 
operated at a temperature of ~8 mK with a spread of ~one mK. A routine energy cali- 
bration is performed before and after each subset of runs, which lasts about two weeks, 
by exposing the array to two thoriated tungsten wires inserted in immediate contact with 
the refrigerator. All runs where the average difference between the initial and final cali- 
bration is larger than the experimental error in the evaluation of the peak position were 
discarded. During the first cool down 12 of the 5x5x5 cm"^ and one of the 3x3x6 
cm"^ crystals were lost, due to disconnection of a few of the thermalizers which allow 
the transmission of the electric signals from the detectors to room temperature 115311 . 
Since the active mass was of ~30 kg and the energy resolution excellent, we decided to 
continue running for a few months before warming up the array and attack the problem, 
which has been now fully solved. The data presented here refers to this first run and to 
the new data obtained in about two months with a second run where the contacts of only 
two bolometers were lost. The sum of the spectra of the 5x5x5 cm^ and 3x3x6 cm^ 
crystals in the region of the neutrinoless DBD is shown in Fig. IJU The background at 
the energy of neutrinoless DBD (i.e. 2510-2580 keV) is of 0.21 ± 0.02 and 0.19 ± 0.06 
counts/kg/keV/y for the 5x5x5 cm^ and 3x3x6 cm^ crystals respectively. No evidence 
is found at the energy expected for neutrinoless DBD of ^^°Te . By appling a maximum 
likelihood procedure 111521 1156ll we obtain a 90% C.L. lower limit of 1.0 x lO^'' years 
on the half-lifetime for PP{Ou) of this nucleus. The unified approach of G.I.Feldam and 
R.D.Cousins II1541I155 1 leads to a similar result. The upper bounds on the effective neu- 
trino mass that can be extracted from our result depend strongly on the values adopted 
for the nuclear matrix elements. As in our previous paper 1 1531 we consider all vailable 
theoretical calculations Il55l l2lfni58l apart those based on the shell model which are not 
considered as valid for heavy nuclei |[40|. We have also not considered for the moment 
the calculation by Rodin et al. since in the case of ^•^'^Te they are based on the not 
yet established value for the two neutrino DBD lifetime. In particular the adopted value 
of 2.7 xlO^^ years II158I is the largest among all geochemical ones L2, l.,A8|. Taking 
into account the above mentioned uncertainties, our lower limit leads to a constraint on 
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the effective neutrino mass ranging from 0.26 to 1.4 eV and partially covers the mass 
span of 0.1 to 0.9 eV indicated by H.V. Klapdor-Kleingrothaus et.al. Iil59l . 

6.3.1 CUORICINO background measurements 

The total collected statistics are therefore of ~4.76 kg xyear for 5x5x5 cm^ crystals 
and ~0.53 kg xyear for the 3x3x6 cm'^ ones. The corresponding background spectra 
are shown in Fig. |23l The gamma lines due to ^°Co, ''°K and of the ^'^^U and ^^^Th 
chains are clearly visible. 

These lines due to contamination of the apparatus, are not visible in the spectra of 
the single detectors: they appear after summing the 29 detectors, and are a good check 
of the calibration and stability of the detectors during the background measurement. 
Also visible are the gamma lines due to Te activation (^^^Te, ^^^'"Te, ^^^'"Te, ^^smjg 
and ^^''™Te) and those due to Cu activation (^^Co, ^^Co, ^°Co and ^'^Mn) by cosmic ray 
neutrons while above ground. 

The FWHM resolution of 5x5x5 cm'^ detectors at low energy, as evaluated on the 
122 keV gamma line of ^^Co, is ~2.8 keV. The ^osj^ gamma line at 2615 keV - clearly 
visible in the background sum spectrum - is used to evaluate the energy resolution in the 
region of double beta decay, the FWHM is 8.7 keV. 

The 3x3x6 cm"^ and 5x5x5 cm'^ crystal background spectra (sum of all the antico- 
incidence spectra of the detectors) are compared in Fig.|211and Fig.E51 

Here the statistical accuracy is much less, nevertheless the gamma lines, not visible 
in the single detector spectra, are clearly visible in the background sum spectrum. The 
FWHM resolution at low energy, measured on the 122 keV gamma line of ^^Co, is ~ 1.5 
keV. The FWHM resolution on the ^osjl gamma line at 2615 keV is evaluated to be of 
about 12 keV (this value has a large error due to the poor statistical significance of the 
peak). 

6.3.2 CUORICINO background analysis 

Because of the low number of counts, any evaluation of the possible background sources 
necessarily has a large statistical uncertainty. Some preliminary considerations are how- 
ever possible. They are based on the analysis of the sum spectra of the operating small 
and large crystals (the count rates of the single detector spectra are too low to consider 
them separately). In Fig.|24land Fig.|25]a comparison between the background spectra 
of the large and small crystals is shown. General spectral shapes and counting rates 
(when normalised to the mass of the crystals) are quite similar. On the other hand the 
intensities of the gamma lines do not show a clear behaviour: the ^^^Bi 1764 keV (the 
only clearly visible one of the U chain) and the ^osj^ 2615 keV lines (the only clearly 
visible one of the Th chain) seem to scale with the efficiency of the detectors (according 
to the results of a Monte Carlo analysis, the ratio of the detection efficiency of 5x5x5 
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Figure 20: Summed calibration spectrum (^^^Th source just outside the crysotat) from 
all the operating 5x5x5 cm^ and 3x3x6 cm^ crystals. 




Figure 21: Distribution of the single CUORICINO detector energy responses normal- 
ized to 1 kg of Te02. 
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Figure 22: Distribution of the single CUORICINO detector energy resolutions (FWHM) 
at the 208x1 2615 keV line. 




Figure 23: Summed background spectra from the operating 5x5x5 cm'^ and (natural 
abundance) 3x3x6 cm^ crystals. 
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Figure 24: Comparison between the background of the 5x5x5 cm^ crystals and that of 
the natural 3x3x6 cm'^ crystals in the gamma region. 

cm^ and 3x3x6 cm^ crystals is ~3). The ''^K line at 1460 keV has the same intensity 
per unit mass on the two kind of detectors. The ^°Co lines at 1 173 and 1332 keV, on the 
contrary show a higher intensity per unit mass in small crystals. Moreover they undergo 
a larger reduction by the anticoincidence cuts with respect to other 7-lines. 

In Fig. EH a comparison between the background spectra of CUORICINO crystals 
and the background measured in MiDBD-II (20 crystal array second run) is shown. 

In Table [13] the counting rates in different energy regions are reported. It is clearly 
evident that an improvement was obtained in the 2-3 MeV region (a reduction of ~30% 
in the counting rates per unit mass either of the small and of the large crystals) despite 
the reduced internal lead shield used in CUORICINO (~1 .2 cm thickness of roman lead) 
with respect to MiDBD-II (~3 cm thickness of roman lead). An even better reduction 
was obtained in the 3-4 MeV region where background is mainly dominated by surface 
contamination contributions. 

With higher statistical accuracy it will be possible to analyze the reasons for this 
variation of the counting rates and also for the differences between small and large 
crystals. To do that, it is indeed necessary to disentangle the different sources of back- 
ground by considering the counting rates in smaller regions as that between the two 
gamma peaks at 2448 and 2615 keV (^'^'^U and ^^^Th chains respectively) - where we 
should have contributions by both degraded a's from surface contaminations and ^os 
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Figure 25: Comparison between the background of the 5x5x5 cm^ crystals and that of 
the natural 3x3x6 cm'^ crystals in the alfa region. 



I m 




im 2000 3000 4000 5000 6000 7000 8000 1000 2000 3000 4000 5000 6000 7000 SOOO 



Energy [keV] Energy [keV] 



Figure 26: Comparison between the background spectra of CUORICINO crystals and 
the background measured in MiDBD-II experiment. 
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Figure 27: Summed background spectrum (5.29 kgxy) from all the operating crystals 
in the region of neutrinoless double beta decay of ^'^^Te (Q-value=2528.8 keV). 
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counts/keV/kg/year 


1-2 MeV 


2-3 MeV 


3-4 MeV 


4-5 MeV 


MiDBD-II 
3x3x6 natural 
5x5x5 


3.21±0.08 
3.29±0.11 
3.25±0.09 


0.61±0.04 
0.38±0.04 
0.41±0.03 


0.29±0.02 
0.24±0.03 
0.17±0.02 


1.88±0.06 
0.78±0.05 
0.55±0.03 



Table 13: Counting rates per unit mass in MiDBD-II and in CUORICINfO. 



7's from bulk contaminations - and the region just above the ^^^T\ 2615 keV line where 
background sources should be limited to surface contaminations. Moreover, it should 
be stressed that the effect of disconnected detectors in reducing the efficiency of the an- 
ticoincidence cut (a fundamental tool to study and localize background sources) is not 
negligible. 



Table 14: Estimate of the relative contributions of the different sources responsible for 
the background measured in CUORICINO. 



Source 


208T1 f3(3{0u) region 


3-4 MeV region 


Te02 ^^^U and ^^^Th surface contamination 


10± 5% 


20± 10% 


Cu ^^®U and ^'^^Th surface contamination 


~15% 50± 20% 


80± 10% 


^^^Th contamination of cryostat Cu shields 


-85% 30± 10% 





By using the same background model resulting from the MiDBD analysis, and tak- 
ing into account that we changed both the detector structure and the internal set-up of 
the cryostat (the roman lead shield and the cryostat copper radiation shields) we evalu- 
ated its consistency for CUORICINO. Because of the poor statistical significance of the 
data so far collected, the analysis has to be considered as preliminary (it was in particu- 
lar limited to the 5x5x5 cm^ detectors for which the statistical accuracy is better) even 
if its results already tend to identify the most important sources of the CUORICINO 
background besides giving a first quantitative guess of their relevance. The percentage 
result of the measured background in the DBD region, is summarized in table [141 It 
should be noticed that CUORICINO data seem to indicate a shallower depth for the 
crystal surface contaminations relative to MiDBD. A better determination of the actual 
depth and density profile of this contamination will be available when a fully operating 
CUORICINO and a better statistical accuracy will allow a more significant analysis of 
the multiple events. 
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Figure 28: Comparison between spectra obtained with different linearization methods. 
The appearence of clear alpha structures with the power law metod is evident. 

6.3.3 Analysis of the CUORICINO background above 3 MeV 

The following analysis tries to identify the background sources responsible of the CUORI- 
CINO counting rate in the region above 3 MeV. The analisys is based on both the co- 
incidence (not used in the previous analyses) and the anticoincidence spectra collected 
with 5x5x5 cm^ detectors during the first run of CUORICINO and calibrated with a 
new, more reliable, techinque in the alpha region. The single spectra have been cali- 
brated and linearized in energy with a new technique that yields quite better results in 
the alpha region: both the gamma lines of the ^^^Th source measurements and the alpha 
line of ^^°Po clearly visible in the background spectra are used, assuming a power law 
dependence for the relashionship between pulse amplitude and particle energy. The an- 
ticoincidence sum spectrum of 5x5x5 cm^ detectors calibrated in this way is compared 
in Fig.|28lwith the analogous spectrum used in the previous analysis and calibrated only 
with gama lines (assuming here a polynomial relashionship between pulse amplitude 
and particle energy and extrapolating this calibration in the alpha region). 

As already mentioned in previous discussions, alpha and beta particles produced 
by environmental radioactivity can give important contributions to background above 3 
MeV. Because of their low range these particles can only be located either in the crystals 
or on the surfaces of the materials directly facing the crystals: the copper holder or, less 
likely, small components placed near the crystals or on the crystal surface such as the 
NTD thermistor, the Si heater, their gold wires, the glue and other smaller parts. Monte 
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Figure 29: Comparison between Monte Carlo and CUORICINO anticoincidence (top) 
and coincidence (bottom) spectra in the case of the Te02 crystal surface contaminations 
(A ~ 1 yum) specified in the figure. 



Carlo simulations of the coincidence and anticoincidence sum spectra of 5x5x5 cm'^ 
detectors produced by radioactive contaminations of the above mentioned elements have 
been obtained using our GEANT4 based code. These spectra have been compared with 
the measured spectra in order to identify the actual contaminations responsible of the 
measured background. 

As a first step, bulk and surface contaminations of the crystals have been consid- 
ered. Both these contaminations give rise to peaks (centered at the transition energy of 
the decay) in the anticoincidence spectrum, but in the case of bulk contaminations these 
peaks are gaussian and symmetric while in the case of surface contaminations the peaks 
show a low energy tail. On the contrary, the only possible sources contributing to the 
coincidence spectrum are the surface contaminations of the crystals. The general shape 
and the structures appearing in both background spectra strongly depend on the depth 
and density profile assumed for the contamination (an exponentially decreasing func- 
tion characterized by a depth A: p(x) = Ae^'*'^ yielded the best results). The results of 
the comparison between Monte Carlo simulations and measured spectra allow to prove 
that the alpha peaks of the anticoincidence spectrum of Cuoricino are due to a surface 
(and not bulk) ^^^Th and ^^^U contamination of the crystals with the only exception of 




Figure 30: Comparison between Monte Carlo and CUORICINO anticoincidence (top) 
and coincidence (bottom) spectra when a Th surface contamination of the detector cop- 
per holder (green line) is added (pink line) to the crystal conatminations considered in 

fig.ESl 

the ^^°Po line which has to be attributed to a bulk contamination of the crystals. The 
shape of the peaks in the anticoincidence spectrum and the shape of the coincidence 
spectrum is fully accounted for when a surface contamination depth of the order of 1 
jjim or less is assumed. As shown in Fig.|251the identified crystal contaminations yield 
a satisfactory explanation of the coincidence spectrum. Some contribution to the con- 
tinuous part of the anticoincidence spectrum seems however still missing. In order to 
explain this continuum, a further source, not contributing to the coincidence spectrum, 
has to be considered. Such a source can be a ^^^Th or ^^^U contamination on the copper 
surface, as shown in Fig. EOl The depth of the contamination should be of the order of 
~5 yum (a deeper contamination would produce too high gamma peaks while a thinner 
contamination would give rise to structure in the anticoincidence spectrum) with a total 
activity in the first 1 .5 /im of the order of 10^^ g/g either in ^^^U or ^'^^Th or both. 



6 PREVIOUS RESULTS OF MIDBD AND CUORICINO 



83 



This analysis shows the importance of continuing the CUORICINO background 
measurement in order to improve the statistical accuracy of the above described results. 
These results, on the other hand, show that with a substantial reduction of the copper 
and Te02 crystal surface contamination levels the sensitivity goal for CUORE (i.e. a 
background counting rate of the order of 0.01 counts/keV/kg/y or better) can be reached. 

6.3.4 Preliminary measurements of the Copper surface contaminations 

According to the suggestions of the LNGS Scientific Committee we have investigated 
various possibilities to identify and measure the surface contamination of the detector 
Copper holders which is expected to represent the main contribution to the CUORE 
background level in the neutrinoless double beta decay region. We would like to remind 
that the bulk contamination of Uranium and Thorium in copper has been found to be 
less than a part per triUion, while the analysis of the MIBETA and CUORICINO data 
indicates a surface contamination by two to three orders of magnitude larger. Degraded 
a and (3 particles from this source have been found to represent a substantial contri- 
bution of the background in the region of neutrinoless double beta decay. Two main 
measurement procedures were considered: 

• dissolution of different superficial layers of Copper in ultra pure acid and analysys 
of the resulting solution with Inductive Coupled Plasma Mass Spectroscopy (ICPMS) 

• laser ablation of a number of distributed spots on the Copper surface followed 
by ICPMS analysis of the dissolved material. We did that in Canada, but the 
sensitivity seem still low ( 10-8 g/g) 

Other approaches like Total Reflection X-Ray Fluorescence (TRXRF) and Secondary 
Ion Mass Spectroscopy (SIMS) are also going to be considered, even if they allow to 
investigate only the very superficial layer of the Copper samples. 

The most powerful method was found by us to be a High Resolution (HR) ICPMS 
with which preliminary results were obtained in collaboration by the Gran Sasso Labo- 
ratory (M.Balata and S.Nisi) and the Mass Spectroscopy group of the European Center 
of Ispra (P.Tricherini). The procedure has been the following: samples of Copper were 
attacked with ultrapure nitric acid with different attacking times in order to extract the 
material from different thickness. The corresponding solutions were then analyzed in 
ISPRA with a HR-ICPMS device characterized by a sensitivity up to 10-15 g/g of U 
and Th in the solution. The sensitivity in g/g of the dissolved Copper was lower by 
about three orders of magnitude due to the small copper concentration in the solution. 
The preliminary results of a measurement carried out on a sample of the same Copper 
used for the CUORICINO detector holders is shown in Fig. 2. It can be seen that the 
average concentration of Uranium and Thorium strongly decreases when considering 
copper layers of 1.48 and 7.5 micron respectively (only a limited number of dissolved 
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layers of different thickness could be analyzed in the preliminary measurement). This 
trend of surface contamination as a function of depth is in excellent agreement with that 
evaluated from our CUORICINO background model (comparison between Monte Carlo 
simulations and actual background spectra). 

6.3.5 CUORICINO vs. CUORE 

While being a self-consistent experiment, CUORICINO finds its place in the CUORE 
project as a test facility intended to verify the technical feasibility of CUORE. Indeed the 

good results obtained as far as the detector performances are concerned prove that the 
realization of CUORE could be more straightforward than expected. There is still some 
RSzD that is required to improve detector reproducibility in order to reduce the spread 
in the main parameters that characterise detector performances (as the pulse height, the 
energy resolution and the detector signal rise and decay times). Nevertheless the good 
results obtained in CUORICINO prove that the increase in the number of bolometers 
and in the total mass of the array do not substantially affect the quality of the experi- 
ment: the MiDBD experiment with 3 times fewer detectors, and a mass 6 times lower 
than CUORICINO, had a similar energy resolution in the background sum spectrum. 
For what concerns background evaluation, as stated above CUORICINO is an excellent 
tool to predict the expected background for CUORE, when the different geometrical 
structure of the two setups is taken into account. Indeed CUORE will have a tightly 
packed structure where the operation of the detectors in anticoincidence will allow a 
strong reduction of background. Moreover the lead shield designed for CUORE will be 
optimized in order to practically cancel the background coming from outside. Such an 
optimization was not possible in CUORICINO which had to be housed in an already 
existing cryostat characterized by a limited space. Indeed to mount CUORICINO in the 
same cryostat that was used for MiDBD, the internal Roman lead shield thickness was 
necessarily reduced thus increasing the background due to 7 rays from the refrigerator 
materials. This will obviously be not a problem for CUORE. The CUORICINO back- 
ground results, although preliminary, are very promising. They demonstrate both the 
effectiveness of the procedures applied during its preparation and our knowledge of its 
main background contributions. The detailed study of CUORICINO background (i.e. 
with better statistical accuracy) will give more insight into background sources. 
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7 CUORE project 

The CUORE detector will consist of an array of 988 Te02 bolometers arranged in a 
cylindrical configuration of 19 towers of 52 crystals each (Fig. I^TT). The principle of 
operation of these bolometers is now well understood (sec. |5l) and since Te02 is a di- 
electric and diamagnetic material, according to the Debye Law (l25t . the heat capacity of 
a single crystal at sufficiently low temperature can be so low that also a small energy re- 
lease in the crystal can result in a measurable temperature rise. This temperature change 
can be recorded with thermal sensors and in particular using Neutron Transmutation 
Doped (NTD) germanium thermistors. These devices were developed and produced at 
the Lawrence Berkeley National Laboratory (LBNL) and UC Berkeley Department of 
Material Science II160I . They have been made uniquely uniform in their response and 
sensitivity by neutron exposure control with neutron absorbing foils accompanying the 
germanium in the reactor II161II . 

The Te02 crystals are produced by the Shanghai Institute for Ceramics (SICCAS), 
China and they will be the source of 750 g Te02 crystals for CUORE II162II . 




Figure 31: The CUORE detector (left), one of the 19 towers (right). 

A single CUORE detector consists of a 5x5x5 cm^ single crystal of Te02 that 
acts both as a detector and source. The detectors will be supported in a copper frame 
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in the 19 tower configuration shown in Fig. |3T1 The frame, and dilution refrigerator 
mixing chamber to which it is thermally connected, forms the heat sink, while the PTFE 
(Polytetrafluoroethylene or TEFLON) stand-offs provide the thermal impedance which 
delays the re-cooling of the bolometers. The bolometers operate at ~ 10 mK. 

As described above (sec. 16.31) . the excellent performance of the 5x5x5 cm^ crys- 
tals to be used in CUORE has been already proved by CUORICINO ifT^fT^ which 
corresponds to a tower of length and number of floors slightly larger than in CUORE. 
We adopted this solution to make use of the full experimental volume available in the 
dilution refrigerator previously used for the MiDBD experiment. 

The CUORE detector then will be the final version of a three step development com- 
prising: MiDBD, CUORICINO, and finally CUORE. This, plus the fact that CUORE re- 
quires no isotopic enrichment, (the isotopic abundance of ^^°Te is 33.8%) puts CUORE 
ahead of other options of truly next generation /9/3(0z/) experiments. The technology, 
though novel, is developed and to a large degree proven (sec.|5]and|6l). 

7.1 Single CUORE detector 

A single CUORE detector is a 5x5x5 cm'^ single crystal of Te02 grown with ultra- 
pure Te02 powders and polished on the surfaces. Crystals of Te02 have a tetragonal 
structure, and are grown along the (001) axis. The two axes normal to this axis are 
crystallographically equivalent, a fact relevant to their use in the search for solar ax- 
ions discussed previously (sec i4.2.2'b . The surface hardness is not the same for all sides, 
which complicates the crystal polishing. We have shown that repeated thermal cycling 
does not damage the crystals, as in the cases of crystals of other tellurium compounds, 
or those of tellurium metal. The Debye temperature of the Te02 crystals was specially 
measured for the CUORE project as 232 K II164I . This differs from the previous value 
of 272 K in the literature [d65l . The specific heat of the Te02 crystals was found to fol- 
low the Debye law down to 60 mK; the heat capacity of the 750 g crystals is 2.3 x 10~^ 
J/K extrapolated down to 10 mK. 

The NTD thermistors are attached by epoxy to the crystal and are operated in the 
Variable Range Hopping (VRH) conduction regime with a Coulomb gap 111661 11671 
11681 . One of the most significant parameters characterizing the thermal response is the 
sensitivity A, already introduced in section|5]and defined in equation (l28t . For our NTD 
thermistors this parameter ranges between 7 and 10. The resistance behaviour follows 
the relation: 

R = Roexp{To/Ty; 7 = 1/2. (31) 

(see the analogous behavior for the resistivity in equation (1271 ). 

The VRH regime occurs in Ge when it is "doped" close to the Metal to Insulator 
Transition, which is ~ 6x 10^^ atoms / cm^. This is achieved by thermal neutron radi- 
ation in a nuclear reactor (see sectior l7?2b . The sensitivity parameter A, depends on the 
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neutron irradiation dose. Therefore, each thermistor must be characterized at operating 
temperatures, as described in Sec. 15.21 and, for Si thermistors, in Ref.| 169|. 

It is very important to optimize the neutron irradiation exposure and to make the 
exposures as uniform as possible. It is not possible to evaluate the thermistor mate- 
rial directly from the reactor because of the long half life of ^^Ge (11.43 days). A 
delay of several months is required to see if the Ge needs more exposure. To circum- 
vent this difficulty, the Ge material is accompanied by foils of metal with long-lived 
(n, 7) radioactive product nuclides. Accordingly, the neutron exposure of the Ge can be 
determined accurately, and uniformity of exposure is achieved. This technique was de- 
veloped recently by the Lawrence Berkeley National Laboratory group of the CUORE 
Collaboration II161I . Following the neutron exposure and radioactive decay period, the 
NTD germanium is first heat treated to repair the crystal structure then cut into 3x3x1 
mm'^ strips. The thermistors are glued to the Te02 crystal by 9 spots of Araldit rapid, 
Ciba Geigy (now Novartis) epoxy, of 0.4 to 0.7 mm deposited on the crystal surface 
by an array of pins. The height of each spot is 50 fim. This procedure was found to be 
reasonably reliable and reproducible in the MiDBD experiment 11701 . The heat conduc- 
tance of the epoxy spots was measured in Milan and the phenomenological relation was 
found to be ~2.6 x 10^^ (T[K])'^ watts per degree kelvin per spot. The stabilization of 
the response of bolometers is crucial because of the unavoidable small variations in the 
temperature of the heat bath that change the detector gain (and consequently deteriorates 
the energy resolution; fig.l^^. 

This problem is successfully addressed by means of a Joule heater glued on to each 
crystal. The heater is used to inject a uniform energy in the crystal and the thermal 
gain is monitored and corrected off-line (see Sec.fTSt. The heaters are Si chips with a 
heavily doped meander structure with a constant resistance between 50 to 100 kTl. They 
are manufactured by the ITC - IRST company in Trento, Italy. Electrical connections 
are made with two 50 fim diameter gold wires, ball bonded to metalized surfaces on 
the thermistor. The gold wires are crimped into a copper tube, which is inserted into a 
larger one forming the electrical connection, and avoiding low temperature solder which 
contains ^^°Pb and traces of other radioisotopes. The larger copper tube, ~14 mm long 
and 2 mm in diameter, is glued to the copper frame that supports the crystals. This 
tube is thermally connected to the frame but electrically insulated. The mounting of the 
Te02 crystals is crucial to detector performance, and must fulfill a number of sometimes 
contradictory criteria: 

1. the crystals must be rigidly secured to the frame to prevent power dissipation 
by friction caused by unavoidable vibrations, that can prevent the crystal from 
reaching the required temperature and can produce low frequency noise; 

2. the thermal conductance to the heat sink (copper frame) must be low enough to 
delay the re-cooling of the crystal, following a heat pulse, such that the pulse 
decay time (re-cooling time) is much longer than the rise time; 
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Figure 32: Effects of the system thermal instabilities on the worsening of the detec- 
tor energy resolution (above). Comparison between calibration spectra before (black 
curve: AEfg^^*'^ = 30 keV) and after (red curve: AEfg?^^^'^ = 6 keV) the stabilization 
procedure. 
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3. however, the heat conductance must be high enough to guarantee efficient cooling; 

4. the frame must compensate for the differences in thermal expansion coefficients 
of the various materials used; 

5. and finally, only materials selected for low radioactivity can be used. 

For CUORE, only two materials will be used, copper and PTFE; they can both be ob- 
tained with very low levels of radioactivity. Copper has a thermal conductivity and spe- 
cific heat high enough to be an ideal heat bath, and has excellent mechanical properties; 
it machines well, and has excellent tensil, tortional and compressional strengths. 

PTFE is used between the copper frame and the crystals. It has low heat conductance 
and low heat leak II171II . It compensates for the differences between coefficients of 
thermal expansion of copper and of Te02 . The heat conductance between the crystal 
and heat sink was measured in a similar configuration, but with a contact area larger by 
a factor of 2 lfT72, l. 

The measurements of the temperature dependence of specific heats and thermal con- 
ductivities is crucially important for developing a model that explains detector perfor- 
mance and that can be used to improve thedetector design (seclsj). 



7.2 The NTD thermistors 

The key to the success of CUORE is the production of large numbers of nearly iden- 
tical neutron-transmutation-doped (NTD) Ge thermistors. These are made by neutron 
transmutation doping of ultra-pure Ge in a nuclear reactor to obtain the proper character- 
istics of resistance and A, the variation of resistance with temperature. Melt-doped Ge 
crystals cannot achieve the necessary uniformity due to a variety of dopant segregation 
effects. The only technique available for producing such uniform doping is NTD. In 
typical applications, the neutron absorption probability for a 3 mm thick wafer of Ge is 
small, on the order of 3 %, leading to a very homogenous, uniform absorption process. 
The most important aspect of this process is that ™Ge transmutes into Ga, an acceptor, 
and ^''Ge transmutes into As, a donor, the primary active dopants in NTD Ge. In this 
process, one places the Ge in a nuclear reactor where the following reactions take place: 

™Ge (21%) +n ^ ^^Ge {ar = 3.43 ± 0.176, = 1.56) 

^^Ge ^ ^^Ga (ti/2 = 11 Aday) Acceptor (32) 



Ge (36%) +n ^ ^^Ge {ar = 0.51 ± 0.086, = 1.0 ± 0.26) 

^^Ge ^ ^^As (ti/2 = 83mm) Donor (33) 



^^Ge (7.4%) + n ^^Ge (ay = 0.16 ± 0.0146, an = 2.0 ± 0.356) 

^1/2 



^^Ge ^^5e {ti/2 = 38.8/ir) DouhleDonor (34) 
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where ctt and an refer to the thermal and epithermal neutron capture cross sections, 
respectively. 

Since the doping level of the Ge needs to be on the order of 1 x 10^'' atoms/cm'^, a 
very high flux reactor, such as that at the University of Missouri or the Massachusetts 
Institute of Technology, is necessary to do the doping in a reasonable time. Even more 
important is the stability of the neutron flux and the energy distribution as measured by 
the Cd ratio. A large Cd ratio (>200) is desirable. In addition, a minimal fast (>5 MeV) 
flux is necessary to eliminate contaminants such as '^H, ^^Zn, and ^^Ge produced by 
fast neutrons. These radioactive contaminants degrade the detector operation because 
of internal radioactive decay and therefore increased background. 

In order to determine both the thermal and epithermal {^r) neutron fluences 
that the Ge wafers are exposed to during the reactor irradiations, we developed a method 
using multiple metallic monitors that are irradiated with every set of Ge samples. For 
each monitor foil, the amount of radioactivity produced depends on the sum (ar^r + 
o'r^r)- Each monitor provides one equation in two unknowns, the solution to which 
can be depicted as a line on the graph of thermal flux vs. epithermal flux. For mul- 
tiple elements, multiple lines can be drawn on this same graph, and they will have a 
common intersection at the solution of values of thermal and epithermal flux that sat- 
isfy all conditions. Essentially, this is a graphical solution to several equations in two 
unknowns. Elements with significantly different ratios of ot/ctr provide significantly 
different slopes to the lines and a better determination of the intersection. Fe and Zr are 
a good combination since they have significantly different otIor and long half-lives. 

Once the thermal and epithermal fluences have been determined, the concentration 
of dopants produced in the germanium can be calculated. The significant quantity in- 
volved in thermistor performance is the net dopant concentration which is equal to the 
concentration of Ga atoms minus the concentration of As atoms minus twice the con- 
centration of Se atoms. Unfortunately, to measure the thermal performance of the ther- 
mistor, one needs to wait for the decay of the activation product ''^Ge (1 1.4 day), which 
requires approximately one year. Following this one year of decay, some of the NTD Ge 
material was made into thermistors and sent to the University of Florence, for testing 
at low temperatures. Samples were cooled to 30 mK and resistance vs. temperature 
measurements were made over approximately the 300-30 mK range. These curves were 
compared tos similar curves from other batches of NTD Ge as shown in figure l34l 

For the various batches we ran (NTD-33-38) the data can be combined to show the 
variation of the temperature performance vs the calculated net acceptor concentrations. 
These results are illustrated in Figure|35l The three points on the right side of the data set 
were the results of two irradiations, and the middle of these three points has properties 
nearly the same at those of NTD-31, currently being used in CUORICINO. 

The problem of producing large numbers of identical thermistors is quite compli- 
cated. These thermistors, nearly identical to NTD-31, were made with an iterative pro- 
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Thermal and EpiUiermal Flux Solution for multiple monitors 
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Figure 33: Graph depicting the solution for thermal and epithermal flux for a set of three 
monitors, each with different ratios of thermal and epithermal cross-sections. 
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Figure 34: Resistivity of various batches of NTD Germanium as a function of \I\/T. 
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Figure 35: Thermistor performance as a function of the calculated net acceptor concen- 
tration. The first six points (from the left) involved only one irradiation (and show a 
relatively large spread in deviations from the line) while the last three points are follow- 
ing the second irradiation and have a very small spread from the line. 
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cess whereby the second doping was calculated following a measurement of the first. 
This is a very inefficient process for producing large numbers of thermistors. Fortu- 
nately, the same method used to make the thermistors can be used to analyze the results. 
Neutron capture on ^^Ga produces ^^Ga (14 hour half life) and ''^As activates to ^^As 
(26.3 hour half life). Both these isotopes emit gamma rays of high enough energy and 
intensity to make precise measurements of both Ga and As doping levels. Selenium 
cannot be measured by Instrumental Neutron Activation Analysis (INAA) because the 
product nuclide, ''^Se is stable. Similarly, Secondary Ion Mass Spectroscopy (SIMS) 
and other mass spectroscopic techniques can be used to determine the dopant concen- 
trations on a much shorter time scale. Two, or even three, irradiations will be required 
to produce thermistors with properties close enough to be used in the large array of 
CUORE. These techniques will be further refined in order to allow faster and more 
uniform production of the thousands of thermistors required for CUORE. 

The plan to make these irradiations more predictable is to use the NTD facility at 
the MIT reactor with multiple irradiations on a single large batch of Ge wafers. Each 
wafer will have its own monitor, and each will be separately analyzed to determine the 
appropriate subsequent irradiation to bring them to the right temperature characteristics. 

Previous measurements carried out in the Gran Sasso Laboratory [ 173 1 have shown 
that the residual acrtivity of the NTD thermistors becomes fully tolerable in an experi- 
ment with thermal detectors already a few months after irradiation. 

7.3 Modular structure of the CUORE detector 

As stated earlier, the CUORE array comprises 988 Te02 bolometers, grouped in 247 
modules of 4 bolometers each. These are arranged in 19 towers of 52 crystals each, or a 
stack 13 modules high. The towers are assembled in a cylindrical structure as shown in 
Fig. El 

The structure of each four-detector module is shown in Fig.|36l The four crystals are 
held between two copper frames joined by copper columns. PTFE pieces are inserted 
between the copper and Te02 , as a heat impedance and to clamp the crystals. There is 
a 6 mm gap between crystals with no material between them. 

The four detectors are mechanically coupled; some of the PTFE blocks and springs 
act simultaneously on two crystals. Tests of this and a similar structure, but with 3x3x6 
cm'^ crystals, clearly demonstrate that the CUORE technology is viable. The signal to 
noise ratio was improved by replacing the 3 x 1.5 x 0.4 mm^ Ge thermistor by one with 
dimensions 3x3x1 mm^, which significantly improved the DMC shown in Fig.|37j and 
by utilizing cold electronics Fig.EHI A further noise reduction was finally obtained by 
suspending the detector by means of a stainless steel spring. 

Reproducibility was tested on the MiDBD array, which had a sufficient number of 
detectors (20) operating simultaneously. The reproducibility is satisfactory as shown in 
the load curves in Fig.|39l Four-detector modules of 50x50x50 mm^ crystals were also 
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Figure 36: A four detector module. 
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Figure 37: Improvement in the Detector Merit Curve by increasing the thermistor size: 
3x3x1 mm'^ thermistors (stars) vs. 3 x 1.5x0.4 mm^ (diamonds). 
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Figure 38: Noise Power Spectra without damping suspension (1), with damping sus- 
pension only (2) and with dumping suspension + cold electronics (3). 

successfully tested, with pulse amplitudes spanning an interval from 50 to 150 /iV/MeV 
at ~100 operation point, in agreement with the thermal model. The full width at 
half maximum FWHM is ~1 keV for low energy gamma peaks, and ~5 to 10 keV at 
2.6 MeV. A stack of thirteen, 4-detector modules, will be connected by two vertical 
copper bars. The details of the wiring of the towers is still under development; however, 
the current plan is for 50 fim diameter twisted pairs of nylon coated constantane wire 
running along the vertical copper bars. On the top of each tower there would be 52 
twisted pairs for the thermistors and four for the heaters. Two more twisted pairs will be 
required for diagnostic thermometers, one on top and one on the bottom of the tower. 

The CUORICINO detector presently taking data at LNGS consists of one of the 52 
detector towers with the exception of the two planes containing the 3x3x6 cm^ crystals 
discussed earlier. CUORE is an array of 19 towers, each of which will be suspended 
independently from a large cylindrical copper plate, thermally connected to the mixing 
chamber of the dilution refrigerator (DR). In this manner, tests of CUORICINO will be 
tests of the CUORE array. The array will be suspended from the mixing chamber by a 
vertical spring to decouple the detector from the vibrations of the dilution refrigerator. 
Tests were performed with piezo-electric accelerometers and vibrators to optimize the 
suspension (see Fig.l38t. As a result a pendulum mechanical suspension was designed 
for the entire CUORE array. 
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Figure 39: Fifteen thermistor load curves for the 20 crystal array. 

The heavy shielding of CUORE will present non-trivial problems in the energy cal- 
ibration with radioactive sources. The presently considered option is to use radioactive 
metal wires encapsulated in PTFE tubes. The tubes would be placed vertically between 
the towers, and free to slide along a fixed tube to the appropriate calibration point. A 
vacuum tight sliding seal will allow the part of the tube containing wire to be inserted 
and extracted from the cryostat. 

7.4 Location 

CUORE will be located in the underground hall A of Laboratori Nazionali del Gran 
Sasso (LNGS, L'Aquila - Italy) at a depth of 3400 m.w.e. where the muon flux is 
reduced to ~3x 10^^ /i/cm~^/s and the neutron flux to ~10^^ n/cm^^/s. After the ap- 
proval of CUORE by the Gran Sasso Scientific Committee its final location just near the 
CRESST installation has been decided as shown in Fig.l40l 
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Figure 40: Final location of CUORE in hall A of the underground Laboratori Nazionali 
del Gran Sasso. 
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8 Cryogenic system 

The CUORE bolometers will operate at temperatures between 7 and 10 mK. This will 
require an extremely powerful dilution refrigerator (DR). At these temperatures, the 
cooling power of a DR varies approximately as T^. Estimates were made of the parasitic 
power the detector and DR would receive from: heat transfer of the residual helium gas 
in the inner vacuum chamber (IVC), power radiated from the 50 mK shield facing the 
detector, and from vibrational energy (microphonic noise). The estimated value is ~1 
W at 7 mK, using reasonable values for the residual gas pressure and achievable surface 
quality for radiation transfer. The resulting estimate for the radiation contribution was 
negligible. CUORE will utilize a similar system to that of the Nautilus Collaboration 
that cools a 2-ton gravitational antenna. That system experienced a parasitic power of 
10/iW from unknown sources. 

The CUORE detector will be cooled by a He^/He^ refrigerator with a cooling 
power of 3 mW at 120 mK. Refrigerators with the required characteristics are techni- 
cally feasible. One example is the DRS-3000 DR model constructed by the Kamerling 
Onnes Laboratory in Leiden. The unit is shown in Fig.|4ll inserted in the dewar. One 
important design feature is the 50 mm diameter clear access to the mixing chamber to 
allow a rod, suspended from an external structure, to suspend the detector array to min- 
imize vibrations from direct connection to the mixing chamber. The temperature of the 
rod will be controlled in stages along its length by flexible thermal contacts to avoid 
vibrations. 

The dewar housing the DR will have a jacket of liquid nitrogen (LN) to avoid the 
need of superinsulation which is not free of radioactivity. The system is designed with 
several tubes passing through the LN bath, the liquid He bath, and the IVC, to allow 
refilling with cryogenic liquids, and for sufficient feed through for electrical connec- 
tions. Liquefiers will provide constant liquid levels for the two baths over long periods 
of operation. The He that evaporates from the main bath and from the IK-pot will be 
recovered in a buffer tank and re-liquefied with a helium liquefier with a capacity of 
several tens of 1/day. The gas will be cooled to ~70 K, then to 10 K, by a two stage 
Gifford-MacMahon cycle. This will be followed by a Joule-Thompson cycle that lique- 
fies it and injects it back into the main bath. The nitrogen liquefier, on the other hand, 
will require only a single stage Gifford-MacMahon cycle to cool a cold finger above the 
LN bath. It will condense the evaporated N2 and let it fall back into the bath. 

This complex cryogenic system will require a constant monitoring system. It must 
be capable of surveying the operating conditions to allow prompt intervention to correct 
malfunctions rapidly. In addition, the operating conditions of the entire cryogenic sys- 
tem must be recorded for future correlation with instabilities in the experimental data 
(section [TT31). 

Besides housing the CUORE detector, the dewar will be able to support an effective 
system of roman lead shields maintained at the low temperature close to the detectors 
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Figure 41: CUORE cryostat and shielding. 
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Figure 42: Temperature instability as observed in a typical MiDBD measurement. 
(FigEland section lO). 

8.1 Development of a Pulse-Tiibe-assisted cryostat 

The experience collected by the proponents in running CUORICINO, and by the Milano 
group in running the previous 20 detector arrays, shows that a constant (or very slowing 
varying, of the order of 1 %/day) level of liquid helium in the main bath of the cryostat 
helps in achieving a steady response of the detectors. Any variation of the main bath 
level translates into small changes of the mixture flow and therefore of the refrigerator 
cooling power. These changes are harmful for the stability of the bolometer output and 
could spoil the effectiveness of the correction procedure (usually referred as "stabiliza- 
tion"). Detectors need to be operated in a DC-coupling mode because their baseline 
level reflects the temperature changes and is used for the detector response stabiliza- 
tion ("Fig.l42land section fT3t. A fundamental point is that the stabilization procedure 
is effective only if the detector operation point (and therefore the baseline level) is kept 
inside a narrow environment of the initial value. This cannot be achieved in presence of 
relatively large or fast changes of the liquid helium level. 

The basic design described above for CUORE, is based on the successful operation 
of CUORICINO and MiDBD in which the helium level was kept almost steady by 
means of a helium re-liquefier. This system is however based on a twenty-year-old 
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technology and has exhibited many weak points during the multi-annual operation of 
the bolometric arrays, such as unpredictable variations of the liquefaction efficiency 
and microphonic noise introduced by the cryogenic heads of the GM cycle. We believe 
that these problems could be even more serious in the CUORE cryogenic system, due to 
the much larger boil-off of the CUORE cryostat (see section [HJl. For this reason we will 
also consider the opportunity to rely on a traditional liquefier which liquefies externally 
the helium gas coming from the main bath and the 1 K pot through recovery lines. In 
this case, an active stabilization of the holder temperature should compensate for the 
base temperature variations induced by the relatively fast decrease of the liquid level 
and the traumatic changes induced by the frequent main bath refills. 

However, we will study (at the R&D level) an alternative cryogenic configuration 
which could allow us to keep the 4 K stage temperature constant for the whole dura- 
tion of the experiment. This solution is based on a new commercial device provided by 
the modem cryogenic technology: a Pulse Tube (PT) cooler. PT is a novel type of me- 
chanical cooling system exhibiting extremely low vibration performance not reached by 
conventional Gifford-McMahon or Stirling coolers. PT coolers consist of a compressor, 
a rotary valve and a cold head. Their advantages are: 

• low maintenance, due to the absence of moving displacers and cold seals; 

• low noise, due to the lack of moving parts in the cold head, which assures a very 
quiet system; 

• low cost, due to its simplicity, which allows the PT to be manufactured at less cost 
than common cryo-refrigeration systems. 

Commercial two stage PT's provide two cold points, typically at 80 K and 4 K, 
with a cooling power of the order of 0.5 W at 4 K. In this framework, a liquid-helium- 
free medium-power dilution refrigerator was installed in September 2003 at the Insubria 
University. This refrigerator does not use cryogenic liquids at all: a PT cools the 80 K 
and the 4 K shields, while the '^He-'^He mixture is condensed by a Joule-Thompson 
expansion after thermalization at the 4 K stage. This refrigerator will be a precious 
test-bench for this new technique, in terms of reliability and noise level. In the event 
of convincing performance of this system, we will design a refrigerator for CUORE 
endowed with an array of PT's with the aim to stop or reduce the main bath boil-off or, 
as an extreme solution, to get rid completely of cryogenic fluids. This solution would 
solve much of the stability problem at the root. 

8.2 Shielding requirements 

Double-beta decay experiments, as well as dark matter searches and searches for rare 
events in general, require deep underground locations and ultra low radioactive envi- 
ronments. This latter requirement is usually accomplished by placing severe constraints 
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on the selection of the materials used and by the realization of proper shielding sur- 
rounding the detectors. In the case that complicated structures like dilution refrigerators 
are used, the problem of their shielding and construction using only low radioactive 
contamination materials is complex. Moreover, since in MiDBD and CUORICINO a 
non negligible contribution to the background rate in the P(3{0u) region was recognized 
as originating from environmental radioactivity, and radioactive contaminations of the 
cryostat structure (sections |6l and \T4h . the design of effective shields (specially those 
directly surrounding the detector) is crucial to reach the sensitivity goal of CUORE. 

Part of the bulk lead shielding will be placed inside of the cryostat, and part outside 
(Fig. |4TT) . This will shield the detector from environmental radioactivity and from ra- 
dioactive contaminations of the dewar structure, also reducing the total amount of lead 
required. A An layer of ultra-low background lead will constitute 3 cm thick walls of 
the cubic structure of the array. This layer will be Roman lead whose ^^°Pb activity was 
measured to be less than 4 mBq/kg [ 174|. As for CUORICINO, the dilution refrigera- 
tor will be constructed from materials specially selected for low levels of radioactivity 
(section [T?b. Nevertheless, these levels might be higher than can be tolerated by the 
sensitivity requirements. The main goal of the inner lead shield will be therefore of 
cancelling their dangerous contributions to the CUORE background level. 

The top of the detector array will be protected by two Pb layers of ~1 m diameter 
and 10 cm thickness each, with a 10 cm diameter central bore to accommodate the 
copper cold finger that supports the detector and the narrow neck of two radiation shields 
of the refrigerator that are at temperatures of 50 mK and 600 mK. The layer close to 
the detector will be made of high quality lead with an activity of 16 Bq/kg of ^^°Pb. 
The upper layer will be made of modem lead with an activity of ^^°Pb of 150 Bq/kg. 
Another layer of lead, 17 cm thick, and with a dimaeter of ~40 cm will be placed 
directly on the top face of the detector. It will be constructed from low activity lead 
of 16 Bq/kg of ^^°Pb. This configuration is designed so that the minimum path to the 
detector from the IVC and the dilution unit is 20 cm of lead. Finally, outside the dewar, 
there will be two 10 cm thicknesses of lead, 16 Bq/kg of ^^°Pb for the inner layer, and 
200 Bq/kg for the outer layer. The lead shield will be surrounded with a 10 cm thick 
box of borated polyethylene that will also function as an hermetically sealed enclosure 
to exclude radon. It will be flushed constantly with dry nitrogen. 

The entire dewar, detector, and shield will be enclosed in a Faraday cage to exclude 
electromagnetic disturbances that also constitute a source of background, albeit at low 
energies, important to dark matter and solar axion searches. The addition of a muon 
veto surrounding the entire structure will be also considered. 
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9 Required laboratory 

The entire CUORE setup will be installed underground inside a proper building which 
will have to guarantee room for all the setup parts (e.g. cryogenics, electronics, shield- 
ings), for a controlled area to be used during assembling procedures and for all normal 
monitoring activities. 

The hut in which the CUORE experiment will be housed will have therefore to 
satisfy the following main requirements: 

• the cryostat and the detector have to be sustained by means of two independent 
structures, in order to avoid any possible transmission of mechanical vibrations to 
the bolometers from the cryogenic facility. 

• the dilution refrigerator has to be kept at a suitable height from the floor to give the 
possibility of an easy dismounting of the thermal and radioactivity shields when 
necessary. 

• the area where the cryostat and the front-end electronics are located has to be 
shielded from electromagnetic interferences by means of a Faraday cage. 

• the area where the detector will be assembled must be a clean room in order to 
avoid radioactive contamination of the array in the mounting phase. 

Following these requirements we made a schematic design of the hut which will be 
divided in three levels as shown in Fig.|43l At the ground floor a platform of reinforced 
concrete, with a central aperture of about 3 m diameter, will be placed. In this aper- 
ture a lifting platform, able to raise at least 50 tons, will be realized with four couples 
screw-nut, driven by four motors simultaneously controlled. The radioactivity shield- 
ing will be placed on the lifting platform; in this way it will be easy to pull down the 
shieldings whenever is necessary to open the cryostat. On the ground level, the pumps 
of the cryostat and the compressors for the liquefiers (or for the pulse tubes) will also be 
placed. 

The dilution refrigerator will be suspended at about 3.5 m from the ground in order 
to guarantee an easy dismounting of the thermal shields (Fig.l43t. The concrete platform 
at the ground floor will be the base for two sets of three pillars each: the first set will 
go up to the floor of the second level and will support a plate to which the cryostat 
top flange will be anchored; the second set will go up to the ceiling of the second 
level and will bear the support structure of the detector. An appropriate design for the 
mechanical suspension for the CUORE array is under study. We have however seen that 
it is relatively easy to realize a suitable device since the only unavoidable requirement 
is the achievement of a low natural vibration frequency of the system (below 0.5 Hz). 



REQUIRED LABORATORY 




Figure 43: Scheme of the CUORE setup installation. 
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The first level will consist of two separated sections: a clean room (at least class 
100), where all the CUORE detector parts will be assembled and stored, and a control 
room for the cryostat. 

The second level will also be separated in two sections. The first one, containing 
the top of the cryostat, the support structure and the front-end electronics, will be sur- 
rounded by a Faraday cage. The second will be used for the data acquisition system. 

9.1 Cooling water requirements 

Depending on the cryostat design, the thermal power required to cool the pumps and 
compressors can vary between 40 and 60 kW. The water cooling circuit of the LNGS 
underground laboratories could in principle easily accomplish to our needs but our ex- 
perience in past years tells us that this circuit has to reach better reliability and purity of 
the water. Therefore we will install a water chiller of adequate cooling power that in the 
future can work as a backup solution to cover possible failures of the laboratory circuit. 

9.2 Power requirements 

In order to guarantee a continuous operation of the experiment, all the equipment has 
to be powered by an Uninterrupted Power Supply (UPS) that should have a power in 

the range 70-100 kVA, depending on the setup design. Only the water chiller can be 
powered by a normal power supply, provided that water circulation pumps are under 
UPS and the water circuit has a buffer big enough to resist to power interruptions for 
the same amount of time the UPS does. 

9.3 Safety 

The CUORE experiment doesn't use dangerous materials and the only safety require- 
ments concern therefore the use of cryogenic liquids. Furthermore the amount of liquid 
nitrogen and liquid helium contained in the cryostat will be of the order of 100 liters and 

then, also in case of an abrupt evaporation the amount of gas released to air is limited 
to few hundreds of cubic meters. Of course, in the case of a cryostat using pulse tubes 
(instead of cryogenic liquids) this issue would also disappear. 

9.4 Space requirement 

The space required to allocate CUORE is lOx 10 m^ and we consider it sufficient for the 
installation of th experiment. Moreover a few square meters outside the hut are required 
in order to accomodate the water chiller (it should be placed in a bypass outside the 
main hall not too far from the hut) and the ventilation system of the clean room (close 
to the hut). 
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10 Electronics 

Although CUORICENO has many fewer channels with respect CUORE, the design of 
its electronics has been developed with the aim of fully satisfying the requirements 
(electrical specification, space occupation, cost, etc.) of CUORE. Therefore, the de- 
scription of the front-end set-up we are proposing for CUORE is very close to the actual 
front-end of CUORICINO. Very few modifications are needed to upgrade the set-up of 
CUORICINO to that of CUORE. 

10.1 General Description 

With the CUORICINO experiment two front-end readouts are being investigated: a con- 
figuration consisting of 38 channels completely operated at room temperature ("warm 
electronics": WFE) and a system of 24 channels characterized by an additional first 
differential unity gain buffer stage at low temperature ("cold electronics": CFE). 

The connection between the thermistors and the first stage of the electronics is via 
a twisted pair of 1 m (CFE) and 5 m (WFE) in length. The main reason for these two 
prototypes is to assess the effect of the microphonic noise of the connecting wires on 
the detector energy resolution at low energy, in order to achieve the maximum possible 
sensitivity at threshold. Since other factors dominate the detector resolution at higher 
energies, an equivalent performance is expected there for the two set-ups. 

Regardless of the configuration considered, the design of the electronic system must 
satisfy many specifications. Among them, the four major issues can be summarized as 
follows: 

1. It should be characterized by a small noise contribution, in the frequency band of 
the detector signal (~one to few tens Hz), generated by the biasing circuitry and 
the preamplifier. 

2. It should accommodate the broad spread of manufacturing characteristics that is 
typical of the bolometer/thermistor modules. In CUORICINO, for example, the 
optimum bias current of some modules is as low as 150 pA, while is as high as 
300 pA for others. As a consequence the signal amplitude exhibits a correspond- 
ing spread. 

3. It should enable in situ measurement of the bolometer DC characteristics. 

4. The design should incorporate a high level of remote programmability in order 
to avoid mechanical interventions in the vicinity of the bolometers during the 
measurements. 




Figure 44: Schematic diagram of the front-end electronic circuit operated at room tem- 
perature (WFE). 

In the following subsections we shall give a brief discussion of the circuitries able 
to address these issues and presently used for the CUORICINO readout. An extended 
description of the CUORICINO electronics can be found in reference 1 175 1. 

Among some of the important tasks addressed by the readout electronics which 
are not reported here, a particular attention was given to the setup power supply and 
the shielding against the electromagnetic interferences. A very precise, multi-output, 
pulse generator instrument was finally developed to allow a continuous calibration of 
the energy conversion gain of the whole array (stabilization). 

10.2 Thermistor biasing and detector wiring 

In Fig. HU and Fig. US are shown the two schematic diagrams of the readout solutions 
for CUORICINO. As can be easily seen, they differ by the presence of the cold buffer 
stage in one of them (Qi and Q2 in Fig. l45t . From left to right of both figures we find 
the detector, its biasing network, the first stage of amplification (warm/cold) and so on. 
Any block characterized by a capital P as the first letter of its label is equipped with at 
least one remotely programmable functionality. Thermistor biasing is accomplished by 
applying a programmable voltage across a series combination of the thermistor and two 
load resistors. 

For the set-up of Fig.|44lthe load resistors, PRl, are also at room temperature. The 
analysis and minimization of the contribution to the low frequency (LF) noise from large 
resistors under bias conditions, is the subject of an independent report L176|. Under 
computer control, the biasing resistors can assume values in the range 11-54 Gfi. This 
facilitates the measurement of the bolometer's I-V characteristics at large bias. The 
impedance at the input of the electronics is sufficiently high that parasitic conductance 
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in the circuit board is an issue that demands great attention. Precautions that we took to 
suppress parasitic conductance were: 

• segregation of the two high impedance detector nodes by a circuit board layer; 

• strategic placement of guard traces around the solder junctions of high-impedance 
component leads; 

• maintenance of a minimum spacing of ~5 mm between high-impedance traces; 

• placement of slots at critical locations on the circuit board to avoid conductance 
through minute amounts of soldering residue that sometimes escapes the standard 
cleaning procedures. 

Altogether, these solutions yielded excellent results. 

In the set-up of Fig. |45lthe pair of load resistors is also put close to the buffer 
location, at low temperature. The rest of the biasing set-up is unchanged, but the 
value is 54 G^l and cannot be switched to 1 1 G^l. 

The application of the bias voltage to the load resistors is accomplished by the pro- 
grammable bias attenuator (PB-i-, PB-), an especially designed Digital to Analog Con- 
verter (DAC) characterized by an optically isolated digital input. It can generate the 
voltages PB-i- and PB- of both Fig.|44|and Fig.|45l symmetrically with respect to ground 
(PB-=-PB-i-). The programmable bias attenuator can also provide the additional fea- 
ture of inverting the polarity of the bias voltage. This feature is very useful, since by 
reversing the bias voltage during DC characterization, it is possible to account for the 
thermoelectric effect operating in the leads that connect the front-end electronics to the 
thermistors fTTl]. To simplify measurements and improve accuracy, an additional bus 
PAUX on Fig.|44|and Fig.|45lcan be connected to any selected channel to perform the 
voltage reading with a very precise multimeter, independently of the acquisition system. 
The biasing network is equipped with a heavy low pass filter that reduces to a negligible 
level whatever noise originates from the power system or from the network itself. 

In CUORICINO, the wiring of all the electrical elements at low temperature has 
been realized by means of shielded twisted cables (superconducting coaxial cables 
with a CuNi and 100 /im diameter NbTi twisted pair). Five stages of thermalization 
cut the thermal conductance of the cables themselves. The five stages are located at 
4.2 K, 1.2 K, 120 mK, 50 mK and 5 mK plates and are composed of sapphire windows 
(10x20x1 mm'^) that are sputtered with gold on the bottom, and have 12 1 mm wide 
tracks (gold-platinum) on the top, where the twisted cables are soldered. 

We are studying the possibility of upgrading this set-up with a less expensive so- 
lution. First, the sapphire windows are now going to be replaced with flexible Capton 
windows with thinner and shorter Cu tracks. A simple pressure on their center against 
the cold plate allows an effective thermalization of these windows. The next step will be 
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Figure 45: Schematic diagram of the front-end having the first stage, the differential 
buffer Qi,Q2,at cold. 

the realization of very thin Cu or CuNi tracks sputtered on a flexible Capton substrate. 
This new approach seems to be feasible thanks to the small radioactive background 
measured for Capton. 

10.3 Preamplifier and Cold buffer stage 

The design of the preamplifier (DVA in Fig.HUand Fig.H3l satisfies four main require- 
ments: 

1 . Because the typical source impedance and signal characteristics of a bolometer, 
the voltage amplifier is the conventional approach Ill78i . Its performance is not 
superior to a charge or current configuration II179I . but it allows floating inputs 
which make biasing naturally compatible with DC coupling. 

2. The long leads (in the WFE option) connecting the preamplifiers with the bolome- 
ters have the drawback of acting as antennas for microphonic noise. We ad- 
dressed this problem by making the preamplifier input differential. Compared 
with a single-sided input the amplifier noise power is doubled, but the trouble- 
some (common mode) microphonic noise originating in the leads can be effec- 
tively suppressed. Moreover, the differential configuration suppresses interchan- 
nel crosstalk. 

3. The signal bandwidth of our bolometers extends from 1 Hz to few tens Hz. We 
then require that the preamplifier noise, both series and parallel, be minimized 
in this frequency interval. The CUORICINO preamplifier is characterized by a 
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voltage noise floor of 3 nV/^/Hz, that becomes 4 nV/V^ 2; at 1 Hz. Its noise 
current is 4 ikl\f}iz. 

4. Monitoring the baseline of each bolometer, and correlating the baseline shift with 
the shifts in the overall system gain, requires a high thermal stability of the pream- 
plifier. A special circuit has been developed that allows the drift to be maintained- 
below 0.2 II V/C. 

With the aim of minimizing the microphonic noise due to the length of the leads, 
we have realized the set-up of Fig. |45l It's main feature is that the first buffer stage of 
amplification operated at low temperature. The box housing the preamplifiers is in fact 
anchored to the 4.2 K plate of the dilution refrigerator, next to the detectors. The buffer 
stage is then warmed in order to reach the optimum temperature for Silicon, around 
110 K. Each cold stage is a differential buffer (Qi and Q2) composed of two Silicon 
JFETs that work in the source follower configuration flSOT. In this way the length of 
the leads that connect the detector to this additional buffer stage is reduced to about 1 m, 
while the remaining wires towards the external stage at room temperature are driven by 
the low output impedance of the differential buffer. The cold buffer stage features a 
negligible current noise, a floor noise of about 2 vlNI\/1Iz and a noise of about 10 
nV/V^ 2; at 1 Hz. Each cold buffer, together with the load resistors, is located on boards 
that accommodate 6 channels. The boards are fabricated on a substrate made of PTFE 
25% reinforced fiberglass. This solution achieves two goals. First, the small thermal 
conductance of the PTFE minimizes the injected thermal power needed to increase the 
JFET temperature from 4.2 K to the optimum 1 10 K. Second, the radioactivity level of 
PTFE is lower by a factor of about two with respect to pure fiberglass. In addition, the 
printed circuit boards have been realized with a thickness of only 0.5 mm thus halving 
their mass. The boards are suspended by means of (low thermal conductivity) nylon 
wires inside thermally shielded copper boxes which can contain two of them. Although 
a bit more complex than WEE, the CFE option has the advantage of offering an effective 
noise reduction (see Fig.l38b. 

For both configurations, the input offset is the sum of the contributions from the 
asymmetry of the JFETs and the bias voltage of the NTD thermistor. The JFET asym- 
metry may contribute up to ± 20 mV to the offset. During normal operation, thermistor 
biasing contributes no more than 20 mV, but for load curves acquisition the bias volt- 
age must be adjustable to both polarities and to a magnitude as high as |50my|. The 
offset corrector circuit (POA in Fig.|44land Fig. l45t is designed to compensate input 
offsets from -80 mV to -1-8O mV, under the control of the digital section. The POA was 
implemented with a circuit solution that minimizes any noise contribution, especially at 
LF II175I . Once fired by the computer, the POA adjusts the output offset of the selected 
readout channel to a default value of 1 V. By request, this value can be changed remotely 
to an arbitrary level. 
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An AC (DC) test signal can be injected directly to the DVA inputs in order to test the 
front-end of any selected channel and the load resistors of the WFE channels. During 
this operation, a programmable system of relays disconnects the detector (or the cold 
stage) from the DVA. 

10.4 Programmable gain amplifier 

The outputs of the preamplifier drive the differential inputs of a programmable gain 
amplifier, (PGA in Fig. HU and Fig.H31). The amplifier output to the preamplifier input 
ratio (voltage gain) can take values in the range 220-5000 V/V. The differential outputs 
travel on an analog bus to the filter circuit PGF, which is followed by analog to digital 
conversion of the acquisition system. 

10.5 Antialiasing filter and analog trigger 

In CUORICINO, optimal filtering algorithms are applied to the acquired detector sig- 
nals (sec. fT3b . To ensure adequate frequency response we have implemented an an- 
tialiasing filter (PGF in Fig. HU and Fig.H3t at the downstream end of the analog signal 
processing. This is an active six-pole Bessel low-pass filter, which yields a roll-off of 
120 dB / decade. In order to adapte the useful bandwidth to each bolometer module, the 
cutoff frequency can be adjusted from the computer. 

The input of the trigger circuit (PTRS in Fig. |44l and Fig. l45l) is AC coupled to 
the output of the antialiasing filter. The core of PTRS is an analog amplifier. The 
gain of PTRS is adjustable from the computer. The AC coupling produces a negative 
excursion of the trigger circuit input that, if uncompensated, would result in a shift of 
the baseline. A baseline restoring section II181II has been added to PTRS to suppresses 
this effect. Pulses from a bolometer may exhibit a trailing undershoot because of the 
coupling of the inductive part of the bolometer's dynamic impedance to the parasitic 
capacitance in shunt with it. After differentiation by the AC coupling, the undershoot of 
large pulses has the potential to re-trigger. An RC network 1 182| within PTRS cancels 
the effect of the undershoot and thus suppresses re-triggering. Its time constant has been 
made adjustable from the computer. The final element in the trigger circuit is a two- 
pole Bessel low-pass filter. The cutoff of this filter tracks the cutoff of the antialiasing 
filter PGF, but is higher in frequency by 4 Hz. A single trigger-Bessel module card 
accommodates three channels. 

10.6 Calibrating Pulse Generator System 

Very long measurements need an accurate energy calibration. A radioactive source can- 
not be left or put (on request) close to our detectors, because of its unavoidable contri- 
butions to the dead time and background level. Our solution is based on the use of both 
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a radioactive source positioned close to the array (typically once a month, just outside 
the cryostat) and an always-present very stable electrical calibrating system (heater), 
which inject signals at a faster rate. The radioactive source fixes the energy scale of the 
detector while the heater effectively corrects gain variations (sec.fT3t. 

The electrical calibrating system is composed of a small heater resistor glued on each 
crystal. Across any heater resistor a small voltage pulse is developed which dissipates a 
small power, that generates a signal very similar to that of an impinging particle. 

The accuracy of our calibrating system depends on the stability of the heater resistor, 
whose characteristics are described in 1 151 1, and to the electrical pulse generator. To this 
aim we have developed and are using a very stable pulse generator having the following 
features: 

1 . It is modular and provides 4 independent outputs for each card used. 

2. Each output generates square pulses that can be completely programmable in am- 
plitude and width, and can span a range corresponding to pulse amplitudes equiv- 
alent to a few tens of keV up to more than 20 MeV. 

3. Each pulse occurrence can be driven by an external input (e.g. DAQ), in order to 
allow for a complete tagging. 

4. The stability of the pulses is within 2 ppm in any parameter. 
A complete detailed description of this device is given in II183I . 

10.7 Voltage Supply system 

The voltage supply system is composed of two main stages. The first stage fulfils a pre- 
adjustment of ± 13 V by means of several 75 W switching regulators. Special filtering 
sections are used to suppress any high frequency spike. The second stage outputs pro- 
vide the desired voltages in the ranges -12—1-12 V for the output stage of the front-end 
and in the range -10— I-IO V for the cold stages and the preamplifiers. The ± 10 V out- 
put is used also as the reference for the offset correcting circuits. For this reason a very 
stable (a few ppm/°C drift) and low noise (1.3 yU Vpp in the 0.1 Hz to 10 Hz frequency 
range, at a rated current of 1.5 A circuit) was designed and built [ 184|. A single second 
stage voltage supply supports 24 complete channels. 

10.8 Digital control 

A digital control system (DCS) enables a computer assisted regulation of the most sig- 
nificant analog front-end parameters: the thermistor bias current, the input offset correc- 
tion, the gain, the bandwidth, the application of test signals to the inputs, the calibrating 
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pulse injection to the detectors, etc. During data acquisition the DCS is generally idle. 
In this state even the DCS clock shuts down, thus eliminating the possibility that DCS 
activity will contribute to noise on the analog outputs. 

To minimize elecromagnetic interferences (EMI) and ground loops, the communi- 
cation between the DAQ (sec. fTTT) and the DCS is made with an asynchronous, serial, 
bi-directional optical fiber link bus. 

10.9 EMI Shielding 

We foresee, as is partially done now in CUORICINO, surrounding the entire cryostat 
and the front-end electronics with a Faraday cage able to shield also low frequency mag- 
netic fields. For this purpose the internal walls of the cage will be covered with a sheet 
of 0.1-0.2 mm thickness of SKUDOTECH, a very high permeabiUty metal |185|. A 
sound absorbing material will also cover the internal walls. The EMI shield will be split 
in two sections. A first Faraday cage will surround the cryostat and its lead/polyethylene 
shields. The very front-end, fed by the detector connections (output signals and input 
bias), will be organized in 12 towers (each composed of 6 standard 19 inch-racks) and 
connected to the second stage located in a smaller Faraday cage sharing one of its walls 
with the main Faraday cage. The small Faraday cage will house the antialiasing filter, 
the trigger circuit and the acquisition system optically connected via TCP/IP to the out- 
side world. The choice for a set-up split in two separate rooms allows a distribution of 
the power dissipation over a large space and an effective reduction of the digital circuits 
EMI. 

10.10 Noise considerations 

The above described analog WEE produces a noise of the order of 50 nV (FWHM) 
referred to its input, over the 10 Hz bandwidth of the signal. The 54 G^l load resistors 
produce 80 nV in the same bandwidth so that the combined noise is ~ 100 nV (FWHM). 
On the other hand, for the CFE option, due to the presence of the cold buffer stage, the 
series power noise is increased by a factor of about two, while parallel noise is decreased 
by a factor of at least three, mainly due to the use of the cold load resistors. In this way 
the overall noise contribution becomes about 90 nV (FWHM). In both cases the noise is 
considerably below the 1.0 fiY (FWHM) level observed for one of our best performing 
bolometers. We believe therefore that the analog front end is more than adequate for the 
CUORE needs. 

The two front-end solutions have the same level of noise, but the presence of the 
cold buffer stage, which reduces the high impedance path, should further improve the 
common mode noise due to the vibration of the leads. Unfortunately, there is another 
source of microphonics which is not a common mode: the mechanical friction between 
the detector and its mounting. If this source would result in a noise level much greater 
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than that of the connecting leads, the use of a cold buffer stage would be unnecessary. 
This would allow a reduction of the connecting leads by a factor of two and would 
avoid thermal power injection inside the dilution refrigerator. It should be noticed that 
the reduction of the parasitic capacitance obtained by shortening the high impedance 
path is not a concern in the case of our very slow signals. 
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11 Data acquisition system (DAQ) 

In this section we describe the architecture of the readout and DAQ system for CUORE. 
In particular we report the main requirements that are driving the system design and 
describe the general architecture of the system oulining a few open options. We also 
give the block diagram of the main readout board and finally describe the main DAQ 
system and its software layout. 

11.1 Requirements 

As described in the previous sections, the signals coming from each crystal are ampli- 
fied and shaped by a preamplifier that generates pulses well described by the following 
typical parameters: 

• Voltage signals with maximum peak amplitudes as high as 10 V; the signals can 
be either differential or unipolar; due to the large energy range, the peak amplitude 
can also be as low as few mV. 

• Rise times of the order of 30 ms, falling times up to an order of magnitude slower. 

• Total duration up to several seconds. 

The DAQ system must be designed to sample and read out the signals of about 
1000 channels^ of this type with virtually zero dead time, independent readout among 
channels, self-triggering capability, and, last but not least, optimized cost. 

In order not to spoil the energy resolution over the whole range (from 5 keV up to 
15 MeV), real 16 bit sampling is required. In addition, the pulse must be sampled for 5 
s, in order to get a good measurement of the tail and of the background level. A typical 
sampling rate of 1 kHz seems adequate, but it is highly recommended to allow for a 
programmable sampling speed around this value. With the aforementioned values, the 
typical event size for one channel is 5000x2 = 10 kbytes. 

It is also very important to measure the background level (baseline), both before 
and after the end of the pulse, therefore it is required to have a suitable internal circular 
memory and a logic capable of recovering the pre-trigger signal conditions. 

The choice for real 16 bit sampling poses severe requirements on the noise level 
acceptable at the channel input. Assuming a maximum range of 10 V, the noise level 
must be below 150 //V r.m.s. 

'As an option, it has been proposed to add a second sensor to each crystal. This would of course 
double the number of channels, but would have a negligible impact on the architecture proposed here, 
except of course an almost linear increase of the total cost of the system. 
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Being that CUORE is designed to search for very rare events like f3/3 decays or dark 
matter induced nuclear recoils, the requirements for the DAQ system for data throught- 
put are not severe at all. The maximum foreseen data rate is given by the calibration 
procedure, that might produce at most few events per minute per channel. As described 
in section fTL?! this rate does not pose any special requirement on the DAQ system, even 
assuming a complete sampling of the channels without zero suppression. 

11.2 General architecture 

The CUORE DAQ system must be able to digitize, acquire, analyze online, store and 
display the pulses generated by 1000 Te02 bolometers. 

Each channel must be completely independent from the others and the readout must 
introduce zero dead time. 

In order to minimize the cost and allow for a maximum flexibility in the on-line 
analysis and triggering, the system proposed here is mostly based on the large and cheap 
computing resources that are now commercially available, reducing the complexity, and 
therefore the cost, of the readout boards as much as possible. 

The general structure of the system is shown in Fig.H^ It is composed of a set of 
personal computers, a network infrastructure, a PCI based VME interface and 5 or more 
VME crates with custom designed digitizing boards. 

We propose to design and build a simple digitizing VME board with full 16 bit 
dynamic range (possibly 24), large internal circular buffers, optional self-triggering ca- 
pability with a programmable threshold and digitizing speed programmable in the range 
from 0.01 KHz up to 10 KHz. 

Each board can probably handle 12 channels, to allow reading 240 channels per 
VME crate (5 VME crates^ would then be required with 1000 channels). 

In this structure the digital boards are just pure digitizers, that continuously convert 
the input analogue signal into a stream of digital words; when a trigger occurs, either 
generated by the remote front end electronics or internally generated by the discrimina- 
tor, the words stored in the internal buffers are copied to another memory that can be 
accessed by the main DAQ software via VME bus. The structure of the board guarantees 
that the channels remain independent and fulfils the zero dead-time requirement. 

Each board issues an interrupt signal on the VME bus whenever a trigger occurs; in 
this way the main software may work either in interrupt or in polling mode. 

The large and cheap computing power available nowadays yields a large flexibility 
in the possible on-line filtering and processing of the data. Full on-line analysis is 
possible with this structure, allowing a fast and efficient monitoring of the experimental 

^Should the two sensor options be adopted, the number of required crates would be 10; additionally, 
we might succeed in designing a 16 channel digitizing board that would reduce the number of VME crates 
to 4 or 8. The general structure of the system would not be changed at all. 
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Figure 46: Block diagram of the proposed DAQ system. The digitizing boards are 
housed in 5 VME crates, each controlled and read by a dedicated Linux PC that uses a 
VME-PCI interface to access the VME bus. All computers are connected via a network 
(normal 100 Mbit ethemet or fast optical link if required) to a set of consoles that provide 
graphical user interface and event display. A custom designed software runs on all 
computers. 
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Figure 47: Base design of the digitizing board. In this "minimal" version the board sim- 
ply has the function to convert the input signal into a stream of words that are transferred 
to the DAQ computers via VME bus. 

The next two sections give a somewhat more detailed description of the digital board 
and of the computing and software architecture. 

11.3 Readout board 

We are considering two possible block diagrams for the digitizing board. The first is the 
one sketched in the previous section and shown in Fig. |47| 

In this design, the digitizing board is very simple and most of the job is left to the 
software that runs on the DAQ computers; in fact, the board is a pure free-running 
digitizer, with the minimum amount of logic needed to recover the data from a few 
hundred ms before the trigger until several seconds after the trigger. The trigger can 
be either external, i.e. generated by the front end electronics and sent to the board via 
a NIM or TTL signal, or can be internally generated with a simple discriminator. To 
achieve the required resolution (real 16 bits), we will probably put two 12 bit ADCs in 
parallel, with a precise amplifier in front of the least significant ADC. In this way 24 
bits will in fact be available. The amplifier gain will be chosen to leave a large overlap 
between the two ADCs, allowing a precise intercalibration. 
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Figure 48: Optional enhanced design of the digitizing board. In this version the board 
is equipped with a DSP that can perform additional computing functions like event 
selection or digital filtering. 

When a trigger occurs, the FPGA logic will transfer the data (a programmable time 
window around the trigger signal) to a Dual Port Ram that can be accessed via VME bus. 
A VME interrupt will be also generated when this transfer is completed, allowing the 
DAQ computer to start read-out. It will be possible to use this board either in interrupt 
mode or in polling mode. Block transfer D32 will be supported too, allowing very fast 
transfer. The Dual Port Ram will be large enough to keep several minutes of data at 
normal trigger rate during calibration. 

We are also considering a second approach, a little more complex from the point 
of view of hardware design, that includes a DSP in the digitizing board (see Fig. B5b . 
With this solution, the DSP could be used to perform a first level data filtering and 
implement "smart" trigger conditions in order to reduce the data rate. In our opinion, 
with the data rate foreseen for CUORE, it will be possible to perform all this higher level 
functions using the computing power of the DAQ system, allowing for more flexibility 
and reducing the cost of the system. Therefore we prefer the scheme shown in Fig. HT] 
as our baseline design and leave the second design as an option. 
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11.4 Computing system and DAQ layout 

The system shown in Fig.H^ with the baseline option for the digitizing board shown in 
Fig.|47l requires a significant amount of computing power because most of the functions 
are delegated to the software. 

As a baseline design, we propose that each VME crate will house 20 boards, each 
of them capable of digitizing 12 channels. Therefore we need 5 VME crates to read 
1000 channels. Each VME crate will be connected to a PC via a PCI- VME interface. 
These 5 PCs will be connected via network to another set of PCs that will perform event 
building, online analysis, graphical display, system control and monitor. 

The network infrastructure will be rather simple, with a single switch star-connected 
to all PCs. Normal 100 Mbit ethernet is probably more than enough for CUORE but it 
would be easy to implement this structure using 1 Gbit optical link without any archi- 
tectural modifications. We think, therefore, that this very standard network structure is 
adequate. 

We propose to use Linux based PCs for all computers and to use root as a software 
environment. The software running on each computer will be custom designed while 
the consoles and the graphical interfaces will be root based. A possible WEB interface 
to the system is foreseen, though network security problem will probably limit its use. 

11.5 Slow control 

During a measurement, all the parameters of the apparatus must be monitored and regis- 
tered. This is particularly required in long running time experiments where the acquired 
data can be "contaminated" with problems related to experimental system conditions. 
In particular in the case of cryogenic detectors many inconveniences can arise from the 
cryogenic apparatus. The experience gained with MiDBD and CUORICINO detectors 
showed us the fundamental role played by the temperature instabilities and by the vibra- 
tions induced by the malfunctioning of the cryostat and liquefier. In order to correlate 
these frequently unavoidable problems it is very important to continuously check the 
main parameters of the cryogenic system (e.g. pressures, mixture flows, cryogenic liq- 
uid levels and others). The cryogenic apparatus must be equipped with transducers able 
to transpose these data in a voltage value to be acquired by the DAQ during the exper- 
iment. This environmental data acquisition system can be implemented directly in the 
DAQ or be designed independently and correlated off-line with the experimental data 
by means of a time synchronization between the two systems. The choice of the best 
approach can definitely be taken later on; the control system and DAQ designs have in 
fact little intereference and can evolve independently. At the moment we are not plan- 
ning an active automatic control of the apparatus. This could be designed later, when a 
better definition of the cryogenic and detector apparatuses will exist. On the other hand, 
we would like to implement a human remote control(e.g. WEB interface) in order to 
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maintain the apparatus in the best measurement condition without the necessity of the 
constant presence of a shifter inside the underground laboratory. 
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12 Construction procedures 

The only possibility for the CUORE collaboration to satisfy the sensitivity requirements 
will be to use severe procedures for the construction of the experimental setup. Accord- 
ing to the experience gained in CUORICINO, this will imply rigorous material selection 
and construction procedures. 

12.1 Material Selection 

One of the most critical points of the CUORE experiment will be the selection of all the 
materials used in the construction of the detector. Since the selection criteria must result 
in the minimum possible contribution to the background level, it is obvious that this will 
be accomplished only if the starting quality of the materials will be sufficiently good 
from the point of view of radio-purity. We can divide the materials needed for CUORE 
in three main categories: i) materials for crystal growth; ii) materials for the detector 
structure and iii) materials for the shielding structures. Material selection is important 
also for other parts of the experimental setup (e.g. the cryostat) but in principle these 
parts can be heavily shielded so the required limits on their radioactive contaminations 
are less stringent (section lOb . 

12.1.1 Materials for crystal growth 

Te02 is a molecule composed of only two atomic species: Te and O. In practice if we 
want to produce a good stoichiometric material, the oxidation of the Te metal must be 
realized in few well defined steps. In each step the selection of the compound used 
will be very critical. The first material to be chosen is obviously Te. The analyzed 
metallic Te for the CUORICINO production showed limits below the sensitivity of the 
HPGe detectors installed in LNGS (apart from the well known contamination in ^^°Pb 
). Te02 powders were also measured at the LNGS low activity underground laboratory. 
Contamination values are reported in table [TSl The presence of ^^°Po can be ascribed 
either to a genuine ^^°Po contamination or to ^^°Pb . Since Pb contaminants should be 
limited in the CUORE experiment, we will try to investigate this point in detail. Two 
acids are usually used for metallic Te oxidation: Nitric and Hydrocloric. The ratio used 
in this preparation is: 1 kg of Te, 1.3 kg of HCl and 2.5 kg of HNO^ II186l . So we 
need to check with the same accuracy the radio-purity of the acids and of the Te metal. 
For CUORICINO we have measured the two acids involved in the crystal production 
using ICPMS with quite good results: HCl showed a contamination of 1.3 ppt of ^'^^U 
and <1 ppt of 232Th, while HNO3 showed 13.1 ppt of ^^^U and 112.3 ppt of ^^'^Th. 
One of the main problems in CUORE will be the systematic control of these values 
during all the crystal production operations. There is another very critical point at this 
stage: a large quantity of water (few liters per kg) is used to "clean" the Te02 powder 
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just after their production and before the growth procedure. This operation can be of 
course very dangerous because of the possible presence of radioactive contaminations 
in the water that can be transferred to the powder. The impact of these contaminations 
could explain the higher contamination observed in the Te02 powders with respect to 
the starting Te metal. The growing procedure usually reduces the contamination lev- 
els of the final crystal with respect to the bare material. However it is usually difficult 
to quantify the actual reduction in the produced crystals. On the other hand, an in- 
crease of the number of growing cycles - although undoubtely useful - would be very 
expensive. The last critical point in the preparation of the crystals will be the surface 
treatment made in order to assure a good contact between crystals and semiconductor 
chips (thermistors and silicon resistances). The data obtained with MiDBD showed that 
a relevant contribution to the background in the double beta decay region could be iden- 
tified with the presence of surface contaminations on the crystals. In the production of 
the early 3x3x6 cm'^ Te02 crystals, this consideration was strongly supported by the 
measurement of the lapping powders used by the factory. In CUORICINO we carried 
out the final lapping procedures using a specially selected AI2O3 powder. This showed 
a radioactive contamination few orders of magnitude lower with respect to the previous 
ones (tab.fTSt. However, because of the large spread in the grains dimension of this pol- 
ishing powder, only roughly smooth surfaces could be obtained. Alternative approaches 
are in principle possible (e.g. surface acid attack, ion beam sputtering, etc.); the careful 
test of these methodologies and the development of powerful dignostic techniques will 
be considered in the early CUORE R&D phase (see sectionfT6b. 

12.1.2 Materials for the detector structure 

The experience gained in MiDBD and CUORICINO indicates that all the structure of 
the CUORE detector should be mainly made of Cu. Only minor components will be 
built with different materials. In the case of Cu, we have indications that very low 
limits for the bulk contaminations of the U and Th chains can be reached in specially 
selected productions. We have measured a number of samples for CUORICINO and 
found that the Oxygen Free High-Conductivity Cu is definitely better than the stan- 
dard electrolytic Cu. Moreover some specific productions, made by our colleagues in 
Germany [ 187 1, have shown very low limits in radioactive contaminations, at levels suf- 
ficient for CUORE sensitivity requirements (see section fT4ll. We must plan and control 
the production of Cu for CUORE directly with the foundry in order to check the con- 
tamination of the material immediately. Moreover the Cu must be stored underground 
immediately after production to avoid cosmic ray activation. This poses severe condi- 
tions on the physical size of the produced Cu in order to minimize the time needed by 
the mechanical factories for the production of the final pieces of CUORE. 

Another important material in the detector structure is PTFE. The spacers holding 
the detector in the right position are made with this material and therefore the radioac- 
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Table 15: Radioactive contamination measurements for CUORICINO materials. All 
values are in mBq/kg. AI2O3-I is the original lapping powder used at SICCAS while 
AI2O3-2 is a certified 0.7 fim Sumitomo powder. 
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tive contaminations of PTFE can be very dangerous for the background level of the 
detector. For organic materials there are two different problems: first it is necessary to 
have very low levels of intrinsic radioactivity, second it is crucial to control the exposure 
of these materials to air radioactive pollution (typically Radon). In general the radioac- 
tivity of organic compounds like PTFE, Nylon etc can be very low, as demonstrated by 
the Borexino collaboration II IB 811 . Concerning the diffusion of Radon it is necessary to 
maintain the selected materials in a very clean environment and, to minimize the pos- 
sible radioactive contamination present in the volume next to the surface, the external 
part must be mechanically removed before the realization of the spacers. 

The other materials involved in the detector construction are characterized by very 
small masses. Also in this case however, a good radiopurity level is required in order 
to avoid any possible dangerous contribution to the background level. In that sense, 
the selection of the epoxy glue (connection between chips and crystals), of the gold 
bonding wires (electrical connections), of the electrical pins, etc. must be carefully 
controlled. The indications obtained from our previous tests, show very good levels for 
epoxies while are non conclusive for gold wires and pins. All these components will 
be measured again to improve their contamination limits; alternative approaches (e.g. 
neutron activation) will also be considered. 

The last critical part of the detector assembly is represented by the wires used for 
the electrical connection between the detector and first thermalization stage on the inner 
lead shield (they are not shielded). It is very difficult, at the moment, to define their 
actual contribution to the background level since we still lack a defined design of the 
final electrical connection. As a matter of fact we will choose the best among a number 
of different wire types (constantan, superconductor, etc) and will then design the wiring 
system using only previously selected radio-clean materials. 

It should be stressed that all the materials must match two different requirements: 
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they must be characterized by low radioactive background and, at the same time, ensure 
a good performance in a cryogenic environment. Of course, is not obvious that to sat- 
isfy these very different conditions unavoidable difficulties might not arise in correctly 
defining the material for the specific application. 

12.1.3 Shielding materials 

The radioactive shields will be made of lead. To minimize the contribution from the 
lead itself, we need to build different concentric shields made with lead characterized 
by different contents of ^^°Pb . It should be stressed that the presence of this nuclide can 
contribute to the background spectra with the bremssthrahlug photons produced by the 
beta decays of the ^^°Bi. The internal lead shield, to be placed as close as possible to 
the detector at cryogenic temperatures, will consist of ancient roman lead that contains a 
level of ^^°Pb lower than 4 mBq/kg [ 189|( section lOt . The intermediate lead shield can 
be realized with commercial lead with a ^^°Pb contamination of the order of few Bq/kg. 
The external part will be constructed with standard commercial lead that normally has 
contaminations of the order of few tens of Bq/kg in ^^°Pb . All the above requirements 
were already satisfied by the lead used in CUORICINO. 

12.1.4 Other components of the experimental setup 

There are some parts of the cryogenic environment that need careful design, not only 
from the point of view of the cryostat performance, but also for their possible contribu- 
tions to the radioactive background. In particular, although the upper part of the cryostat 
will be shielded with respect to the detector volume with a heavy layer of lead, its most 
massive components should be constructed with materials characterized by an accept- 
able level of impurities. We can indicate a rough limit for these materials and reduce 
at a minimum their required amount. Of course, the requirements connected to pro- 
ducing a well performing cryogenic apparatus must also be taken into account. Some 
possible materials that satisfy both requirements have already been found during the 
last few years for CUORICINO. This job will continue in the future with the systematic 
selection task needed for the specific application to CUORE. 

12.2 Construction environment 

The Te02 crystals will be fabricated at the Institute for Ceramics in Shanghai (SICCAS), 
China, from low-radioactivity ore. Our prior experience with similar crystals fabricated 
for MiDBD and CUORICINO indicates that the final stages of lapping and polishing 
must be done in a clean, radioactivity-free environment, recognizing that these detec- 
tors are the most-sensitive detectors available for traces of radioactive materials. Two 
possible options can be envisaged for the choice and set-up of this environment: a) a 
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conservative option, which follows strictly the CUORICINO experience; b) an innova- 
tive option, which foresees a substantial use of out-sourcing, therefore requiring the de- 
velopment of new environments or the modification of existing environments presently 
not suitable for the delicate operations involved in the CUORE assembly. The main 
points of the two options are outlined below. 

(a) In the conservative option, the crystals will be grown in China, cut in a cubical 
shape and lapped using alumina powder down to a grain size of about 15 fim. 
The control of the residual activity in the lapping pads and in the alumina will be 
performed (see subsection 112. II) . The lapped crystals will be delivered by ship 
in order to limit exposure to the high-energy hadronic component of the cosmic 
rays and stored underground. The final most delicate polishing operation will be 
accomplished directly in the Gran Sasso laboratory, possibly underground or at 
least keeping each crystal outside the underground laboratory just for the time 
required for polishing. A special environment has to be developed for this pur- 
pose, consisting of a clean room where several polishers can be installed in order 
to proceed in parallel with many crystals. A radon-free atmosphere consisting 
of synthetic air should be realized in the clean room. Presently, in the hall C of 
the underground Gran Sasso laboratory, there exists a facility which with some 
modifications, could have the required features: it is the combination of clean 
rooms developed for the BOREXINO experiment. If available at the time of the 
CUORE construction, this would be an ideal environment for the polishing pro- 
cedure. Otherwise, a similar structure must be built. After polishing, each crystal 
will be stored underground in a radon-free container. Once the polishing opera- 
tion is completed, in the same clean room where it was realized, the thermistors 
will be attached to the Te02 crystals. Finally, the crystals will be inserted in their 
copper frames and the thermistor wires connected to the read-out pads. Each 4- 
detector module will be stored independently in a radon-free container and moved 
to the cryostat area. The final assembly of the modules will be realized directly in 
the lower level of the cryostat room, that will have clean room features. As far as 
the whole detector will be completed, it will be connected to the dilution unit of 
the refrigerator and will be ready for the final wire connections. 

(b) In the innovative option, a large fraction of the detector assembly work will be per- 
formed in China, from where the crystals will be shipped in the form of completed 
4-detector modules in special radon-free containers. A clean room environment 
(class 1000 or better, with radon-abatement) must be set up either at SICCAS 
or at some other facility under the supervision of LBNL scientists. These same 
scientists will develop procedures to be used for the production of the crystals. 
The final polishing operation will be performed by skilled Chinese personnel, 
but oversight, either in the form of an employee on site or remote access, must 
be maintained. Finally, these crystals should be instrumented and packaged into 
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elementary modules, i.e. groups of 4 crystals, frames, thermistors, and heaters be- 
fore leaving the clean area. Even this operation will be performed by supervised 
Chinese manpower, after an appropriate training period. Once all the elemen- 
tary modules have been completed, they will be shipped to the underground Gran 
Sasso laboratory as quickly as possible to reduce cosmic ray activation. The final 
assembling will be performed in Gran Sasso as in option (a). 

An intermediate solution between the two options is possible. The Chinese part 
of the work in option (b) could be terminated at the final polishing phase, while the 
delicate operation of thermistor attaching and 4-module realization would be realized 
in the Gran Sasso laboratory. In this way, no polishing set-up is required in Europe, 
but the detector assembly operations will be entirely accomplished underground in the 
Gran Sasso laboratory, where skillful physicists and technicians can more easily check 
the mounting procedure. 

An important point to be considered for the choice between the two options is the 
total exposure to cosmic rays of the crystals from their growth to their final storage un- 
derground. Of course, this exposure will be much shorter in case of option (a). An esti- 
mate of the production rate of radio-nuclides generated by nuclear reactions induced by 
high-energy hadrons in the crystals has been already done using the COSMO program 
(see section [T?l). Cosmogenic ®°Co in Te02 crystals could be a limiting contribution 
for CUORE so that exposure to cosmic rays must be limited as possible. Monte Carlo 
simulations have shown that a total exposure period of 4 months to cosmic rays would 
give a contribution to the PP^Ou) CUORE background below the sensitivity requirement 
(see sectionfT4b. Moreover, a recent experimental determination of the ^°Co production 
cross section by protons on Te (0.63±0.15 mb for 1.85 GeV p II190I ) results in a strong 
disagreement with COSMO (which yields a much larger one) so that the allowed expo- 
sure period could be much longer. An independent activation measurement carried out 
by exposing a 5 1 . 1 g Te02 crystal to a beam of ~ 10^'^ 24 GeV protons has found a cross 
section of 1.6 mb instead than 2 mb as predicted by COSMO. Since however the effect 
of cosmic rays is expected to be dominant in the region of 1-2 GeV a new activation 
measurement with protons of 1.4 GeV is planned at the ISOLDE facility of CERN. A 
similar measurement is planned at Berkeley. Their reults will help in the final decision 
about the option to be chosen. 

12.3 Low background technology transfer to the crystal growth com- 
pany 

In either case, many critical operations in terms of radioactivity control must be accom- 
plished in China, while more will be required in the case of option (b). Therefore, it is 
mandatory that low background concepts will be transferred to SICCAS in order to get 
satisfactory results as for crystal radiopurity. 
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- The operations which precede crystallization involve the use of chemical prod- 
ucts. In particular, metallic tellurium must be transformed into Te02 powder: 
there are indications that U and Th impurities are added in this phase, that must 
be carefully checked, both in terms of procedure and of used reagents. In particu- 
lar, radiopure reagents must be selected by the CUORE collaboration (see section 
112.11) and provided to the Chinese company. 

- There are indications that the Pt crucible used for crystal growth contaminates the 
crystal with the isotope ^^°Pt, which is responsible for an alpha line at 3.25 MeV, 
very dangerous if degraded in energy. Therefore, new crucibles (possibly in Ir) 
should be developed in collaboration with S ICC AS. 

- The materials to be used for lapping and polishing (see section [TTT]), such as pads 
and powders, must be carefully checked and provided to SICCAS. This operation 
is even more critical in case of option (b) of section fTT2l Some experience has 
been gathered during CUORICINO construction, where special pads and powders 
down to 0.5 fim grain size were selected (see section imi) . This selection must 
be continued until a complete and satisfactory identification of the polishing ma- 
terials is accomplished. The materials will be delivered to China and used under 
the supervision of members of the CUORE collaboration. 

- In the event that option (b) of section fTT^ is selected, the copper frames, the ther- 
mistors, the heaters and all the tools necessary for the elementary module assem- 
bly must be shipped to China in special containers protecting them from exposure 
to radon and other contaminating agents. The Chinese personnel involved in the 
module assembly will be trained to fabricate and store the modules in the proper 
conditions in terms of radiopurity. 

12.4 CUORE assembly procedure 

We list below the main steps which precede and constitute the CUORE assembly pro- 
cedure, indicating the typical time necessary to accomplish each step, evaluated on the 
basis of the CUORICINO experience and discussions with SICCAS technicians. Under 
the general assumption that the total assembly procedure needs to last no longer than 
12 months, these indications will be used in the subsequent section to estimate the man- 
power needed, the amount of work that needs to be done "in parallel" and to satisfy the 
CUORE construction schedule. 

- Crystal production. Performed entirely in China. Estimated production rate: 
about 60 crystals/month. Total production time: 18 months. In this phase the 
crystals are produced with lapped surfaces, down to a ~ 15 fim roughness. 
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Crystal polishing. Performed in China or in the Gran Sasso laboratory. The expe- 
rience gathered with polishing crystals in Como for the CUORICINO experiment 
shows that 3 persons can polish 16 crystals/week with two polishers. In order to 
reach a higher rate of the order of 50 crystals/week, as required in order to remain 
on schedule with the complete assembly, 3 groups will have to work in parallel. 

Thermistor preparation. Performed in Berkeley (thermistor development and pro- 
duction) and in Como / Florence (bonding and characterization). The thermistors 
will be all ready for bonding and characterization 3 years after the start-up of 
the CUORE project. The bonding rate is about 20 thermistors/day. Two persons 
working in parallel with two bonding machines (already available) will bond 200 
thermistors/week. Total bonding time: 5 weeks. The characterization will be per- 
formed in parallel on samples bonded in the first week and should not require 
additional time. 

Heater preparation. Performed in Trento and in Milano. Eighteen months will 
be necessary to produce 2000 heaters, assuming a reasonable 50% yield. The 
bonding time will be 10 weeks at the same conditions mentioned in the previ- 
ous point. The characterization, that cannot be by sample, will have a rate of 50 
heaters/week. Therefore, about 1 year is needed for heater selection. Their pro- 
duction and characterization can take place during crystals growth, because no 
R&D is required for them. 

Heater and thermistor gluing. Performed in China or in the Gran Sasso labora- 
tory. The CUORICINO experience shows that 2 persons can glue thermistors and 
heaters with a rate of 15 crystals/week, using the special tools developed for this 
aim. 

Frame preparation. Performed in the Gran Sasso laboratory. Once the frames 
and the other copper elements for the typical 4-element module are machined 
and cleaned (see sections mi and fT2b. a preparation work preceding the module 
mounting must be performed, consisting in assembling the Teflon supports, the 
read-out pads and in wiring locally the frame. It is reasonable to assume that 2 
persons can prepare 20 frames (= 1 tower)/week. 

Crystal insertion in the frames. Performed in China or in the Gran Sasso labora- 
tory. The CUORICINO experience shows that 2 persons can insert crystals in 2 
planes/day. 

Connection of the gold wires. Performed in China or in the Gran Sasso laboratory. 
The CUORICINO experience shows that 1 person can wires in 1 plane/day. 
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- Assembly of the towers. Performed in the Gran Sasso laboratory. The CUORI- 
CINO experience shows that 2 persons working together can form a tower of 10 
elementary modules in about one week. 

- Assembly of the cube. Performed in the Gran Sasso laboratory. A reasonable as- 
sumption is that 2 persons can assemble the 19 towers into the CUORE cylinder 
in about 1 month, provided that special tools for towers manipulation and move- 
ment have been realized previously. Of course, they will need to work without 
anyone else around the cryostat, and this must be taken into account. 

- Cryostat wiring. Performed in a laboratory to be defined and in the Gran Sasso 
laboratory. In the case of CUORICINO, the preparation of a complete set of wires 
for the tower was made by 2 persons and took around 1 month. This operation can 
be made almost completely on the bench, before the construction of the CUORE 
cube, by several groups in parallel, so as to prepare all the wires in 6 months. An 
additional time of about 1 week/tower is required to lay the wires in the cryostat. 
This operation should be made by 2 persons working together. Unfortunately, 
this work cannot be done in parallel, since no more than 2 people can work on the 
cryostat simultaneously. Therefore, another 6 months should be considered for 
laying all the wires in the cryostat. 

- Final wire connection and check of the contacts. Performed in the Gran Sasso 
laboratory. Provided that the wires in the cryostat have been installed and that a 
clever method has been developed for connecting low thermal conductivity thin 
wires (i.e. special low radioactivity connectors), this job should be perfomed in 
about 1 month by 2 persons working simultaneously on the cryostat. Of course, 
this job is the only one that cannot be done in parallel with any of the other oper- 
ations. 

The proponents would like to stress that in the previous points no improvement in 
terms of simplification of the procedures was considered with respect to the CUORI- 
CINO experience. There are several innovations (larger crystals and consequent reduc- 
tion of the total channels, simplified methods for thermistor- absorber coupling, simpli- 
fied methods for wiring, etc.) which could speed up the assembly procedure. Some of 
these innovations are discussed in section [T51 

On the other hand, if the method for the abatement of the surface radioactivity pre- 
sented in section[T6lshould be adopted, the assembly procedure will present an increased 
complexity and might imply some delay. 




Figure 49: Gantt diagram illustrating the CUORE detector fabrication schedule. 
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12.5 Work breakdown and construction schedule for the CUORE 
detector 

On the basis of the previous section, it is possible to attempt a work breakdown and a 
construction schedule for the construction of the CUORE detector. The various items 
described in section [TT?I can be reported synoptically in a Gantt diagram (see fig-HU), 
taking care of the two following points: 

1) The assembly of the CUORE cube must be completed in about one year. Consid- 
ering that two years are necessary for the growth and preparation of the crystals, 
this means that CUORE will be cooled down for the first time three years after the 
growth of the first crystal. It is not feasible to span over three years the 4-detector 
module fabrication, along with crystal growth, since this would conflict with the 
time necessary to produce and clean the copper elements and with the Research 
and Development program. 

2) In order to concentrate the detector assembly in one year, and taking into account 
the requirements in terms of time and personnel for the realization of the vari- 
ous items, some of the tasks reported above need to be accomplished by more 
personnel working in parallel. For example, the operation of chip gluing will be 
performed by three groups of two people each, since only in this way it will be 
possible to accomplish this task in about three months. 

Once again, we stress that this construction schedule refers to the CUORE base- 
line. Innovations (consider for example Nd-based crystals or surface sensitive detec- 
tors) would inevitably modify the schedule in a way which is difficult to predict now, 
and might imply some delay. However, the reported analysis shows that the construc- 
tion of the CUORE detector in its conventional form is feasible in one year, but at the 
price of an extensive use of manpower. Notice that in some phases of the assembly, 
35 people must work simultaneously and full time on the detector preparation. Given 
the size of the CUORE collaboration and the limited availability of technical staff, it is 
hardly conceivable to realize this formidable task without a substantial use of external 
resources. 

12.6 General time schedule of the CUORE project 

As synoptically reported in the Gantt diagram of Figure |50| the general time schedule of 
the CUORE project requires approximately 5 years after approval. This summarizes the 
time sequence and duration of all the various tasks described in the previous sections. It 
should be stressed that the deadlines of many tasks are mandatory since often other jobs 
depend on them. However, most of the activities foreseen for the first 18-24 months 
can run in parallel without much interference. Therefore, most of the required R&D 
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activities could proceed during this period without resulting in sizeable delays on the 
project schedule. This is particularly true (as specified by the dashed filled area in 
fig.lSOb for the improvement in the surface (mainly copper) cleaning procedures. 
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Figure 50: Gantt diagram illustrating the CUORE project time schedule. 
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13 Off-line analysis 

The main goal of the off-line analysis is the extraction of the relevant physics informa- 
tions from the large amount of raw data recorded by the DAQ system. This task includes 
both the raw detector pulse analysis (e.g. amplitude evaluation, noise rejection, gain in- 
stability and linearity correction) aiming at producing reliable n-ple's (a proper number 
of parameters fully describing each bolometric pulse) and energy spectra, and the fol- 
lowing multi-dimensional analysis aiming at obtaining the sought physics results (e.g. 
PP(Oh') or DM interactions). The identification of the single background sources is also 
an important goal of the second-level analysis. 

13.1 First-level analysis 

As described in a previous section, when the output voltage of one detector exceeds the 
trigger threshold, the acquisition system records a number of converted signal samples. 
The acquired time window (~few sec) must fully contain the pulse development in order 
to allow an accurate description of its waveform. The existence of a pre-trigger interval 
just prior to the production of the pulse ("baseline") guarantees that a small fraction of 
the number of acquired samples can be used to measure the DC level of the detector 
(which can be directly correlated to the detector temperature). 

The fact that for each trigger an entire waveform is sampled and recorded implies 
that off-line analysis of a large body of data will require a significant effort. Neverthe- 
less, this effort will be justified because of the useful information that can be extracted 
from the signal waveform. The following are important goals for the (first-level or pulse) 
analysis (FLA): 

1. maximization of the signal to noise ratio for the best estimate of the pulse ampli- 
tude. This is accomplished by means of the optimum filter (OF) technique II191I : 

2. correction of the effects of system instabilities that change the response function 
of the detectors (gain "stabilization"); 

3. rejection of the spurious triggered pulses by means of pulse shape analysis; 

4. identification and rejection of radioactive background pulses by means of coinci- 
dence analysis. 

The OF technique is frequently used with bolometers to evaluate the amplitude of 
a signal superimposed on stochastic noise. This algorithm has proven to provide the 
best estimate of the pulse amplitude under general conditions. Relative to a simple 
maximum-minimum algorithm, this technique allows the evaluation of the signal ampli- 
tude with much higher efficiency resulting in an effective improvement of the detector 
energy resolution. The following information is needed to implement the OF technique: 
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the detector response function (i.e. the shape of the signal in a zero noise condition) 
and the noise power spectrum. Once these are known, the OF transfer function is easily 
obtained and used as a digital filter for the acquired pulses. The role of the OF transfer 
function is to weight the frequency components of the signal in order to suppress those 
frequencies that are highly influenced by noise. The amplitude of the pulse is then eval- 
uated using optimally filtered pulses. Effective pulse-shape parameters can be deduced 
from the comparison of the filtered pulses with the filtered response function (e.g. rms 
difference after proper time synchronization). 

In processing the data off-line, the following parameters are evaluated and recorded 
to disk for each digitized pulse (n-ple): 

1. the channel number i.e., the number of ADC channel that exceeded the trigger 
threshold; 

2. the absolute time at which the pulse occurred with a precision of 0.1 msec; 

3. the OF amplitude i.e. the amplitude of the optimally filtered signals 

4. the baseline, obtained by averaging a proper number of points from the pre-trigger 
interval. Since the detectors are DC coupled, it provides a direct measurement of 
the detector temperature at the creation of the signal; 

5. the signal rise and signal decay times; 

6. the pulse shape parameters, obtained by comparing the acquired pulse with the 
expected response function of the bolometer after OF or adaptive filters. A further 
powerful technique is based on the use of artificial neural networks (ANN); 

7. the pile-up fraction. Pile-up is usually efficiently rejected by pulse-shape analysis 
even if this technique can't identify the rejected pile-up events. In order to im- 
prove the pile-up rejection and quantitatively evaluate its rate (e.g. for short-time 
coincidence analysis), the Wiener-filter algorithm is implemented II192II . 

The next PLA step is the gain instability correction. The OF amplitudes are cor- 
rected to reduce or cancel the effects of system instabilities responsible for the variation 
of the ratio between the energy E deposited into a given crystal and the amplitude A V of 
the corresponding electrical pulse. According to equation (l29t of our very naive detec- 
tor model there are three instabilities that can modify the ratio A V/E (where V=V;,*G is 
the output voltage given by the product of the bolometer voltage and the electronics 
gain G): i) a variation in the electronic gain G, ii) a variation in the bias Wtoi and finally 
iii) a variation in the temperature T;, of the crystal. The electronic system is designed 
to guarantee a stability of G and WTot within 0.1% . It is however, much more difficult 
to maintain stability within 0.1% of the detector temperature on long time scales. At a 
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temperature of 10 mK this would require maintaining the temperature of 1000 crystals 
to an accuracy of 10 K for a period of several days. 

To overcome this problem, and as already mentioned in previous sections, a silicon 
resistor glued to each crystal is used as a heater to produce a reference pulse in the 
detector. It is connected to a high precision programmable pulser that produces a fast 
voltage pulse every few minutes dissipating the same amount of energy (E^e/) into the 
crystal each time. These voltage pulses mimic pulses produced in the crystal by par- 
ticle interactions and are used to measure the value of the ratio A V/E. Any variation 
of the amplitude of the reference pulse is due to variations of the A V/E ratio. The OF 
amplitude of the reference pulse is therefore used to measure, every few minutes, the 
actual value of A V/E while the baseline of the reference pulse provides the contem- 
porary measurement of the value of T. A fit is then used to obtain the values of A V/E 
as a function of temperature. Therefore, in this step of the off-line analysis, the OF 
amplitude of each pulse is corrected according to the given value of A V/E(T[t]) for 
the detector temperature at which the pulse has been generated. The effectiveness of 
this technique has been proven in the MiDBD experiment, where a typical temperature 
fluctuation over a day ranged from a few tenths to ~100 /iK . After correction, these 
fluctuations were reduced to less than 1 K [ 151 1 (Fig.l32b. 

Pulse shape analysis is very useful in rejecting spurious signals produced by micro- 
phonics and electronic noise. A confidence level is determined for each pulse shape 
parameter and for the rise and decay time of each pulse. Signals falling within these 
intervals are defined as "true" (or physical) pulses, while signals having one or more of 
their parameters outside of the relevant interval are rejected as noise. The use of more 
than one pulse shape parameter results in better reliability of the rejection technique. 

The linearization of the detector response is critically important for energy calibra- 
tion. The final step in data processing is the conversion of the OF amplitudes into energy 
values. The naive bolometer model previously used assumes linearity; several parame- 
ters depend however on the crystal temperature, rendering the corresponding equation 
non-linear. Accordingly, the relation between A V and E will be obtained periodically 
by the use of radioactive calibration sources. The ratio AV/ E will be measured for sev- 
eral gamma lines, and the data will be fit to the model previously described, but taking 
into consideration the fact that the bolometer resistance and the crystal heat capacity are 
temperature dependent. This will provide the calibration function of E as a function of 
AV, that will then be used to convert the OF amplitudes into energy values. 

Finally, the close packed array of the crystals will allow the rejection of events that 
left energy in more than one single crystal. This will be particularly useful in rejecting 
a very high energy gamma rays that enter from outside of the array. The small angle 
Compton scattering in a single crystal can mimic a double beta decay event, a dark mat- 
ter scattering, or a solar axion. The probability that this photon would escape the rest 
of the array without a second interaction is small. In addition, background events from 
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radioactivity within the structure of the array will also have a significant probability of 
depositing energy in more than one crystal. This will also be true for high and interme- 
diate energy neutrons. In the final stage of off-line analysis these coincidence events can 
be identified from the data which will contain the detector number, signal time, pulse 
energy, and pile-up parameter. 

13.2 Second-level analysis 

Since the main goal of CUORE is the search for /?/?(0i/) of ^^°Te , the same will also 
hold for the second-level analysis (SLA). In CUORICINO, most of the tools of this 
analysis are powered by the experience gained by the Milano group during the previous 
experiments on /5/9(0z/) of various nuclides. Since many institutions participating in 
CUORE have a similar experience, we are planning to merge the various analysis tech- 
niques of the different groups to apply them to double beta decay processes in order to 
obtain a single powerful standard analysis procedure. Nevertheless, it will be essential 
to be able to guarantee a safe cross-check of the results, so that independent analyses 
from the various groups will be encouraged. 

Detection efficiencies for the various transitions will be evaluated by means Monte 
Carlo simulations specifically designed for the CUORE detector. An updated collection 
of the most recent calculations for nuclear matrix elements will be kept and updated 
(according to the various nuclear models). That will be essential to derive reliable results 
for the effective neutrino mass parameter, the right handed current parameters, or the 
majoron coupling parameter. 

Physics analysis for Dark Matter and axion searches in CUORE will be completely 
designed and realized in the framework of CUORICINO. In order to improve the noise 
rejection near threshold an Artificial Neural Network (ANN) approach will be improved 
(the possibility of an ANN trigger is under investigation). The effectiveness and relia- 
bility of the rejection techniques based on the selection of multi-dimensional confidence 
intervals (strongly depending on the suitability of the chosen parameters) will be tested 
in CUORICINO and methods to calibrate the rejection efficiency will be defined. Time 
stability of threshold parameters will be a major issue. 



14 SIMULATION AND PREDICTED PERFORMANCES 



140 



14 Simulation and predicted performances 

The goal of CUORE is to achieve a background rate in the range 0.001 to 0.01 coun- 
ts/(keV-kg-y) at the /3/5(Oz/) transition energy of ^'^°Te (2528.8 keV). A low counting 
rate near threshold (that will be of the order of ~5-10 keV) is also foreseen and will 
allow CUORE to produce results for Dark Matter and Axions searches. In Sec. ll4.ll we 
present a very conservative evaluation of the background attainable with CUORE; this is 
mainly based on the state of the art of detector design and of radioactive contaminations. 
This rather pessimistic approach in background evaluation is the only one that presently 
guarantees a reliable prediction. However CUORE construction will require about five 
years (sec. 112. 61) : there will therefore be time enough for R& D dedicated to background 
reduction as is foreseen in most next generations experiments. 

14.1 Background simulations 

Radioactive contamination of individual construction materials, as well as the laboratory 
environment, were measured and the impact on detector performance determined by 
Monte Carlo computations. The code is based on the GEANT-4 package; it models the 
shields, the cryostat, the detector structure and the detector array. Even smallest details 
of the detector apparatus (copper frames, screws, signal wires, NTD thermistors, etc.; 
Fig. BTt and of the cryogenic setup were taken into account. 

Table 16: Bulk contamination levels (in picograms per gram) used in the simulation for 
Te02 , copper and lead. 



Contaminant 




238^ 40 


K ^^"^Ph 


60Co 


TeOs 


0.5 


0.1 


1 10 /iBq/kg 


0.2 /iBq/kg 


copper 


4 


2 


1 


10 yuBq/kg 


Roman lead 


2 


1 


1 4 mBq/kg 





16Bq/kglead 


2 


1 


1 16Bq/kg 






It includes the propagation of photons, electrons, alpha particles and heavy ions 
(nuclear recoils from alpha emission) as well as neutrons and muons. For radioactive 
chains or radioactive isotopes alpha, beta and gamma/X rays emissions are considered 
according to their branching ratios. The time structure of the decay chains is taken into 
account and the transport of nuclear recoils from alpha emissions is included. 

The considered background sources are: 

1 . bulk and surface contamination of the construction materials from the ^^^f/, ^^^T/i 
chains and ^^K and ^^°P6 isotopes; 




Figure 51: Details of the CUORE single modules accounted for in the Monte Carlo 
simulations 

2. bulk contamination of construction materials due to cosmogenic activation; 

3. neutron and muon flux in the Gran Sasso Laboratory; 

4. gamma ray flux from natural radioactivity in the Gran Sasso Laboratory; 

5. background from the /3/5(2z/) decay. 

In the next sections these background sources and their contribution to CUORE 
background are discussed. It will be shown that, in a very conservative approach where 
the CUORE-tower mechanical structure is assumed identical to the structure used in 
CUORICINO, considering the worst possible condition for bulk contaminations (i.e. 
all the contamination equal to the present 90% C.L. measured upper limits (Table fT6b) 
and assuming for surface contamination a reduction of about a factor ten with respect 
to CUORICINO, the CUORE background will be ~0.007 counts/(keV-kg-y) at the 
PP^Ou) transition and ~0.05 counts/(keV-kg-d) near threshold. 
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14.2 Bulk contaminations 

The main contribution to background from bulk contaminants comes from the cryostat 
structure (cryostat radiation shields), the heavy structures close to the detectors (the 
copper mounting structure of the array, the Roman lead box and the two lead disks 
on the top of the array) and from the detectors themselves (the Te02 crystals). The 
radioactivity levels used in the computations for these materials are given in Table fT6l 
These levels are nearly equal to the best upper limits obtained for the radioactive content 
of these same materials as shown in Table flTl 

All the values reported in this table are 90% upper limits. Indeed in all cases no 
evidence of the presence of radioactive contaminants was obtained for the material ex- 
amined. The contamination levels of Te02 reported in Table fTTI were obtained from the 
data collected in the first run of CUORICINO with the 5x5x5 cm^ crystals. Copper 
||T93 1 and lead HI 941 contaminations were determined through low activity Ge spectrom- 
etry. These data represent the best evaluation reported in literature for the radioactive 
contamination of Te02, copper and lead. From a comparison between Table fTTI and 
Table fT6l it is clear that the contamination levels assumed for the CUORE simulation 
are perfectely compatible with the presently available upper limits. The only exception 
is that of ^°Co in Te02 crystals, whose presence as a cosmogenic contamination of the 
crystals will be discussed later ("sec. 114.41) . 

Concerning small (mass) components of the CUORE detector, which we have not 
considered so far in the simulation (e.g. NTD Ge thermistors. Si heaters, glue layers, 
pins, wires and soldering material,etc.), a very careful selection according to their con- 
tamination is planned for the future. In some case we already have upper limits on their 
contamination that make their contribution to the CUORE background negligible; in 
other cases further measurements will be required. 

The results of the Monte Carlo simulations (see Fig. l52b using the contamination 
levels discussed here are given in Table flSl for the PP{Ou) decay and the low energy 
(10-50 keV) regions. The assumed threshold is 10 keV and only values obtained after 



Table 17: Available 90% C.L. upper limits for bulk contaminations of Te02, copper and 
lead (levels in picograms per gram if not differently indicated). 



Contaminant 


method 


232Th 


238^ 


40k 


210pb 


60Co 


Te02 


bolometric 


0.7 


0.1 


1. 


100 /uBq/kg 


1 ;uBq/kg 


Copper lfT93l 


Ge diodes 


5.6 


2 


0.3 




10 yuBq/kg 


Roman lead 


Ge diodes 


50 


30 


2 


4 mBq/kg 




Low act. lead lil94J 


Ge diodes 


3.4 


2.7 


1.7±0.3 


23. 4± 2.4 Bq/kg 


18± 1 /iBq/kg 
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Figure 52: Simulated spectra for bulk contaminations of the Te02 crystals, the Copper 
structure, the Roman Lead shield and the outer Lead shield. Each spectrum is obtained 
by summing the simulated anticoincidence spectra of all the CUORE detectors. 

requiring an anti-coincidence between detectors are indicated. These values have to be 
considered as upper limits on the possible contribution of bulk contaminations to the 
CUORE background; they prove that with the already available materials, background 
levels lower than ~4x 10^'^ counts/keV/kg/y in the /3/3(Oz/) decay region (see Fig l53b 
and ~3x 10^^ counts/keV/kg/d in the low energy (10-50 keV) region (see Figl54b are 
assured. 

14.3 Surface contaminations 

Surface contaminations contribute to the background only when they are localized on 
the crystals or on the copper mounting structure directly facing them. As learned in 
MiDBD and CUORICINO (sec.|6l), the presence of even a very low level of radioactive 
impurities on the surfaces of the detectors can produce a non negligible contribution to 
the /3/3(0z/) background level. Unfortunately this kind of contamination, that is impor- 
tant mainly for detectors without a surface dead layer (as is the case of bolometers). 



14 SIMULATION AND PREDICTED PERFORMANCES 



144 



X :9 

o o 
^ m Z 
^ T3 c/> 

o §: E S E 

CD O O ^ ^ 

I- O DC O C/D 



J I L 



\ a 




o o 
^ [i/§5I/A3>I/sjuno3] 



Figure 53: Double beta decay region of the simulated spectra for bulk contaminations 
of the Te02 crystals, the Copper structure, the Roman Lead shield and the outer Lead 
shield. Each spectrum is obtained by summing the simulated anticoincidence spectra of 
all the CUORE detectors. 
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Figure 54: Dark matter region of the simulated spectra for contaminations of the 
Te02 crystals, the Copper structure, the Roman Lead shield and the outer Lead shield. 
Each spectrum is obtained by summing the simulated spectra of all the CUORE detec- 
tors after the anticoincidence cut. 
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counts/keV/kg/y 


1.6x10-3 


1.5x10-3 


7.0x10-^ 


3.8x10-3 


dark matter region 










counts/keV/kg/d 


2.3x10-2 


9.6x10-^ 


5.0x10-^ 


2.4x10-2 



Table 18: Computed background in the PP^Ou) decay and in the low energy regions for 
bulk contaminations in the different elements, the Cu structure accounts for the detector 
mounting structure and the 50 mK shield. 

is poorly studied. Lacking any data coming from direct measurements of the typical 
impurity levels present on Te02 and copper surfaces, we rely on the results obtained in 
MiDBD and CUORICINO. 

These indicate that both a surface contamination of the crystals and of the copper 
surface of the mounting structure is present. In the case of Te02 , we know that the 
contamination is mainly ^^^u and that its presence is strictly connected to the kind of 
surface treatment undergone by the crystals. In the CUORICINO 5x5x5 cm^ crystals 
this contamination produces a background counting rate in the /9/5(0z/) region of about 
(4± 3) X 10-2 counts/keV/kg/y (table fT4ll. Concerning the copper surface contamina- 
tion, we still miss a clear indication of its origin and identity (i.e. we are not yet able to 
distinguih between contributions generated by different sources), but we have an eval- 
uation of the /3/3(Oz/) background counting that can be ascribed to them (whatever they 
are, ^^^U or 232Th): ~0.1± 0.05 counts/keV/kg/y. 

According to a Monte Carlo simulation of the CUORE detector, based on the CUORI- 
CINO contamination levels, we obtained a background counting rate in the PP{Ou) re- 
gion (after the anticoincidence cut) of about 1.6x 10-^ counts/keV/kg/y and 5.8x 10-^ counts/keV/kg/y 
for the Te02 crystals and the copper structure respectively |20l The goal of CUORE is 
to reduce the surface contribution by a factor at least 20 with respect to this evalua- 
tion, obtaining a background coming from surfaces of about 3x 10-^ counts/keV/kg/y. 
Improvements of a factor of ~ 1.5-2 are expected simply by possible reductions of the 
copper mounting structure dimensions (the Monte Carlo simulation we used so far refers 
to an identical mechanical structure for the CUORE and CUORICINO single detector 
modules). A reduction by a factor of at least ten of the copper (Te02 ) surface contami- 
nation is therefore the first milestone of CUORE. 

The expected shape of the background produced by surface contaminations in CUORE 
(after a reduction by a factor 20 of both the crystal and copper surface contaminations 
but assuming the same mechanical structure of CUORICINO) is given in Figures |55] 
and |56l The corresponding background contributions in the two regions of interest are 
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summarized in Table fT9l 

As discussed in Sec. 112.11 the MiDBD crystal surfaces were contaminated, prob- 
ably during the polishing procedure, due to the use of highly contaminated powders. 
Polishing powders with a much lower radioactive content are however commercially 
available 1 195| and have already been used for CUORICINO (tab. [T5[ obtaining a clear 
improvement in the surface contamination. Complete over all control of the surface 
treatments undergone by the crystals after growth, and the use of radiopure substances 
should guarantee an even better result than that obtained in CUORICINO. 

A similar situation holds for copper. In MiDBD and in CUORICINO the copper 
surfaces were treated with an etching procedure 119611 optimized in order to reduce 
impurities on surfaces before the sputtering process. 

This procedure significantly improved the surface quality of copper and reduced its 
surface contamination. However it was not optimized from the point of view of back- 
ground. The use of low contaminated liquids (water and inorganic acids are available 
with ^^^U and ^^^Th contamination levels lower than 0.1 pg/g) in a low background 
environment will allow a sizeable improvement of the surface contamination of copper. 

This (i.e. similar to CUORICINO but with much more care to prevent radioactive 
contamination) and other improved cleaning procedures (sec. flSb are under study and 
will represent a consistent part of the CUORE R&D activities during the first two years 
after approval. Fast diagnostic methods to measure the surface contamination levels of 
Copper and determine the effectiveness of the adopted surface cleaning procedures are 
also under investigation. Preliminary results indicate that high resolution ICPMS could 
satisfy our requirements 119711 . Sensitivities achievable with this technique for ^'^^U 
and ^^^Th are in the range of 10~^^~^^ ppt of solution under investigation. Since the 
maximum solution concentration can be 2x lO^'^, the expected sensitivity on the solute 
contamination (Copper) would be in the range 1-10 ppt of ^^^U and ^^^Th, i.e. enough 
for our requirements. From the above discussion in fact, the surface contamination 
measured in CUORICINO (assuming a 1 fim copper layer) is of the order of ~1-10 ppb 
of ^'^^U and ^'^^Th. Of course radio-purities comparable with the instrumental sensitivity 
are required for liquids (mainly acids and water) to be used in the solution preparation 
(i.e. copper surface removal). Liquids satisfying this requirement are however already 
commonly available. 

14.4 Cosmogenic contribution 

Cosmogenic activation is produced by cosmic rays when the crystals are above ground 
during fabrication and shipping of the crystals from the factory to the underground lab- 
oratory. To determine the activation we used COSMO, a code based on computed cross 
sections, to estimate the type and amount of radionuclei produced by cosmic rays on 
Te02 II198I . The radionuclei produced by the activation of tellurium are mostly tel- 
lurium isotopes (A=121, 123,125,127) as well as ^'^'^Sh, ^^^Sb, ^°Co and tritium, these 
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Figure 55: /9/3(0z/) region of the simulated spectra for the surface contaminations of the 
Te02 crystals and of the Copper structure. Each spectrum is obtained by summing the 
simulated spectra of all the CUORE detectors aftre the anticoincidence cut. 




Figure 56: Dark matter region of the simulated spectra for the surface contaminations 
of Te02 crystals and of the Copper structure. Each spectrum is obtained by summing 
the simulated anticoincidence spectra of all the CUORE detectors. 
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Cylindrical structure 



Simulated 



232 



Th 



238 



u 



60, 



Co 



Sum 



TeOs 
Cu Box 
Cu Bars 
Cu Frames 
Cu 50mK Shield 
Roman Pb Shield 
Sum 



212±23.6 
393±46.3 
29.7±4.08 
163±21.5 
31.9±2.21 
604±24.4 
1430±61.5 



1200±19.5 



14.5±5.11 
84.8±28.3 
27.1±5.12 
233±33.7 



360±44.6 1200±19.5 



1420±31.0 
4784±54.3 
56.8±6.55 
397±39.9 
31.9±2.21 
604±24.4 
2390±74.5 



Table 19: Computed background (after the anti-coincidence cut) in the (3(3{0u) energy 
region for bulk contaminations in the different elements for cylindrical structure in units 
of 10-6 c/keV/kg/y. 



Element 


Contamination 


Contribution to thePP{Ou) region 




Bq/cm^ 


c/keV/kg/y 


Te02 


1.9E-08 ± 9.8E-09 


1.6E-02 ± 8.4E-03 


Copper 


4.9E-08 ± 1.9E-08 


5.8E-02 ± 6.6E-03 


Total 


3.4E-()8 ± 1.3E-()8 


7.4E-()2 ± l.lE-()2 



Table 20: Estimated upper contribution to the CUORE /3/3(0z/) region from surface con- 
taminations obtained by using the surface contamination levels evaluated for CUORI- 
CINO and assuming an exponential density profile with A=l /im for Te02 crystals (^^^U 
) and A=5 //m for Copper (^^^u and ^^^Th ). 
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last three being of more concern because of their long half-life ( Sb: beta decay of 2.7 
years, end-point energy of 767 keV, ^°Co: beta decay of 5.27 years end-point energy 
2823 keV and ^H: beta decay of 12.3 years, end-point energy of 18 keV). 

We have studied, toghether with SICCAS who produces the crystals, and on the basis 
of our present knowledge of the cosmic rays production rates, a possible time schedule 
for the crystal growth and shipping to Gran Sasso that can guarantee the required low 
level of ^°Co (see Fig l57t . In the case of CUORE, the control on crystal production 
will be severe. A single 750 g crystal can be grown in about two months, while to 
grow the 988 CUORE detectors will require 18 months. Once grown the crystals will 
be shipped to Italy and stored underground, therefore their total exposition to cosmic 
rays could be limited to about 4 months. The total induced activities remaining after 
2 years underground have been estimated [1981 and the consequent contribution to the 
detector counting rate was deduced by a MC simulation. The radionuclei that contribute 
to background throught their decay are: 

• the long living ^°Co isotope in the /3/3(Oz/) region, with an activity of ~0.2 /iBq/kg, 
while a minor contribution is due to the isotopes ^^^'^Ag and ^^''Sb whose activity 
is 4 times lower and fast decreasing with time; 

• in the dark matter regions the long lived nuclei of '^H and ^^^Sb with an activity 
of ~7 yuBq/kg for the former and of ~15 /iBq/kg for the latter. 

The influence of ^°Co in CUORE background was already considered in the evaluation 
of the contributions due to bulk contaminations of the crystals, while the contribution 
of '^H and ^^^Sb to the dark matter region is completely negligible if compared to the 
intrinsic background from all other sources. 

As already mentioned in section 112.21 a recent experimental determination of the 
^°Co production cross section on Te by 1.85 GeV protons has been obtained at LNBL | l^UH . 
The measured value (0.63±0.15 mb) is in strong disagreement with the value used by 
the COSMO program. A similar disagreement, even if at a lowerextent, has been found 
on the basis of an activation measurement carried out at CERN at an energy of 24 GeV, 
which is however less important for cosmic ray activation. The contribution to CUORE 
/3/3(Oz/) background by cosmogenic ^°Co discussed above, could be therefore overesti- 
mated and the allowed exposure period for Te02 crystals could be consequently much 
longer. Further measurements of the ^°Co production cross section are therefore of cru- 
cial importance and we plan to have, in the near future, other direct measurements by 
irradiating Tellurium samples with proton beams of suitable energies. 

14.5 Underground neutron, fi and 7 interactions 

As mentioned earlier, neither contributions from underground cosmic muons nor neu- 
trons have been taken into account in detail in the estimation of the background. How- 
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ever, the following simplified arguments will serve to have an approximate idea of their 
contribution. The depth of the LNGS (3500 m.w.e) reduces the muon flux down to 
~2x 10~^ cm~^s~^, but a further effective reduction could be obtained with the use of 
an efficient (99.9%) active veto for muons traversing the CUORE setup in order to tag 
possible events associated with them. The muon-induced contribution to the background 
is therefore expected to be negligible. 

On the other hand, the heavy shieldings surrounding the CUORE detector will sub- 
stantially reduce the event rate due to environmental radiation of various origin (neu- 
trons and photons), environmental radioactivity (natural decay chains U/Th, ^^°Pb, ^°K, 
. . . ), as well as muon interactions in the surroundings rock or in the shielding itself. 

Neutrons may constitute a worrisome background for the dark matter experiment 
because for appropriate neutron energies (few MeV) they can produce nuclear recoils 
(<100 keV) in the detector target nuclei which would mimic WIMP interactions. Sim- 
ple kinematics implies that in the case of tellurium, neutrons of 1(5) MeV could elas- 
tically scatter off tellurium nuclei producing recoils of energies up to 31 (154) keV. In 
general, one considers neutrons of two origins: from radioactivity in the surroundings 
or muon-induced. Depending on the overburden of the underground site (i.e., depend- 
ing on the muon flux), muon-induced neutrons are produced, at lesser or greater rate, 
both inside and outside the shielding. They are moderated (according to their energies) 
by the polyethylene/lead shield (when produced outside) or tagged by the muon veto 
coincidence (when produced within the passive shielding). 

In the case of external neutrons (from the rocks, from fission processes or from 
(n,a) reactions, as well as neutrons originated by muons in the walls of the under- 
ground site), the environmental neutron flux has been measured in LNGS. The result 
is of ~1 X 10~^ cm~^s~^ for the thermal component, ~2x 10~^ cm~^s~^ for the epither- 
mal and ~2x 10~^ cm~^s~^ for energies over 2.5 MeV. They are fairly well moderated 
by the polyethylene and eventually absorbed or captured. 

We have carried out a Monte Carlo simulation of the propagation of neutrons through 
the 10 cm thick borated polyethylene shield of CUORE. The result is that the neutron 
induced event rate on the entire energy range (from threhsold to 10 MeV) is much 
lower than the contribution due the bulk contamination of crystals. Neutrons produced 
by muon interaction inside the shielding materials are very scarce and they can be ef- 
ficiently tagged by the muon veto. We have estimated that for the LNGS muon flux 
(2.5x 10~^ /x/(cm^ s)), muons would produce in the CUORE shielding of polyethylene 
(10 cm) and lead (20 cm)), about ~0.04 neutrons/day in the polyethylene shield and 
~25 neutrons/day in the lead shell. So, independent of the mechanism used to reject or 
tag the events associated to neutrons, their rather small number is expected to play a sec- 
ondary role in the total background compared with other main sources of background. 

A preliminary evaluation of the influence of the environmental 7 background in 
Gran Sasso resulted in a negligible contribution for the PPiOu) region and a contribution 
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Figure 57: Simulated spectrum for the assumed °°Co bulk contaminations of Te02 crys- 
tals. The spectrum is obtained by summing the simulated anticoincidence spectra of all 
the CUORE detectors. 
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similar to that of bulk contaminations for the dark matter region. 

A more complete and detailed study of the background rates from external sources 
for CUORE is underway and will be used for the optimization of shieldings and muon 
veto. 

14.6 Two neutrinos double beta decay background 

Using the present upper limits for the ^^°Te I3l3{2u) half-life, the unavoidable back- 
ground produced by the P/3{2i') decay in the dark matter region is lower than 10~^ co- 
unts/(keV-kg-y) and is completely negligible in the PP{Ou) region. This is true because 
of the relatively good energy resolution of the Te02 bolometers. 
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15 Innovations in the CUORE structure 

The CUORE structure described up to now is based on the results obtained in CUORI- 
CINO and on the experience gathered with the test-runs performed in the hall C cryostat 
on the 5x5x5 cm^ crystals. We have been very conservative in the design of every 
CUORE element. However, we foresee the requirement for a significant activity of Re- 
search and Development during the next few years, in parallel (but with lower priority) 
with the realization of the CUORE baseline structure outlined so far. In fact, there are 
indications that several aspects of the CUORE detectors can be improved substantially 
with respect to the present performance (for instance, single element reproducibility) 
and some of them need to be improved dramatically (mostly, the background level due 
to surface contamination that could spoil CUORE, and a large part of its scientific value, 
if not controlled adequately). 

In this section, we outline the R&D work planned for CUORE and discuss the possi- 
ble upgrading of the present structure. We would like to stress that all this R&D activity 
can be carried out in Hall C without any interruption of the presently running CUORI- 
CINO experiment. 

15.1 Innovation in the single module 
15.1.1 Optimization of the thermistor sensitivity 

As described in section ITTTl our current NTD thermistors (belonging to series #31) have 
a To value of about 3.0 K. This parameter has never been optimized and was obtained 
only in the attempt to get a tolerable resistance value between 10 and 15 mK (not higher 
than a few hundreds of M^7). In the two versions of the MiDBD arrays, and in the 
CUORICINO experiment, we had however to heat the refrigerator base temperature up 
in order to achieve these resistance values. The naturally achieved base temperature 
would have lead to too high resistances, of the order of few GVl. There are two reasons 
why lower Tq (and therefore lower resistivities at equal temperatures) would probably 
imply better detector performance, in spite of the lower thermistor sensitivity: 

• it would be possible to operate detectors either at lower temperatures (therefore 
with higher signals) or at the same temperature as now, but with lower resistances 
(therefore with lower noise); 

• it is well known that lower To thermistors have better electron-phonon coupling 
(see for instance [il99in . favoring therefore a more efficient heat transmission to 
the thermistor electrons with consequent higher pulse amplitudes. 

An R&D program, lasting about one year, will therefore define the optimum To for 
CUORE thermistors. Several thermistor species with To's ranging from 2 through 3 K 
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will be produced by the Berkeley group and test detectors will be produced in order to 
select the best doping level in time for the production of the final thermistor batch. 

15.1.2 Innovative thermal couplings in detector construction 

A set of glues will be selected to improve the coupling between NTD Ge thermistors 
and energy absorbing crystals. Methods will be studied to deposit controlled amounts of 
glue in a reproducible way, in the form of spots and/or a uniform layer. Measurement of 
the low temperature thermal coupling provided by the various glues will be performed, 
in order to define the most conductive coupling. Of course, all these cryogenic tests 
must be accompanied by a careful analysis of glue radioactivity, so as not to introduce 
dangerous radio-impurities. Special NTD Ge thermistors with "pedestals" and "legs" 
will be developed in order to reduce the surface touching the crystals. This reduction 
will allow the use of a glue layer at the interface in spite of the large difference in thermal 
contraction coefficients between glues and crystalline solids. 

In addition, the proponents intend to perform tests in a thin film laboratory on eu- 
tectic coupling between NTD Ge thermistors and energy absorbing crystal. The most 
promising couplings will be realized in small scale prototype detectors and compared 
with the standard 9-araldyte-spot coupling, mentioned in section 1531 

Innovative methods to hold and couple the absorbing crystals to the thermal bath 
will be studied, such as special TEFLON supports, nylon threads, simple gluing etc., 
with and without additional calibrated thermal links. Prototype small scale detectors 
will be studied. 

Finally, the proponents will perform tests of real- size CUORE modules in LNGS 
hall C dilution refrigerator with the innovative couplings determined in the prototype 
small scale detectors. If new couplings, and a related procedures for their series realiza- 
tion, will be defined within three years since the start of the CUORE project, they will 
be introduced in the final CUORE structure. 

15.1.3 Modification of the basic crystal size 

There is clear experimental evidence that a crystal mass increased by about a factor 
two should not spoil basic detector performance. At present, the energy resolution at 
2.6 MeV is higher by a significant factor (typically a factor ~5) than that expected from 
the detector noise. Therefore, an increase of the detector mass by a factor two should 
affect the intrinsic detector resolution because of the consequent reduction of the sig- 
nal amplitude, but should influence only marginally the energy resolution in the Double 
Beta Decay region. Of course, this argument does not apply to the search for Dark Mat- 
ter and Solar Axions, where threshold, and therefore signal-to-noise ratio at low ener- 
gies, is the crucial parameter. However, in a more Double-Beta-Decay-oriented version 
of the CUORE project, an increase of the crystal size (with equal total mass) could bring 



15 INNOVATIONS IN THE CUORE STRUCTURE 157 



Table 21: Most favorable /3/3{0h') candidate nuclei. 
Compound Isotopic abundance (%) Transition energy (keV) 



^^CaFa 0.0187 4272 

^^Ge 7.44 2038.7 

i°°MoPb04 9.63 3034 

ii6CdW04 7.49 2804 

i30TeO2 34 2529 

isoNdFg 4.64 3368 



some advantages, mainly of practical and economical nature (fewer electronic channels, 
shorter and easier assembly operations). The proponents have in mind in particular ele- 
mentary modules based on 6x6x6 cm^ crystals, whose mass is larger by a factor 1.73 
than present crystals' one. A possible configuration consists of 16 towers arranged in a 
4x4 matrix. Each tower will be a stack of 9 modules. As in the present configuration, 
each elementary module will consist of 4 crystals held in common frames. The total 
mass would be very similar to the baseline CUORE structure (namely, 756 kg, assum- 
ing the same density that provides 750 kg for the present structure), and an almost cubic 
shape for the array would be preserved. Another approach consists of giving up the 
4-detector arrangement of the single module (which might not be testable in the present 
refrigerators because of its size) and to move to stacks of single crystals. Even in this 
case, the favored arrangement is a 8x8x9 (height) three-dimensional crystal matrix. In 
both solutions, the total number of channels would be 576 instead of 1000. 

Negotiations for the order and purchase of 6x6x6 Te02 crystals are in progress. 
Two such larger crystals have been already ordered and should be tested in Gran Sasso 
in early 2004. 

A final decision on the crystal size should be taken not later than 1 8 months after 
CUORE start: this should leave sufficient time for the series production of the crystals. 
That implies a hard R&D work on 6x6x6 Te02 crystals in the first year. 

15.2 Possible insertion in CUORE of crystals containing compounds 
of other nuclear candidates for PP{Oiy) 

As pointed out before many compounds containing /9/5(0z/) candidate nuclei could be 
inserted in form of single or multiple crystals in the structure of CUORE. Some of these 
compounds which look as particularly favourable from the thermal and nuclear point of 
view are listed in Table 1211 

All these compounds have been tested as absorbers of thermal detectors. These tests 
were successful, but some difficulty was found for NdFs and NdGaOa, due to the para- 
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magnetic properties of Neodimium. A particularly good candidate would be a crystal 
of ^^CaF2. In fact a "scintillating bolometer" where light and heat were simultane- 
ously detected was already operated |200|. Operation of these detectors in coincidence 
mode between heat and light pulses would be very effective in the suppression of the 
PP(Ou) region background due to degraded alpha particles. Unfortunately the very low 
isotopic abundance of '^^Ca makes this experiment fiscally impossible, unless new pro- 
cesses for enrichment of this isotope be discovered. Other scintillating bolometers like 
those with CdW04 or MoPb04 can obviously be envisaged. We would like to specify 
here that also Te02 crystals have been found to emit light [201 J. Their yield is however 
very low (about 50 eV per MeV) which make them inefficient for the operation of a 
scintillating bolometer. This approach would be in any case much less effective than 
the active shield with Ge or Si bolometers, which will be considered later. This shield 
would in fact suppress not only the surface activity due to alpha particles, but also the 
surface activity due to beta decays. 

A dedicated R%D work is foreseen for the Neodimium case, since ^^°Nd is probably 
the best candidate to search for /?/5(0z/) decay |Qj, thanks to the high transition proba- 
bility due both to favorable nuclear matrix elements and to the high transition energy 
(3.37 MeV). Unfortunately, ^^°Nd natural isotopic abundance is low (5.6 %) and no 
dielectric diamagnetic compound of Nd is known to exist. The first problem can proba- 
bly be solved thanks to innovative techniques of isotopic enrichment, such as ICR (Ion 
Cyclotron Resonance) method [202|, while the second one poses a challenge for the 
bolometric technique. The main difficulty, as already mentioned, is due to the magnetic 
properties of Nd: all the Nd compounds compatible with growing large crystals present 
a magnetic ordering implying inevitably high specific heat at low temperatures. How- 
ever, the chance to realize a working Nd-based bolometer relies on our ignorance. The 
thermal models which attempt to explain how a large-mass, phonon-mediated bolome- 
ter works are not completely convincing. In particular, there are indications that a non- 
thermal component of the signal, due to "high" energy phonons, may be important (see 
section |5l last paragraph). In this case, the specific heat could be a scarcely relevant 
parameter: the right question would be rather how long a non-thermal phonon lives in 
a crystal with a magnetic ordering and which chance it has to transmit its energy to the 
NTD thermistor. We have decided to follow a strictly experimental approach and to 
plan bolometric tests on commercially available crystals containing Nd. We have per- 
formed tests on existing NdFa (a paramagnetic material with huge Shotckey peaks at 
about 70 mK) and NdGaOs (which presents a more promising antiferromagnetic order- 
ing below ~1 K) crystals in the hall C dilution refrigerator. In particular tests carried 
out with two cubic crystals of NdGaOs of 1.3 and 3 cm side have not allowed to reach, 
after weeks of cooling, temperature lower than ~40 mK even if with different time pro- 
files. These temperatures would not be sufficient to produce bolometers with acceptable 
energy resolution. 
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In about two years from the start of the CUORE project a final answer on the fea- 
sibility of a Nd-based bolometer must be given. In the event of a positive response, 
a part of the CUORE array (probably not more than ~ 10 % in terms of mass) could 
consist of Nd-based crystals. The Nd section of the array could be even more sensitive 
to Majorana neutrino mass than the Te section, in spite of the lower mass. In this case, 
CUORE would be a formidable two-fold experiment, capable to cross check within 
itself a positive result in the search for /3/3(Oz/) decay. 

15.3 Development of surface sensitive Te02 elements 

Perhaps, the most important objective of the CUORE R&D program consists of devel- 
oping methods to control and to reduce at a negligible level the radioactive background 
coming from surface events. This background is presently observed in the large mass 
bolometers of the CUORICINO experiment and was observed in the previous 20 de- 
tector arrays. Monte Carlo simulations (see section [1431 show that the dominant back- 
ground contribution in the PP^Ou) region of interest (around 2.5 MeV) is originated by 
passive elements which surround the detectors (in particular, the copper supports for the 
crystals, which work also as a 5 mK thermal bath for the detectors). This background 
is to be ascribed to surface contaminations in Th and U, which emit energy-degraded 
alpha particles able to deposit their energy at the detector surfaces, populating the 2000- 
4000 keV spectral region. The presently designed CUORE configuration would not 
allow a substantial improvement of this background (expected to be of the order of 
0.1 counts/(keV kg y)) unless sizeable improvements in the reduction of the surface 
contributions would be obtained. 

The objective of a specific R&D program has to be therefore the realization of ele- 
ments whose surface background is reduced down to 0.001-0.001 counts/(keV kg y) in 
order to satisfy the CUORE sensitivity requirement. This result can be achieved in three 
ways: 

• to improve dramatically the quality of the surface treatment of copper (this option 
is not discussed here, see section [161 and should be pursued in parallel with the 
two following ones; 

• to realize single CUORICINO and CUORE elements with a radically innovative 
structure, in which the copper surface facing the detectors is reduced by a factor 
at least 100-1000; 

• to realize bolometers able to distinguish an event originated at the surface from 
one originated in the bulk. 

As far as the second item is concerned, the R&D program consists of identifying 
materials which can replace copper in forming the main part of the detector holder 
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Figure 58: Scatter plot of the pulse amplitudes from Te02 thermistor versus the pulse 
amplitudes from shield thermistor. 



(best candidates are plastics, such as TEFLON, nylon or acrylic materials, achievable 
with very high radiopurity). Detector holders based on these new materials have to be 
developed, first for reduced-scale detectors, and secondly, after selection of the most 
promising solutions, for modules with the final size foreseen for the CUORE experi- 
ment. 

As far as the third item is concerned, the basic idea consists in the realization of 
active shields in the form of thin, large-area, ultrapure, Ge or Si bolometers, by which 
the Te02 crystals will be surrounded so as to get rid of possible events originating in 
the germanium (which could release part of their energy in the main detector) with the 
anti-coincidence technique. A powerful method would be to thermally couple these 
active shields to the main Te02 detector. The Ge or Si auxiliary bolometers would be 
attached at the main crystal providing almost complete coverage without a dedicated 
holder. In this way, a composite bolometer will be realized with multiple read-out, ca- 
pable to distinguish the origin of the event (active shields or Te02 crystal) by means of 
the comparison among pulses coming from the different elements. A degraded alpha 
coming from outside would release its energy in the shield, originating a thermal pulse 
that, seen by the shield thermistor, would be much higher than a pulse corresponding 
to the same energy released directly in the main Te02 crystal. In fig. |58]it is possible 
to appreciate the efficiency of the separation method. A scatter plot is shown, where 
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the amplitude of the pulses collected with the Te02 thermistor is plotted against the 
amplitude of the corresponding pulses collected with the shield thermistor. Two clearly 
separated curves identify the events generated in the shield with respect to the events 
generated in the main crystal. The pulses have been simulated using a multi-stage ther- 
mal model of the composite bolometer. In order to reduce the readout channels, the six 
thermistors of the shields could be connected in series or in parallel. In this way, the 
total number of channels would double. We believe that this is a tolerable complication. 
The size of the auxiliary bolometers would be 50x50x0.3 mm'^ or 60x60x0.3 mm^: 
therefore, the main crystal plus the shields would form a cube with a side only 0.6 mm 
larger than the main crystal alone, preserving the general CUORE structure and render- 
ing the assembling procedure substantially unchanged. This proposed solution has the 
potential to control the problem of the surface radioactivity relying on a technological 
improvement of the detector rather than on a better cleanness of the employed materials. 




Figure 59: Scatter plot of the pulse amplitudes from Te02 thermistor versus the pulse 
amplitudes from shield thermistor obtained with a 12 g Te02 surface sensitive detector. 

In order to test the exposed principles, several small scale prototype detectors were 
realized. In all of them the main absorber consisted of a Te02 single crystal, with masses 
ranging from 12 g (2x2x0.5 cm^) to 48 g (2x2x2 cm^). A Ge or Si single crystal with 
a 1.5x1.5 cm^ surface and 500 /im thickness was glued at a 2x2 cm^ face by means 
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of four epoxy spots (1 mm diameter and 50 /im thick). Neutron Transmutation Doped 
(NTD) Ge thermistors were glued at the main absorbers and at the auxiliary Ge or Si 
bolometers with six epoxy spots (0.5 mm diameter and 50 mm thickness). The NTD 
chip size was 3x1.5x1 mm^ in all cases. The resistance-temperature behavior of the 
chips is parameterized as usual by the formula R{T) = Rqc^'^^^'^^^^^ . The values of 
Rq and To were respectively 2.65 Ohm and 7.8 K. The coupling of the main absorber 
to the heat bath was realized by means PTFE blocks, as in the CUORICINO modules. 
The Ge or Si active layers were exposed to a particles. The source was obtained by 
a shallow implant of ^^"^Ra nuclides onto a piece of copper tape, facing the Ge or Si 
layer. ^^^Ra is an a emitter with a half life of 3.66 d in equilibrium with its a and P 
emitting daughters. The main a lines are at 5.68, 6.29, 6.78 and 8.78 MeV. Two weak 
lines sum up at about 6.06 MeV. The detectors were cooled down in a low power dilu- 
tion refrigerator located in the Cryogenic Laboratory of the Insubria University (Como, 
Italy). The detector thermistors were DC biased through a voltage supply and a 20 
GOhm room-temperature load resistance. Typical operation temperatures were T ~25 
mK, corresponding to a thermistor resistance R ~100 MOhm. The voltage across the 
thermistor was typically ~25 mV. We report here about the results obtained with a 12 
g detector with a Ge active layer. Similar results were obtained with the other devices 
(no remarkable difference was observed between Ge and Si shields). The energy ab- 
sorber was a Te02 crystal with a rectangular shape (2x2x0.5 cm^ dimension). The fast 
pulses coming from the active Ge layer had a rise time of ~6.5 ms and a decay time 
of ~40 ms, while the pulses coming from the main absorber were remarkably slower 
(rise time ~9.7 ms, decay time ~70 ms). The scatter plot obtained is shown in fig.l59l 
The two bands, corresponding to surface and bulk events, are well separated and the 4 
main a lines of the source are clearly appreciable. A powerful separation between the 
two event types can be achieved also by rise time selection. In fig. ISUt a) the rise time 
distribution is shown for pulses from the layer thermistor. The fast surface events are 
cleanly separated by the slow bulk events. A selection of the "fast" events identifies 
clearly the surface event band, as shown in fig. l60l b). The rise time analysis shows 
that the scatter plot is redundant to select surface events. It might be convenient then 
to get rid also of the main absorber thermistor, reducing the read-out channel to one as 
in the conventional case. On the other hand, a detailed analysis will be performed to 
check if pulse shape discrimination can allow to identify surface events using signals 
only from the main absorber thermistor. In this case, the Ge or Si layers would oper- 
ate as signal- shape modifiers and would not require additional thermistors, leading to a 
substantial simplification of the device. This experiment demonstrates dramatically the 
power of the technique we have developed to identify a superficial impact point in low 
temperature calorimeters. 
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Figure 60: Rise time distribution of the pulses from the Ge layer thermistor of a 12 g 
Te02 surface sensitive detector 



15.4 Development of scintillating Te02 bolometers 

Measuring at the same time the light and heat produced in the detector is an attempt 
already exploited in the quest for the background reduction |203 , 204 1. In the DBD, the 
electrons of the decay will produce at the same time phonons and light signals (in the 
region of MeV), while the alpha particles from residual radioactive decays are expected 
to release a negligible signal in light. In the case of WIMP interactions, the signal is 
produced by nuclear recoil, then without scintillation, contrary to the main background, 
given by beta and gamma radioactive decays. In the latter case however, the expected 
signals are in the region of 10 KeV. Light signals in bolometric detectors can be detected 
by a second bolometer, thermally decoupled form the main crystal. An R&D program is 
developed in the framework of the CUORE collaboration to produce such a hybrid de- 
tector based on Te02. A scintillation signal in a Te02 bolometer operated at 20 mK was 
previously reported [205| and has been confirmed by our preliminary measurements. 
The measured light yield reported was quite low (0.67% that of BGO) but sufficient to 
distinguish alpha particles and gammas at energies in the MeV range. The scintillation 
in pure Te02 crystals may be due to exciton luminescence or to the presence of a defect, 
for example V-0 (oxygen vacancy with one trapped electron). Anyway, whatever the 
nature of this scintillation could be, it is clear that a larger light yield will be welcome 
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for an improved discrimination power on particles. This can be done only by doping 
the Te02 crystal with an activator. The crystal growth of doped Te02 crystals is not 
easy [206|. This is why the choice of the activator ion to enter in the Te02 lattice (sub- 
stitutional position is preferred) has to be carefully studied keeping in mind two main 
criteria: luminescence properties of the free ion and host lattice possibility to accept the 
dopant. In our case, a supplementary constraint is imposed by the fact that the dopant 
for obvious reasons, has to have stable isotopes only. On the other hand, it is to be 
noted that in this application practically no restriction is to be imposed on the emission 
spectrum and/or on the decay time. With respect to the host lattice possibility to accept 
the dopant, several parameters have to be taken into consideration: ionic radius, oxi- 
dation number, coordination number, type of outer atomic shell, electronegativity and 
chemical hardness of the ion. A first campaign of crystal growth has already started 
in collaboration with the Research Institute for Solid State Physics and Optics of the 
Hungarian Academy of Sciences using four different dopant activators. We expect to 
get preliminary results at the end of 2004, in order to select the most promising dopant 
for further attempts. 

According to option (a) of section fT2!2l the R&D program about reduction of surface 
radioactivity must be completed at most within two years from the start of CUORE, in 
order to have time to transfer the new detector concepts to the CUORE mass production. 
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16 Cleaning of CUORE structure: crystals, frames and 
ancillaries 

A peculiar property of thermal detectors is that they are active over the entire volume 
and therefore strongly subject to radioactive surface contaminations. This contribution 
dominates the CUORICINO PP^Ou) background (sec. 16.31) and could be a dominant 
contribution also in CUORE (sec.fT4l) The material surface cleaning problem is therefore 
crucial for CUORE success. Notable improvements in the background reduction were in 
fact obtained by applying standard surface cleaning techniques to the detector elements 
(e.g. Te02 crystal lapping) and to the detector mounting structure (copper box and 
frames, teflon stand-ups, contact pins, etc.). In a sense, MiDBD and CUORICINO 
represent a successfuU R&D for CUORE material surface cleaning. In order to reach 
a background level lower than 0.01 c/(keV kg y) however, further improvements are 
required. In particular, all possible surface cleaning techniques must be considered in 
order to choose the most appropriate. Moreover, since cleaning techniques are deeply 
selective, effective contamination analysis methods must also be considered/developed. 
A series of up-to-date cleaning methods are summarized in the following. Each of them 
shows positive and negative features that must be considered and compared in order 
to reach a final decision before the CUORE detector mounting starts. As discussed 
in sections 16.31 and [l4l a factor of ten reduction in the copper surface contaminations 
would be enough to reach the 0.01 counts/(kev-kg-y) /5/3(Oz/) region background level 
in CUORE. The R&D program about surface cleaning (which must include also the 
development of powerful diagnostic methods) must be completed at most within two 
years from the start of CUORE (sec. 112.6b . 

16.1 Surface cleaning techniques 

A surface contaminant is defined as any material deposited on a substrate that interferes 
with any process involving the surface layer up to a certain depth or that affects surface 
properties or surface stability in an undesirable way. Surface contaminants can cover 
the whole surface (e.g. oxide layers or adsorbed hydrocarbon layers), can be limited 
to restricted areas (e.g. particles or fingerprints) or can be dispersed into an epidermal 
layer for a certain depth (e.g. inclusions of foreign material or diffused impurities). 
Recontamination (for example by adsorbed vapors from the atmosphere) means con- 
tamination picked up after the cleaning process in the external processing environment 
and before any further processing begins. 

The variability of the contamination is what mainly affects cleaning reproducibil- 
ity. Therefore, to fail the control contamination means not to control cleaning, and it 
will be expected that this will be the main problem for the cleaning operation of the 
whole CUORE structure. Cleaning is the reduction of surface contamination to an ac- 
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ceptable level in order to have no significant amounts of undesirable material. But what 
are the CUORE requirements to get an adequately clean surface? In the absence of a 
clear answer to this question, the strategy suggested by the proponents is to transfer to 
CUORE the know-how developed at INFN Legnaro National Laboratory (LNL) in the 
topic of surface treatments for Superconducting Cavities either made of bulk Niobium 
or Niobium Sputtered Copper. In addition to that, Padua University and the LNL of 
INFN are currently holding a Master in "Surface Treatments to innovative mechanical 
technologies for the Industry". This experience involves both industries and the mate- 
rial science section of Padua University, and represents the best humus for facing and 
solving "CUORE problems". 

Our goal is twofold: we would like to clean the CUORE copper structure and to 
process all the Te02 Crystals. In the following, we will list all the possible cleaning 
processes that could be significant for our purpose. Then, on the basis of the exposed 
advantages and disadvantages, we will select among them the best processes for clean- 
ing the CUORE structure. In this way we will propose two different surface treatments, 
the former for the copper structure, the latter for the Te02 crystals. 

The scope of reviewing all the possible treatments advisable for CUORE is however 
not only limited to the selection of the best cleaning steps. It is also a warning in getting 
deep in any issue where cleaning is concerned, in order to have the highest possible 
control of contamination. 

Not knowing if the surface is contaminated, indeed, does not mean that the surface 
is clean. Simple operations, normally considered of little importance can be detrimen- 
tal. For instance, it would be absolutely normal to degrease pieces, before cleaning, 
by means of Alkaline cleaners. However it is absolutely mandatory to know that al- 
kaline cleaners generally contain dissolved silicates, carbonates, borates, citrates and 
even molybdate. Residuals of these particles can remain attached in micro- valleys of a 
microrough surface. 

Abrasive cleaning 

The removal of gross contamination by abrasive cleaning includes the use of: 

• Abrasive surfaces: sandpaper, emery paper, steel wool, scotch-brite, scouring 
pads, etc. 

• Abrasive powders in a paste or fluid carrier: SiC, AI2O3, diamond, precipitated 
calcium carbonate (CaCOs), talc, kaolin, wheat starch. Cerium Oxide (widely 
used in glass lenses industrial cleaning was avoided for CUORICINO/CUORE 
Te02 crystals because of its high level of radioactive impurities). 

• Impacting particles entrained in a high velocity gas or liquid stream: vapor hon- 
ing, glass bead blasting, liquid honing, grit blasting, sand blasting. 



16 CLEANING OFCUORE STRUCTURE: CRYSTALS, FRAMES AND ANCILLARIES161 



• Abrasives combined with an etchant to provide chemical-mechanical abrasion and 
polishing. 

Abrasive particles can be used wet or dry and with various particle sizes. Particles 

can be accelerated in a high velocity water stream or in a gas stream by using a pressur- 
ized siphon system. In addition to removing gross contamination, the following results 
can be obtained: 

• Abrasive embedding into the surface even up to hundreds of microns 

• Roughening of the surface 

• Introduction of microcracks into the surface of brittle materials 

• Production of unacceptable distortion of metallic surface due to the compressive 
stresses introduced by a king of "shot peening" 

Blasting done with plastic media is less damaging than with harder grits; however 

it can leave organic contamination on the surface). There is however a clean alternative 
to all that: the cleanest mechanical polishing is ice micro-droplets blasting in a high 
pressure nitrogen stream. This system has been conceived for the cleaning of super- 
conducting cavities, but it is only in the prototype stage and any application to CUORE 
should start ab initio. 

Chemical etching 

Chemical etching removes surface material along with strongly adherent contami- 
nants. Advisable procedures before chemical etching are: 

• Alkaline pre-cleaning, in order to get uniform wetting and etching 

• Pickling, i.e. the removal of the large amounts of oxides formed during fabrication 
process. 

• Bright dip, which is nothing else than a Mild Pickling 

It must be kept into account that acid cleaning of metals can have the detrimental 
effect of introducing hydrogen into the surface. If hydrogen embrittlement is a concern, 
the part should be high-temperature vacuum baked after etching. Sometimes chemical 
etching does not remove some constituents from a surface and leaves a smut that must 
be removed by further etching. Whenever other elements in addition to Copper frames 
should be treated, a chemical etching using Etchant Vapor instead of a Liquid Phase 
could be also considered. 
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Figure 61: Zoom of the micro-roughness of a metal surface immersed in an electrolytic 
bath 

Electrochemical polishing 

In an electrolysis cell, the surface of an electrical conductor can be removed by mak- 
ing it the anode of an electrolytic cell and "de-plating" the material. This is electro- 
etching or electrocleaning or electrodissolution. The higher the current density, the 
more roughened the surface (fig. I^TT) . Oxygen is released at the anode: it will react 
oxidizing contaminants on the surface at low anodic potentials. 

Electropolishing is a much smoother process, where the working current density is 
carefully chosen. In case of Cu electropolished in a H3PO4 solution, in a first approx- 
imation, the leveling action is due to the protection of the smooth areas by a deposited 
phosphate material, and to the erosion of the exposed peaks. 

Electropolishing leaves a phosphate film on the surface, which can be removed either 
by leaving pieces in the solution without any applied voltage or by an HCl rinse. 

A typical I-V characteristic for the electropolishing of Copper in Ortophosphoric 
acid solution in the case of planar and parallel face electrodes and with negligible edge 
effects is shown in fig. |62l The following behaviour is observed at the voltages corre- 
sponding to different sections of the curve: 

• Over section Va and Vb the current increases as a linear function of the voltage. 
Material dissolution happens with a too low rate. The process is accompanied by 



] 6 CLEANING OF CUORE STRUCTURE: CRYSTALS, FRAMES AND ANCILLARIES\69 




Figure 62: Typical I-V characteristic for electropolishing 



the evolution of Oxygen bubbles that promote local pittingof the surface of the 
piece by sticking to the anode. 

• Polishing effect is observed between Vb and Vc; the grain structure is brought 
into relief, as long as the process takes place. However roughness levelling and 
brilliant surface are obtained at voltages close to Vc. Here at the end of the plateau 
there is the minimum of the evolution of Oxygen bubbles. Even a minimum 
amount of bubbles can represent a limitation to the roughness smoothing action. 
Migrating toward the top they produce undesired vertical traces depending on the 
solution agitation. 

• At a higher potential, the gas evolution becomes stronger and the surface erosion 
is accompanied by pitting. A better surface quality is obtainable at voltages well 
above Vc, since Oxygen bubbles have not enough time to stick to the surface. 

In the framework of more than 18 years experience in the field of superconducting 
cavities, the LNL of INFN have well-established experience in Copper electropolishing. 
At LNL, indeed, a method for controlling the electropolishing process has been devel- 
oped. We try to minimize parasitic problems (such as gas evolution) by searching the 
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Figure 63: 1-V characteristic curve for a standard electropolishing process (circles). The 
differential conductivity curve dl/dV (continuous line) and the ratio I/V (shaded line) are 
displayed versus voltage. We interpret the first minimum in differential conductance as 
due to edge effects and to the not uniform distance between electrodes 

best working point of the I-V characteristic and by following it during the whole pro- 
cess. The method consists in a computer assisted monitoring of the I-V characteristic. 
The process is driven by voltage. An automatic program displays the numeric derivative 
of I versus V. The working point is chosen as the minimum of such a derivative, i.e. the 
minimum of bath differential conductance that corresponds to the point of maximum 
resistance of the viscous layer. 

In figure |63l we display the differential conductance compared to the ratio I-V and 
together with the I-V characteristic. Some literature approximates the minimum of 
dl/dV with the minimum of I/V. Already from the figure below, it is visible that the 
two minima differ quite substantially in voltage. In any case that of an electrolytic cell 
is a non-linear circuit, hence the solution conductivity is a differential quantity while, 
the ratio I/V has no any physical sense in such a case. 

Our method allows the determination automatically of the best working point, when 
electropolishing. Its great advantage is that, by computer controlling, the operator does 
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not need to know what is the metal being cleaned, and what acid is being used. 

Electroless electrolytic cleaning 

This method exploits the difference in electromotive potentials to remove material 
from one surface and deposit it on another. Such displacement-type electrolytic cleaning 
processes are commonly known for Silver. When immersing Silver in an undiluted 
solution of pure household ammonia contained in an Al or Mg tray, the Silver surface is 
cleaned as the Al or Mg is oxidized. 

Deburring 

Deburring is the removal of the rough edges produced in cutting or shearing. Debur- 
ring is performed by: abrasion, chemical etching, laser vaporization and "flash debur- 
ring" (which uses a thermal pulse from an exploding gas-Oxygen mixture to heat and 
vaporize the thin metal protrusions). 

Solvent cleaning 

Solvents dissolve the contaminant taking it into solution. It is possible to use: i) 
Polar solvents, such as water and water-alcohol mixtures, used to dissolve ionic salts 

which are polar contaminates; ii) Non-polar solvents, such as chlorinated hydrocarbon 
solvents, used to remove non-polar contaminates such as oil. A mixture of solvents 
can dissolve both polar and non-polar contaminates. A solvent's effectiveness is deter- 
mined by the solubility parameter for specific contaminants (i.e. the saturation maxi- 
mum amount of a specific contaminant dissolvable in a specific amount of the solvent). 
Of course the solubility parameter will be temperature dependent. 

Solvent cleaning leaves a surface layer of residue which must be removed. This 
removal can involve: i) water rinsing; ii) an elevated temperature; iii) other solvents 
that displace the surface layer. For example, both for cristals and copper frames; a 
solvent wipe-clean cleaning sequence might be Trichloroethylene-Acetone-Methanol- 
Isopropanol. Good polar solvents for a variety of polar contaminants (such as ionic 
salts) are Water and Water- Alcohol mixtures (typically in a ratio of 1:1). The addition 
of alcohol lowers water surface energy allowing it to penetrate into "hideouts" to remove 
hidden contamination. 

Water and water-alcohol mixtures are not good solvents for non-polar contaminants 
such as oils. 

Chlorinated and chlorofluorocarbon solvents 

Chlorinated hydrocarbon solvents, such as Trichloroethylene, are preferred to hydro- 
carbon solvents for their lower flammability, but there is concern with their toxicity and 
carcinogenic properties. Chlorofluorocarbons are not carcinogenic and can be used in 
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large quantities. They are, however, less effective than chlorinates solvents in dissolving 
contaminants. 

1,1,1-trichloroethane (TCA-CHCCI3, or methyl chloroform) is a common chlori- 
nated solvent widely used in vapor degreasers. TCA has high Permissive Exposure 
Level (PEL), is classed as non- volatile, has low toxicity, but it has a high Ozone Deplet- 
ing Potential (ODP) rating. Hence, in order to meet regulation standards. Vapors must 
be contained in enclosed vapor systems. 

A possible alternative to TCA is Methyl-ene chloride (MEC-CH2CI2). It has a rather 
low PEL, but it has a boiling point of 39.8°C. This makes it more difficult to handle 
but it has the advantage of being applicable for vapor degreasing temperature-sensitive 
materials. It must be kept in mind that MEC is a very aggressive solvent; it can damage 
plastics and rubbers. 

Perchloroethylene (PCE or PERC - Cl2C:CCl2) is particularly suitable for being 
applied to Copper frames after machining. It has a high boiling point (121. 1°C), is 
useful for dissolving heavy greases, and it can contain a large amount of water without 
degrading its solvency power. 

Solvents can be mixed to give synergistic cleaning actions: a typical blend used for 
metal vapor degreasing is the Azeotrope mixture of Freon TF with methylene chloride 
or with ethanol. When using solvents however great attention must be payed to the 
performed operation. Some solvents can react with the surface being cleaned. For 
example. Chlorinated solvents can react with water to form HCl which can react with 
many metals, particularly Al, Mg, Be, Zn , but also Cu and Te to form inorganic salts. 

Liquid CO2 as a solvent 

Liquid CO2 has low surface tension and low viscosity. LCO2 at 20°C, which is below 
the Supercritical Point, is a well-known good solvent for cleaning metals 

Supercritical fluids 

If a gas, such as CO2, is compressed to its "critical pressure", it liquefies to be- 
come a "critical fluid". If it is also heated above its "critical temperature", it becomes 
a supercritical fluid. Critical fluids and supercritical fluids are good solvents for many 
medium-molecular-weight, non-polar or slightly polar organics. Solubility parameter 
for which CO2 can vary from in the gas phase to 10 under high pressure supercritical 
conditions (SCF-CO2- critical point 31°C, 74 bar pressure). Values of 6-8 are typical. 
Supercritical CO2 fluid has the following advantages: 

• it is stable 

• has low toxicity 

• minimal cost 
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• is a solvent for many organic materials 

Semi-aqueous cleaners 

Semi-aqueous clenears refer to solutions of natural or synthetic organic solvents, 
used in conjunction with water in some part of cleaning cycle. Water immiscible semi- 
aqueous cleaners include: terpenes; high-molecular- weight esters; petroleum hydrocar- 
bons; glycol esters. Terpenes are natural hydrocarbons such as: the d-limonene; the 
alpha- and beta-pinenes, derived from citrus and pine oils. They may be as effective as 
the CFCs in many instances, but have a greater tendency to leave residues and suffer 
from the fact that they are slow drying. 

Water miscible semi-aqueous cleaners are low-molecular- weight alcohols; ketones; 
esters; and organic amines. Acetone (CH3COCH3) removes heavy oils quite effectively 
but tends to leave a residue, that should be removed by a methanol rinse. 

Saponifiers, soaps, and detergents 

Alkaline saponifiers (generally silicate and phosphate-based) convert organic fats to 
water soluble soaps. Alkaline cleaners are generally used hot. Mild alkaline cleaners 

(pH of 8-10) often have dissolved silicates, carbonates, borates, and citrates and should 
be used to clean alkaline-sensitive materials, such as Al and Mg. Strong alkaline clean- 
ers, caustics, (pH = 12) may have water, sodium silicate, sodium molybdate and sodium 
fluoroborate. 

Sodium silicate developes electrostatic forces that displace the contaminants while 
depositing a glassy film that prevents recontamination. The glassy material is generally 
removed in the DI water rinse. After using alkaline cleaners, the surface should be 
followed by an acid dip prior to the water rinse to remove alkali salts since alkali salts 
adhere strongly to surfaces and are difficult to remove by "water rinsing". Clean oxide 
surfaces, in particular Te02 monocrystals strongly adsorb hydrocarbons, and detergents 
or solvents normally do not completely remove the hydrocarbons; alkaline or oxidative 
cleaners must be used to remove the remaining hydrocarbons. 

Detergent cleaning is a comparatively mild cleaning technique. In detergent clean- 
ing, the detergent surrounds contaminants, emulsifying them into suspension, without 
actually dissolving the material. This emulsifying action is assisted by wetting agents 
and surfactants which loosen the contaminants from the surface. The most common 
detergents are soaps which are the water-soluble reaction product of a fatty acid ester 
and an alkali (usually sodium hydroxide). In solutions, pH adjusters are used to aid in 
the cleaning action. Generally it is found that basic solutions clean better than acidic 
solutions if chemical etching is not involved. The pH of the cleaning solution is often 
adjusted to be basic, using ammonia or ammonium hydroxide. 

Chelating agents (sequestering agents) keep into solution the normally in-soluble 
phosphates. Glass cleaning solutions use chelating agents such as ethylene diamine 
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tetraacetic acid (EDTA) and citric acid with salts containing hydroxyl and amine substi- 
tutes. 

Ultrasonic cleaning 

Low frequency ultrasonic cleaning relies on the jetting action of collapsing cavitation 
bubbles in contact with a surface to provide a high pressure jet of fluid against the 
surface. Ultrasonic cleaning is often a good way to remove loosely adhering particles 
after a grinding or abrasive procedure and can be used with solvents to remove adsorbed 
contaminants. Ultrasonic jetting is good for removal of large particles but less efficient 
as the particle size decreases into the submicron range. The cavitation bubbles are 
formed by the tension portion of an ultrasonic wave in a fluid media and grow with 
time. The size that can be attained depends inversely on the frequency and the surface 
tension of the fluid. High frequencies (>60 kHz) give smaller bubbles and a higher 
bubble density. The ultrasonic wave is produced by magnetostrictive or electrostrictive 
transducers which can be attached to the fluid-containing tank walls or immersed in 
the fluid in the form of a probe that can concentrate the ultrasonic energy into a small 
area. Typically the transducers operate at 18-120 kHz, at an energy density of about 100 
watts/gal of fluid. The ultrasonic cleaner size can be from 5 gallons for a small cleaner, 
up to very large systems using many transducers. The size of cavitation bubbles in 
the fluid depends on the vapor pressure, surface energy, and temperature of the fluid. 
For example, pure water at 60°C and 40 kHz has a maximum cavitation bubble size of 
about 100 microns if a surfactant is present, the bubble size is smaller due to the lowered 
surface energy. The jet pressure from the collapsing bubble can be as high as 300 psi. 
The cavitation jetting is more energetic for cooler media and when there are no gases in 
the bubble to hinder its collapse. 

Ultrasonic cleaning must be used with care on Te02 crystals, since the jetting action 
can produce high pressures that cause erosion and introduce fractures in the surface. 

Megasonic cleaning 

High frequency (>400 kHz) ultrasonic cleaning does not cause cavitation. Instead, the 
action consists of a train of wave fronts that sweep across a smooth surface producing 
disruption of the viscous boundary layers on the substrate surface by viscous drag. The 
resulting pressure is less than 50 psi and does not hurt fragile surfaces. High frequency 
transducers can be focused to restrict the area of impact and allow lateral fluid flow from 
the area of concentration. 

Megasonic cleaning utilizes high frequency (850-900 kHz) transducers to produce 
non-cavitating pressure waves. The megasonic agitation system is applicable to smooth 
surfaces, particularly for removing particles, but does not work well on configured sur- 
faces where the surface is shadowed from the pressure wave. 
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Wipe-clean 

In some cases the surface can not be immersed in a fluid and must be cleaned by 
wiping. Wiping with a fluid should be done with a moist lint- free cloth or sponge which 
has no extractables when in contact with the wiping-fluid. The wiping motion should be 
a rolling motion such that contamination that is picked-up does not come into contact 
with the surface as wiping proceeds. 

Strippable coatings 

Particles can be removed from surfaces by covering the surface with a liquid poly- 
mer, allowing it to solidify, then mechanically stripping (peeling) the polymer from the 
surface. This technique is used by the optics industry to remove particles from mirror 
surfaces and protect surfaces from abrasion during assembly, and in silicon technology 
to remove particulates from silicon wafers. There are many types of "strip coats," each 
coating leaves different residues on stripping and have differing corrosion compatibility 
with surfaces. Hydrocarbon residues left by strippable coatings can be removed by ox- 
idation techniques. A strippable coating is not necessarily organic, but it can be also a 
sputtered getter film, for example a Niobium, Hafnium or Titanium film sputtered onto 
copper for Oxygen purification. The film is then removed by chemical etching. 

Rinsing 

After any wet cleaning process the surface should be thoroughly rinsed in pure or 
ultrapure liquid, usually water, before allowing to dry. This avoids leaving residues 
on the surface. A common rinsing technique is to use successive rinses (cascading 
rinsing) in pure or ultrapure water until the rinse water retains a high resistivity (e.g., 
>12 megohm). A problem can be the "dragout" of one fluid with the part which then 
contaminates the subsequent fluid tank. For the beginning rinse, a sheeting agent can 
be added that lowers the surface tension of the water and aids in flowing the rinse water 
off the surface. After rinsing, the surface should be dried as quickly as possible since 
the residual water film on the surface will cause particles to stick to the surface and on 
drying the particles will adhere very tenaciously. 

Drying and outgassing 

Displacement drying uses anhydrous fluids, such as isopropanol, anhydrous ethyl 
alcohol denatured with acetone or methanol, or a commercial drying agent. One of the 
best drying techniques is a "vapor dry" where the cold surface is immersed in the vapor 
above a heated anhydrous alcohol sump. The cold surface condenses the alcohol vapor 
which flows off into the sump taking water and particulates with it. When the surface 
becomes hot, condensation ceases and the hot surface, when withdrawn dries rapidly. 
Surfaces can be blown dry using a low velocity dry gas stream. When blowing, a nozzle 
with a 0.2 micron or smaller particulate filter should be used in the nozzle. In addition. 
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when drying insulator surfaces, the gas should be ionized to prevent charge build-up 
on the surface. The gas can be ionized with an electronic (corona), or a laser ionizer. 
But be aware that Electronic ionizers can arc and produce particulates. A high velocity 
jet of gas can be shaped to blow-off a moving surface. The jet should impact the on- 
coming wet surface at about a 30° angle. Hot gas drying or evaporative drying, uses the 
recirculation of hot dry, filtered air over the surface to promote evaporation. This drying 
technique has the problem of "water spotting" if the fluid is not ultrapure. 

Volatile materials from the bulk of the material are removed by outgassing. Obvi- 
ously Te02 monocrystals will not outgas appreciably, while Outgassing is especially 
important for Copper. Since diffusion is required, the time to outgas a material may be 
very lengthy if the diffusion rate is slow and/or the diffusion distance is long. Gener- 
ally metals primarily outgas hydrogen, particularly that taken up during acid cleaning 
or electropolishing. 

After drying, the cleaned parts should be stored and transported in a manner that 
does not unduly recontaminate the parts. 

Reactive cleaning 

Reactive cleaning uses liquids, gases, vapors, or plasmas to react with the contami- 
nant to form a volatile or soluble reaction product. Reactive cleaning liquids are often 
oxidizing solutions. Many acid-based systems can be used as oxidants. Nitric acid can 
also be used as the oxidizing agent. Nitric acid with hydrofluoric acid is used to oxi- 
dize/etch surfaces. Hydrogen peroxide (H2O2) is a good oxidizing solution for cleaning 
glass. [64] Often boiling 30 unstabilized H2O2 is used. Hydrogen peroxide is often stabi- 
lized, which reduces the release of free oxygen. Unstabilized H2O2 must be stored in a 
refrigerator to slow decomposition. Hydrogen peroxide is sometimes used with ammo- 
nium "hydroxide", to increase the complexing of surface contaminants. However the 
decomposition rate of the H2O2 is greatly increased by combination with ammonium 
hydroxide. 

Oxidative cleaning can be performed using chlorine-containing chemicals. For ex- 
ample, a water slurry of sodium dichloroisocyanurate (i.e. swimming pool chlorine) 
which has 63% available chlorine, can be used to scrub an oxide surface to remove 
hydrocarbon contamination. This combines mechanical scrubbing with oxidation and 
improves the cleaning action. 

Anodic oxidation in an electrolysis cell can be used to clean surfaces. The contam- 
inated layer can be removed by anodically oxidizing the surface in an electrolytic cell. 
The oxide layer must be subsequently removed. 

Gaseous oxidation cleaning can be used on surfaces where surface oxidation is not a 
problem. Oxidation is usually accomplished using oxygen, chlorine, fluorine, ozone, or 
Nitric oxide which creates volatile reaction products such as CO and CO2. Reactive gas 
cleaning may use a reaction with a gas at high temperature to form a volatile material. 
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High temperature air fire is an excellent way to clean surfaces that are not degraded by 
high temperature. 

Ozone cleaning 

The use of oxidation by Ozone created by ultraviolet radiation (UV/Ozone clean- 
ing) at atmospheric pressure and low temperature has greatly simplified the production, 
storage and maintenance of hydro-carbon-free surfaces. Historically the UV/0, clean- 
ing process was developed because of the need to clean very delicate quartz oscillators. 
Any attempt to clean the quartz using physical contact caused breakage. The UV/03, 
cleaning technique provided a non-contacting way to clean the delicate quartz plates. 
The short wavelength radiation causes bond scission in the hydrocarbon contaminants 
and generates ozone which reacts with carbon to form volatile CO and CO2. Typical 
exposure times for UV/03 cleaning are from a few minutes to remove a few monolayers 
of hydrocarbon contamination to hours or days or weeks for storage of cleaned sur- 
faces. The UV/03 cleaning technique has the advantage that it can be used as a dry, 
in-line cleaning technique at atmospheric pressure. The UV/03 cleaning technique is 
also useful for cleaning holes (vias) in surfaces. 

Hydrogen cleaning 

High temperature hydrogen or forming gas (90% N2: 10% H2), can be used to remove 
hydrocarbon contamination from a surface by hydrogenating the material and making 
it more volatile. Hydrogen reduction of oxide layers can be used to clean surfaces 
in a furnace environment. However, Hydrogen cleaning can also change the surface 
chemistry. For example, hydrogen firing of Te02 could produce metallic Tellurium 
surface by reducing the Oxyde to metal on the surface. 

Reactive plasma cleaning and etching 

Reactive plasma cleaning uses a reactive species in the plasma to react with the sur- 
face to form a volatile species which leaves the surface at much lower temperatures than 
those used for reactive gas cleaning. The additional requirement on reactive plasma 
cleaning is that it does not leave a residue. Oxygen or hydrogen are often used for 
plasma cleaning while fluorine (from SFq, CF4 , CHF3, C2F6 CaFg, or SFg) and chlorine 
(from HCl, CCI4, or BCI3) are used for plasma etching. The reactive plasma clean- 
ing/etching technique is typically specific and can be used to selectively remove the 
oxide from the surface and then have a low etch rate for the substrate material. 

Oxygen (or air) plasmas are very effective in removing hydro-carbons and absorbed 
water vapor from surfaces. A typical plasma cleaner has plasma generated by an rf 
discharge and the surfaces to be cleaned are in a "remote" or "downstream" location and 
not in the plasma generation region. The surface attains a potential (sheath potential) 
that is negative with respect to the plasma, and ions are accelerated from the plasma 
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to the surface. For the case of a "cold plasma" which has low energy particles, this 
sheath potential will only be a few volts. When the plasma particles are more energetic 
or the electrons are accelerated to the surface, the sheath potential can be tens of volts. 
In addition to being bombarded by ions, the surface in contact with the plasma will be 
bombarded by "activated species", excited species, thermal species, and high energy 
photons (UV and, under some conditions, soft X-rays). Ions and excited species will 
release their energies of ionization or excitation when they impinge on the surface. 

Plasma etchers and strippers typically use more aggressive reactant gases such as 
chlorine or fluorine and are constructed to withstand corrosion and pump the particulates 
that are often formed in the etching and stripping process. 

When plasma is etching copper, a copper chloride (CUCI2) residue is left on the 
surface which can be volatilized by heating to above 200 C. Often mixtures of gases are 
used for etching and cleaning. Oxygen is often added to the fluorine plasma to promote 
the formation of atomic fluorine and to oxidize the surface and thus increase the etch 
rate. One of the most common gas mixtures used to etch is 96 CF4. Helium is often 
added to dilute the mixture and to increase the thermal conductivity of the plasma thus 
reducing the temperature rise of the substrate during etching. 

Plasma cleaning 

Plasmas can be used to clean surfaces in the deposition system in the same manner as 
they are used to condition vacuum surfaces. In some PVD deposition systems which are 
not normally used with a plasma, a "glow bar" or "glow plate" is used as the cathodic 
or anodic electrode of a DC discharge to create the plasma. The larger the area of the 
surface, the better the plasma distribution in the system. Plasma cleaning can be done 
using an inert gas plasma or can use a plasma containing a reactive gaseous species 
to form a volatile reaction product from the interaction of the gaseous species and the 
surface species. 

"Ion scrubbing" of a surface occurs when a surface, which is in contact with an 
inert gas plasma, develops a wall sheath and is bombarded by inert gas ions accelerated 
across this wall sheath. Generally the ion energy is too low to cause surface damage or 
physical sputtering, but is effective in the desorption of adsorbed surface contaminants 
such as water. 

Reactive plasma cleaning/etching can be done in the deposition system in much the 
same way as was described in the "external" plasma cleaning. The surface in contact 
with a plasma containing reactive species develops a negative potential with respect to 
the plasma (self -bias). Ions, along with neutrals and "activated" species, of the reactive 
species bombard the surface producing volatile reaction product either with contam- 
inants (cleaning) or the substrate material (etching). The most common reactive gas 
used is oxygen or air. 

A major concern in any plasma process is to obtain a uniform plasma over the sur- 
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face. Some plasma-generation configurations are more amenable to uniformity than 
are others. The magnetron configuration is one where plasma uniformity is difficult to 
obtain except in certain applications such as passing the substrate through the plasma 
of a planar magnetron. For CUORE magnetron configurations using moving electro- 
magnetic fields could be conceived in order to obtain uniform etching without substrate 
damage. 

The reactive etching/cleaning processes produce volatile species which may be de- 
posited in other parts of the system where there are different plasma conditions. This 
may have a detrimental effect on the gas handling/ pumping system and can be a source 
of particulates in the etching system, that must be periodically reconditioned. 

Sputter cleaning 

Sputter cleaning uses physical sputtering, not chemical reaction, to remove some 
of the surface layer, which includes the contaminates. Sputter cleaning is a kind of 
universal etch since conceptually everything can be removed by the sputtering process. 
However certain types of surface contamination, such as particles and inclusions of 
inorganics compounds, are very difficult to remove by sputtering because of their shape. 
This cleaning process can be easily integrated into the deposition process so as to allow 
no time for recontamination between the cleaning and the deposition process as is done 
in ion plating. Sputter cleaning however could produce detrimental surface damage on 
Te02 crystals, unless very low power and low density plasmas are ignited. For Te02, 
indeed Helium plasmas are foreseen. 

During sputter cleaning, the bombarding gas may become incorporated into the sur- 
face and subsequently released on heating. Sputtering from a plasma environment has 
the disadvantage that gaseous contamination in the plasma becomes activated and can 
react with the surface being cleaned; also sputtered species can be returned to the surface 
by scattering (redeposition) and contaminate surface species can be recoil-implanted 
into the surface. Sometimes this makes sputter cleaning difficult. Ion milling, where 
ion beam sputtering is used to remove surface material, can be done in a vacuum en- 
vironment where the sputtered species are not redeposited on the substrate surface and 
gaseous contamination is rapidly pumped away. 

Ion beam cleaning 

The use of ion and neutral beams allow the cleaning/etching of a surface in a good 
vacuum environment. Energetic ion beams of reactive species can be used to clean/etch 
surfaces and the process is called Reactive Ion Beam Etching. Beams of uncharged 
radicals of reactive species (H, CI, F) can be used to clean surfaces in vacuum. The use 
of energetic inert gas ion beams to bombard a surface concurrently with a molecular 
beam of the etchant gas (Ion Beam Assisted Etching-IBAE) shows enhanced etching 
over either the inert ion bombardment (sputtering) or the molecular beam alone. 
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16.2 Conclusions 

From the analysis of the possibiUties Usted above, we learn that two different cleaning 
procedures must be conceived for Copper and for Crystals. 

For cleaning the Copper structure delivered after fabrication, indeed, we propose the 
following procedure: 

• Mechanical finishing 

• Degreasing 

• Tumbling or vibratory finishing 

• Rinsing 

• Bright dipping 

• Rinsing 

• Electropolishing 

• Rinsing 

• Ultrasound rinsing 

• Chemical etching 

• Copper passivation 

• Rinsing 

• Ultrasound rinsing 

• Drying 

• Vacuum degassing 

• Sputter cleaning 

• Storing in controlled atmosphere environment 

On the other side for cleaning Te02 crystals delivered after fabrication, indeed, we 
propose the following procedure: 

• Alcohol Rinsing 

• Drying 
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• Helium gas low pressure and low power Sputter cleaning 

• Ion beam cleaning 

• Storing in controlled atmosphere environment 

Tests of both procedures will be carried out in the next feature to assess their effec- 
tiveness. 
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17 Enrichment option 

CUORE is the only proposed next generation /3/3(Oz/) decay experiment that plans the 
use of non-enriched material. Indeed, the high isotopic abundance of ^^°Te allows very 
high sensitivities even using the much cheaper natural abundant material. This ob- 
viously implies that the construction of a larger array would not only be technically 
possible but also fiscally acceptable. On the other hand the high natural abundance 
of ^•^°Te allows a much simpler and less expensive enrichment procedure than that for 
other double beta decay candidates. Moreover Te02 crystals made with ^^°Te isotopi- 
cally enriched material have been already tested during the MiDBD experiment with 
results consistent with those obtained for natural crystals, both from the point of view 
of detector performance and detector background. As discussed in section lOl a 95% 
enrichment in ^•^°Te of CUORE with a background of 0.001 counts/(keV-kg-y) would 
imply a 5 year sensitivity F/j of 1 .9 x 10^'' y implying a sensitivity on the Majorana mass 
parameter in the range 8-45 meV. 

We are presently considering a "partial enrichment" option for CUORE which would 
leave unaltered the CUORE structure shown in Fig|3T] Four central columns of CUORE 
would contain 9 inner planes made by 95% enriched Te02 crystals with the upper and 
lower two planes made by natural ones. We would have therefore in total 144 enriched 
and 844 natural crystals, respectively. The sensitivity of CUORE on the Majorana neu- 
trino mass would be almost doubled if the background of the enriched crystal would 
be the same as that of the natural ones. In addition comparison between the spectra 
from the enriched and natural crystals would be precious to understand the origin of the 
background especially if a peak indicative of neutrinoless double beta decay would be 
found. Our plan would be to produce in the USA and test at LNGS three 5x5x5 cm^ 
crystals enriched in ^^°Te plus three natural ones, within two year from the beginning 
of the project to ensure that the thermal and background properties would be the same 
as those of CUORICINO. If this would not be the case the option of CUORE entirely 
made by natural crystals will be adopted. 
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18 Structure of the collaboration 

The CUORE collaboration involves physicists from 14 active groups in Europe and 
United States, plus about 50 technicians and a number of undergraduate students. Five 
Laboratories of international relevance participate to it: the Lawrence Berkeley National 
Laboratory, the Lawrence Livermore National Laboratory, the Laboratori Nazionali del 
Gran Sasso, the Laboratori Nazionali di Legnaro and the Kamerling Omnes Cryogenic 
Laboratory in Leiden. All the specific services of these Laboratories will participate to 
the various technical problems involved by CUORE: cryogenics, chemical, construc- 
tion, material properties, electronics, low radioactivity measurement etc. In addition a 
collaboration has been initiated with the Department of Structural Engineering of the 
Politecnico of Milan, one of the leaders in Europe, for the design and construction sim- 
ulation of the CUORE structure. With respect to the other experiments on the same sub- 
ject, a peculiarity of CUORE is the multidisciplinary aspect required by its construction 
and operation. For this reason we would like to mention the "historical" background of 
the various participating group and their familiarity and internationally acknowledged 
skill in the specific fields of the project. We will also mention how they will collaborate 
on the various tasks. 

1. The Berkeley group, in addition to its well known experience in low energy nu- 
clear physics, has an outstanding skill in the fabrication and operation of NTD 
(Neutron Trasmutation Doped) thermistors obtained by controlled irradiation in 
two USA reactors. Since 15 years they have established a close collaboration with 
the Milano group in the development of thermal detectors. 

2. The Milan, South Carolina and Zaragoza groups have largely contributed in the 
field of double beta decay experiment since more than 30 years. The South Car- 
olina group is concentrating on the issue of other candidate PP{Ou) bolometer 
materials. 

3. The Leiden group is one of the leaders in the world arena in the construction and 
operation of dilution refrigerators and more in general in the field of Cryogenics 

4. the Berkeley, Livermore, Milan (with Como) , South Carolina and Zaragoza 
groups are working since tens of years and often in coUaborarion in the study 
and reduction of radioactive background 

5. the Legnaro group is a leading team in Europe in the field of surface treatment 
and cleaning 

6. the Florence group has since tens of years a deep knowledge in the material prop- 
erties at low and very low temperature, and is presently collaborating with Como 
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and Milano in the characterization of a various samples of NTD thermistors pro- 
duced in Berkeley 

7. The Genoa group was and is involved since many years in the field of data acqui- 
sition 

8. The Milan group is the only one in the world performing searches on double beta 
decay with the cryogenic technique since fifteen years 

All these groups are since the last year collaborating to the development and operation 
of CUORICINO, which is a first step towards the realization of CUORE. 
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19 Cost statement 

In this section we give an estimate of the costs of the basic CUORE design described in 
this proposal with possible sources of these funds. 

Table 22: Breakdown of the total CUORE cost estimates (kEuro) 
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150 


100 


100 


670 




670 


Cleaning 
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50 
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700 
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Shieldings 
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600 




600 
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50 


300 


500 
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850 


Electronics 
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40 


40 


480 




30 


590 


590 




Data Acquisition 


INFN 


10 


30 


155 


155 


170 


520 




520 


CUORE Installation 


INFN/US 






50 


50 


900 


1000 


500 


500 


Contingency 


INFN/US 


300 


700 


700 


600 


600 


2900 


1450 


1450 


Total cost 




920 


2670 


5635 


2855 


1850 


13930 


7440 


6490 



Te02 crystals 

The final cost of the 988 5x5x5 cm'^ Te02 crystals, as agreed some time ago with 
SICCAS if of the order of 4,300 Euro per crystal (comprehensive of Te metal, chemical 
reduction to Te02 , crystal growth and preliminary lapping). This would amount to a 
total of about 4.3 million Euro. 

Shields 

According to the current quotation of normal and low activity lead and the basic 
design described in section lOl the total cost would amount to 500,000 Euro. An addi- 
tional cost of 100,000 Euro would be required for the external neutron shield. 

Electronics 

As described in section [TOl the front-end electronics for CUORICINO was designed 
taking into account CUORE requirements. Few modifications will be therefore required 
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and the cost of each channel (from thermistor bias circuit to DAQ input) can be precisely 
estimated: 590 Euro. This cost would include also a 10% excess of spare channels. 

Dilution refrigerator and cryogenic equipment 

This section would include the dilution refrigerator system, the CUORE cryostat and 
the liquid helium and nitrogen systems. A rough estimate amounts to one million Euro. 

Underground laboratory 

As described in section 13 the underground building will house (besides the CUORE 
cryogenics, shields, electronics and DAQ) a clean room for detector assembling and a 
Faraday cage for front-end electronics. The total cost of these structures would amount 
to 850,000 Euro (200,000 for the Faraday cage, 400,000 for the clean room, 150,000 for 
the laboratory and 100,000 for the UPS's). 

Detector structure and machining 

The cost of the detector structure materials will be dominated by the required amount 
of copper (5 tons) whose cost can be estimated around 50,000 Euro. PTFE and other 
materials would cost 10,000 Euro. The main conribution would come however from 
detector structure machining whose cost is estimated of the order of 500,000 Euro. 

DAQ system 

At present, only a very rough estimate of the cost of the DAQ system is possible. 
In order to satisfy CUORE requirements in fact, a new design is required with respect 
to CUORICINO (secinil- A rough estimate of the total cost amounts to 400 Euro/ch. 
Additional 50,000 Euro will be required for the slow control system. 

Surface cleaning facilities and procedures 

As previously discussed (secfT6t. surface cleaning for CUORE requires an inten- 
sive R&D phase. A final cost for the whole item will be possible only after we'll get 
a final decision. A rough estimate is however possible by extrapolating our CUORI- 
CINO experience. In particular, excluding the cost of cleaning materials (which could 
be extremely variable in dependence of the purity requirements) we can estimate a total 
budget of 700,000 Euro, including both R&D (mainly at INFN LNL) during the first 
1-2 years and facilities/procedures costs (at LNGS). 

NTD thermistors 

In order to produce 2000 NTD Ge thermistors for CUORE, we estimate a cost of 
approximately 600,000 Eu which will be used to pay for the purchase of the Ge ma- 
terial, the reactor irradiations, the determination of the doping levels achieved, and the 
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processing of the NTD material into thermistors. We also estimate the time required to 
do this production to be on the order of three years. 
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20 Conclusions 

It should be abundantly clear from the foregoing discussion that the most important 
remaining issue in neutrino physics is the determination of the neutrino mass scale. 
Neutrino oscillation experiments have clearly demonstrated that: 

1. neutrinos have mass and that their mass eigenstates mix; 

2. the mixing matrix elements have values that imply that next generation /9/9(0z/) de- 
cay experiments should be able to determine the scale of neutrino mass in the case 
of Majorana neutrinos. 

The main advantage of CUORE is that it will utilize Te02 crystals of natural isotopic 
abundance Te even if an enriched option could also be considered at an affordable cost. 
CUORE will offer also the unique possibility to search for /5/5(Oi/) of isotopes other 
than Te. On the other hand, the bolometric technique has yielded energy resolutions 
approaching those of intrinsic Ge detectors and energy resolution is crucially important 
for an actual observation of PP^Ou) and a measurement of \{m^)\ . Important results 
on background reduction have already been obtained in CUORICINO which represents 
also a vaUd test of the CUORE project feasibility . With a little R&D CUORE will reach 
a 5-year sensitivity of ~30 meV on | (rrii,) \, which could decrease to 15 meV if the full 
CUORE R&D program will be successful. 

We would like to add that CUORE is the only second generation double beta ex- 
periment based on a pilot perfectly operating prototype setup (CUORICINO) which 
consists essentially in one of the CUORE towers, thus ensuring the complete feasibility 
of its larger "brother". 

The present and future results of CUORICINO will obviously be precious to antici- 
pate the performance of CUORE. Last but not leat, CUORE is the only experiment with 
an already approved and fixed underground location. 
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During the preparation of this Proposal to which he validly contributed, Angel 
Morales passed away. We lost an outstanding physicist, a devoted teacher and a dear 
friend. 
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